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PREFACE 

The chief objective of the University of Wyoming ERTS-1 
investigation (contract NAS 5-21799) was to develop exper- 
tise in the application of ERTS data in resource studies and 
to evaluate the ERTS-1 data in these applications. The study 
began prior to the ERTS-1 launch in July of 1972 and will 
continue through October, 1975 under contract modification 7 
(10 January, 1975). The original objectives were achieved by 
December, 1974 and results of those studies are summarized in 
this report. These studies emphasized geologic applications 
of the ERTS data but also included hydrologic, agricultural, 
land-use, and ecologic evaluations of the ERTS-1 data. 

Results of these studies were mostly positive. The 
regional , repetitive coverage proved invaluable in several 
applications and the spectral information provided by the 
ERTS multispectral scanner allowed investigators to make many 
subtle distinctions regarding lithologic units, possible min- 
eral deposits, soil moisture, and vegetation condition. 

Major limitations of the ERTS data are its lack of stereo- 
scopic geometry and low resolution. 
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INTRODUCTION 

The Wyoming ERTS investigation, originally designated as 
proposal 364, was combined with EREP investigation (No. 434) 
and was made part of a single coordinated study involving ac- 
quisition and interpretation of data resulting from both the 
ERTS and EREP experiments. The test site was the state of 
Wyoming (figure 1). Initial funding requests were reduced 
substantially as a result of the coordination and also due 
to modification of original objectives. Basically, geologic 
and vegetation studies were retained and certain environ- 
mental and wildlife studies were eliminated. In keeping 
with NASA requests a sincere effort was made to make the 
results of studies available early in the form of special 
reports. Summaries of these special reports are appended as 
a part of this final report. 

LABORATORY AND OPERATIONAL RESEARCH 

In addition to the research projects, other functions of 
the University of Wyoming Remote Sensing Laboratory include 
maintainance of a library of remotely sensed data for the 
State of Wyoming, and demonstrations to other potential 
users the advantages or disadvantages of remote sensing for 
their applied research. The following results were achieved! 
1. A data library has been established which comprises 
NASA ERTS-1 imagery, Skylab photography and other 
EREP data, and aerial photography obtained from 
the Air Force, Army, Soil Conservation Service, 
and private firms. The depository also has a cata- 
log of existing maps on such subjects as land use, 
vegetation, soils, etc. where available. Main- 
tenance of this data file requires about twenty 
hours per week and is a University contribution. 

The data and the technical assistance available 
at the Remote Sensing Laboratory is used weekly by 
the public. Thus, it gives the NASA program wide 
circulation to individuals who may see for them- 
selves the utility of the data. 
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2. Meetings with individuals and groups in State and 
Federal agencies have been held to discuss the 
the possible utility of remotely sensed data for 
agency use . A formal meeting was held in September , 

1972 in the state capitoi with Dr. Short of Goddard 
in attendance and the meeting resulted in joint 
studies and use by individuals of several agencies, 
as follows: 

a. The Department of Economic Planning and Develop- 
ment commissioned a pilot study (emphasis on 
land-use maps) for the Powder River Basin of 
Wyoming (appendix A) . This study led to a joint 
proposal submitted as an ERTS-B proposal for 

a more detailed demonstration project in this 
area and elsewhere in the state. The program 
was not approved by NASA, but the United States 
Geological Survey has recently (September, 1974) 
indicated a desire to fund parts of the program 
in the Powder River Basin of Wyoming and Mon- 
tana o 

b. Wyoming Game and Fish Commission researchers 
have made use of ERTS color infrared images 
and aircraft infrared images for vegetation 
mapping related to game control programs. 

c. The State Highway Department has indicated a 
willingness to allow use of their aircraft 
for programs in remote sensing. These pro- 
grams are pending, but will involve use of the 
highway department Zeiss aerial camera plus 
the three-camera Hasselblad unit designed by 
the remote sensing laboratory for various 
geologic and earth resource studies. 

d. The Soil Conservation Service has cooperated » 

with the Remote Sensing Laboratory by fur- 
nishing new and unpublished soils maps of 
various areas in the state to see if ERTS, 


1 I J 1 1 i 


aircraft, and EREP images can. be used to generate 
soils maps. The maps prepared to date (appendix 
A, p. 72 ) are satisfactory general soils maps, 
but the units mapped do not relate directly to 
the detailed soils classification schemes in 
current use. The Soil Conservation Service 
has also supplied the laboratory with a vast 
collection of enlarged black and white air 
photographs that are extremely useful in deter- 
mining such factors as changing land use from 
the 1940 p s to the present. 

e. The Wyoming Geological Survey has cooperated in 
the remote sensing studies by supplying geologic 
data and information on current and proposed 
mineral developments. The Geological Survey 
hired the first individual trained by the Remote 
Sensing Laboratory and is currently using the 
data in the Remote Sensing Laboratory on an 
operational basis. Examples of maps prepared 
by the State Geological Survey in Cooperation 
with the Remote Sensing Laboratory are in 
appendices B and C. 

f. The United States Geological Survey uses 
remotely sensed data in many of its programs , 
and maintains a file at the Denver Federal 
Center that contains much of the available 
imagery for the Rocky Mountain area. But 
geologists of the Survey have also made some 
use of Wyoming imagery files; and the U.S.G.S. 
is currently cooperating with the Wyoming 
Remote Sensing Laboratory in two projects: 

l) Powder River Basin landforms and land-use 
studies, and 2) in preparation of the new 
state geologic map. A portion of the geologic 
mapping done with ERTS imagery has been used 
to compile the state map where mapping is 
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out-dated or incomplete. Quaternary rocks and 
rocks of Precambrian age are generally poorly 
mapped and the U.S.G.S. has used maps of sand- 
dune areas prepared by the remote sensing lab- 
oratory and has supported mapping of other 
deposits of Quaternary age using ERTS imagery 
and Skylab photographs. 

g. Research scientists at the University are aware 
of remote sensing and use is being made of 
imagery by engineers, agricultural scientists, 
and plant scientists. Perhaps the major use of 
remote sensing has been in the preparation of 
environmental impact statements . The remotely- 
sensed data serves as a tool by which researchers 
in many disciplines can obtain physical and 
environmental data in areas of proposed devel- 
opment. This information can then be applied 
to planning and impact estimation. The data 
available prior to development also represents 
a set of baseline data for change assessments. 

The development of these programs in this relatively 
short period of time show that it is possible to move directly 
from a successful demonstration program to an operational in- 
formation system. Such a system must have broad application 
and be readily accessible to a large number of users* 

DISSEMINATION OF DATA 
Public 

An effort has been made to acquaint the general public 
with the potential of NASA-acquired remote sensing data 
through talks to civic groups, news releases, and through an 
exhibit at the Wyoming State Fair in August, 1974. Non- 
technical presentations dealing with the ERTS-1 program include: 

1. Marrs, R. W. , Lecture on remote-sensing and ERTS 
applications; Laramie Civitans, April, 1974. 


2. Marrs, R. W., and others. Remote Sensing applied to 
problems in energy and the environment; Wyoming State 
Fair; August 19 - 25, 1974 (program previewed in 
several news releases ) . 

3. Houston, R. S. and others. University research in 
energy and environment; a program to aquaint Wyo- 
ming people with University research programs, pre- 
sented in 27 Wyoming cities . Several of these 
sessions emphasized remote-sensing and the ERTS 
program. June - August, 1974, (accompanied by 
radio broadcasts and news releases). 

4. Marrs, R. W., and Gordon, R. C., Remote-Sensing re- 
search at the University of Wyoming; Agriculture 
Extension Conference, December, 1974. 

5. Marrs, R. W., Remote-sensing research at the Univer- 
sity of Wyoming; State Legislators Workshop, Decem- 
ber, 1974, (this same presentation was given at the 
State Fire Meeting, Gillette, December 6, 1973. 

Scientific 

Presentations dealing in greater depth with remote sensing 
and the ERTS program were given before scientific and profes- 
sional audiences: 

1. Houston, R. S., and Short, N. M. , 1972, Geologic 
interpretation of ERTS imagery of western Wyoming; 
First ERTS-1 Symposium, September 29, 1972, Report 
X-650-73-10, NASA Scientific and Technical Informa- 
tion Office, Washington, D. C., p. 58 - 77. 

2. Breckenridge , R. M. , 1973, Glaciation of north- 
western Wyoming interpreted from ERTS-1, Symposium 
on Significant Results Obtained From ERTS-1, March 
5-9, 1973, Report NASA-SP-327, NASA Scientific and 
Technical Information Office, Washington, D. C., p. 

363 - 370. 

3. Breckenridge, R. M. , Marrs, R. W., and Murphy, D. J., 
1973, Remote sensing applied to land-use studies in 
Wyoming; Symposium on Significant Results from ERTS-1 
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March 5-9, 1973, Report NASA SP-327, NASA Scientific 
and Technical Information Office, Washington, D. C., 
p. 981 - 990. 

4. Parker, R. B., 1973, Regional Fractures in the Wind 
River Mountains Precambrian from ERTS-1 Imagery; 26th 
Annual Meeting of the Geological Society of America, 

Rocky Mountain Section, (abs.), p. 1520. 

5. Houston, R. S., 1973, Application of remote sensing 
to assessment of Wyoming* s natural resources, Wyo- 
ming Geological Society Meeting, May 11, 1973, Lara- 
mie, Wyoming. 

6. Houston, R. S. and others, 1973; Application of the 
ERTS system to the study of Wyoming resources - with 
emphasis on the use of basic data products; Third 
ERTS-1 Symposium, December 10 - 14, Report NASA-SP- 
351, NASA Scientific and Technical Information Office, 
Washington, D. C., p. 595 - 620. 

7. Houston, R. S. and Marrs, R. W. , 1974, Some prelimin- 
ary results of the use of S— 190A and S— 190B Skylab 
sensors for photogeologic studies; 55th Annual Meet- 
ing of the American Geophysical Union, Washington, 

D. C., (abs.), EOS p. 344. 

8. Marrs, R. W., 1974, Remote sensing in mineral explora- 
tion; 17th Mineral Symposium of the Eastern Wyoming 
Subsection of the Society of Mining Engineers/AIME , 

May 11, 1973, 15 p. 

9. Marrs, R. W., 1974, Summary of remote-sensing research 
in Wyoming; Soil Survey Work Planning Conference, 
Laramie, Wyoming, May 21, 1974. 

10. Marrs, R. W., 1974, Interpretation Techniques in 

remote sensing; Wyoming Geological Society Meeting, 

May 19, 1974, Laramie, Wyoming. 

Houston, R. S. and Marrs, R. W., Examples of the use 
of standard data products of the NASA, ERTS, Skylab, 
and aircraft programs in the study of Wyoming surface 
resources, American Association of Petroleum Geologists, 
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Rooky Mountain Section Meeting, June, 1974, Casper, 
Wyoming . 

12. Marrs, R. W. , 1974, Interpretation techniques in 
remote sensing; American Association of Petroleum 
Geologists, Rocky Mountain Section, June, 1974, 

Casper, Wyoming. 

13, Houston, R. S., and Marrs, R. W., Use of Skylab 
S-190A and S-190B photographs in Wyoming earth 
resource studies; 20th Annual Meeting of the Ameri- 
can Astronautical Society, August 20 - 22, Los Angeles, 
California. 

OBJECTIVES 

The primary objectives of the 434 project were to set up 
a laboratory of remote sensing and to use the laboratory to 
interpret data obtained from both ERTS and EREP experiments, 
j The objectives as modified, by the negotiations that resulted 

in combination of the ERTS and EREP studies are as follows: 

\ Monitoring 1. Collect field data to document micrometeoro- 

] logical parameters and spectral signatures, 

! for selected ERTS, EREP, and aircraft passes. 

I 

j Geology 2. Evaluate sensor data for application to 

j regional mapping, tectonic analyses, strati- 

f graphic studies, and mineral exploration, 

j This work will include an effort to dis- 

tinguish major rock units, study of region- 
al fracture systems, and attempts to dis- 
l tinguish mineralized zones from unmineral- 

| ized zones . 

I Geothermal 3. Utilize the thermal channel of the S-192 

imager to determine to what extent hot 
springs and other active areas can be iden- 
tified and mapped from spacecraft. (ERTS-1 
did not have a thermal channel, so thermal 
studies have been confined to the study of 
snow scenes and vegetation anomalies as 
related to thermal areas ) . 
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Agriculture- 4 . 
Forestry 


Environmental 5 . 
Quality 


Ecology 6. 


Determine the spectral reflectance char- 
acteristics of various forest and grass- 
land communities in Wyoming, map major 
vegetation types and attempt to identify 
and map seasonal changes in vegetation 
communities ■ This would include an evalu- 
ation of the usefulness of spacecraft 
imagery for determining the effects of 
grazing on plant communities, and evalu- 
ate the EREP data as a tool for mapping 
indicator plants and toxic soils. 

Prepare a base map showing current envir- 
onmental conditions in areas containing 
potential coal strip mines and the fly 
ash distribution from planned coal-fired 
power plants, and map density and extent 
of emission plumes using spacecraft imagery 0 
Attempt to relate the seasonal distribu- 
tion of wildlife to seasonal changes in 
plant communities and ground temperature 
as determined for EREP imagery. 


PERSONNEL . 

The original ERTS proposal was a cooperative proposal 
with Dr. N. Short of Goddard as principal investigator. The 
modified ERTS-EREP proposal placed the responsibility of the 
combined investigation at the University of Wyoming and R. S. 
Houston as principal investigator. Dr. Short, however, remained 
as co-investigator and contributed significantly to the study. 


Individuals who have contributed to the ERTS study are, as 


follows : 

Name and 

Title Affiliation 


Contribution to 
ERTS research 


R. Houston, 

P. I. - 
Geology 


Univ. of Wyo. 


Administration plus geo- 
logic mapping studies 
and mineral exploration 
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Name and 

Title Affiliation 

R. Marrs, Univ. of Wyo. 

Co-I . - 

remote sensing 
geology 


Contribution to 
ERTS research 

Administration and geologic 
studies plus land-use and 
monitoring 


R. Parker, Univ. of Wyo. Lineation and field studies 

Co-I. - 
structure 


D. Blackstone, Univ. of Wyo. 

Co-I. - 
structure 


R. Breckenridge , Univ. of Wyo. Glacial geology and envir- 

onmental geology 

Co-I. - 
geomorphology 


Regional structure and 
field studies 


E . Decker , 

Co-I. - 
heat flow 


L. Borgman, 

Co-I. - 
statistics 


D. Knight, 

Co-I. - 
botany 


Univ. of Wyo. 


Geothermal and computer 
processing 


Univ. of Wyo. 


Statistical methods of 
data processing 


Univ. of Wyo. Vegetation mapping 


H. Fisser, Univ. of Wyo. Range management 

Co-I. - 

plant science 

D. Lane, Wyo. Geol. Petroleum geology 

Survey 

Consultant - 
petroleum geology 
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Name and 
Title 

G. Glass, 

Consultant - 
coal geology 


B. Mears, 

Consultant - 
geomorphology 


R, Surdam, 

Consultant - 
Geochemistry 


P. Huntoon, 

Consultant - 
hydrology 


D. Love, 
Consultant - 

geology 


K. Diem, 

Consultant - 

wildlife 

management 

D. Strickland, 

Consultant - 

wildlife 

management 


Affiliation 

Wyo • Geol . 
Survey 

Univ . of Wyo . 

Univ. of Wyo. 

Univ. of Wyo. 

U • S • Geol . 
Survey 

Univ. of Wyo. 
Univ. of Wyo. 


Contribution to 
ERTS research 

Coal geology 


Regional geomorphology 


E vapor ite lakes 


Groundwater prospecting 


Geologic mapping 


Wildlife management 


Wildlife management 


R. Gordon, Univ. of Wyo. 

Research Asst . - 
plant science 


K. Kolm, 

Research Asst, 
geomorphology 


Univ. of Wyo. 


Range management 


Aeolian deposits and soils 


\ 
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Name and Contribution to 

Title Affiliation ERTS research 

B. Tomes, Univ. of Wyo. Structure and hydrology 

Research Asst. - 
geology 

J. Earle, Univ. of Wyo. Tectonics 

Research Asst . - 
structure 


M. Evans, Univ. of Wyo. Vegetation mapping 

Research Asst. - 
botany 


F. Redfem, Univ. of Wyo. Vegetation spectra 

Research Asst. - 
botany 


W. Lauer, Univ. of Wyo. Photo processing 

Technical Asst. - 
photography 


S. Baskett, Univ. of Wyo. Technical assistance 

Technical Asst . - 
general 


S. Agard, Univ. of Wyo. General geology 

Work/Study - 
Research Asst. 


N. Froman, Univ. of Wyo. 

Work/Study - 
Research Asst, 


General Geology and 
Computer Applications 


K. McCarty, 
Clerical Asst, 


Univ. of Wyo 


Film librarian 


Name and 
Title 


Affiliation 


Contribution to 
ERTS research 


L. Apperson Univ. of Wyo. 

Work/Study - 
Clerical Asst. 


K. Meister, Univ. of Wyo. 

Clerical Asst. 


Clerical 


Clerical assistance 


RESULTS OF THE INTERPRETATION OF ERTS IMAGERY 
Geology 

All geologic tasks were accomplished and some were extended 
to operational procedures that are now applied to other programs 

Geologic Mapping i Neither ERTS or Skylab images 
provide sufficient information for detailed geologic mapping 
because of low resolution. ERTS images are also limited in 
their utility by the lack of stereo-coverage. The ERTS syn- 
optic view is unparalleled however, for regional studies, and 
in areas where vegetation is relatively sparse it has been . 
possible to make reconnaissance geologic maps by establishing 
key lithologic units that can be recognized from space (the 
remote sensing unit) and to use these for mapping purposes. 

ERTS images can be best utilized where dips of formations a 
low and outcrops are correspondingly wide, or where litho- 
logic subdivisions can be recognized by contrasts among forma- 
tion groups. General mapping capability was increased by use 
of enhancement techniques such as color-additive viewing of 
selected ERTS bands, isodensity slicing, snow scenes, seasonal 
vegetation, pseudo-stereo viewing, edge enhancement and image 

masking . 

As part pf our proposal we suggested that it should be 
possible to distinguish major lithologic types in rocks of 
Precambrian age by use of ERTS imagery. Specifically, we 
hoped to distinguish terrane of mafic metamorphic rocks 
(greenstone belts) from that comprising mostly felsic rocks 
(granite-gneiss). This indeed proves to be a valuable 
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application of ERTS imagery. Details of the study are re- 
ported in Appendices D and E. 

Mapping of volcanic te.rrane using ERTS imagery produced 
mixed results. Major lithologic units, such as volcanic, 
sedimentary, and crystalline rocks in the Absaroka area of 
Wyoming, are distinguishable but most of the pyroclastic 
rocks have similar reflectance characteristics . Distinctions 
between the various pyroclastic units could only be made locally. 
An important contribution was the utilization of ERTS color- 
composites to identify intrusive bodies within the Absaroka 
volcanics . Thirty-four of the forty -nine known centers were 
identified using ERTS imagery. Recognition of structural 
features such as calderas, major faults and flow structures 
in the rhyolite plateau area of Yellowstone Park was also 
a very useful aspect of the ERTS study. Details of this 
study are in Appendix F. 

General geologic mapping has been undertaken in a number 
of areas in Wyoming, and the success of such mapping is pri- 
marily related to the skill of the interpreter and to his 
experience in a given area. The lithologic section in Wyo- 
ming comprises crystalline rocks of Precambrian age, dis- 
tinctively colored sedimentary formations of Paleozoic and 
Early Mesozoic age, relatively drab but vegetatively expressed 
sandstones and shales of Late Mesozoic age, and light colored 
sedimentary and volcano-sedimentary rocks of Tertiary age. 
Spectral reflectance measurements were ir.a.de at the Goddard 
Space Flight Center for a representitive group of samples from 
a sedimentary section in Wyoming. Field spectral readings 
were also taken at several type localities. A general review 
of the literature on spectral distinction of rock types and 
comparison of the spectral signatures for Wyoming rocks 
suggests that unique identifications cannot be made by use 
of multispectral data except in local areas and under con- 
trolled conditions (i. e., a general knowledge of the litho- 
logic section and the field conditions). Geologic mapping 
from space is best accomplished by recognizing distinctive 
units on images, identifying these units on the ground, and 


using them as guides to mapping. The ERTS color infrared 
composites are the most useful image product for most geologic 
mapping. The experienced photo-geologist can use these as 
well, if not better, than color photographs of the same reso- 
lution (Appendix E, p.lOl). 

Examples of regional geologic mapping using ERTS are in 
Appendices E and G. This mapping has been checked and supple- 
mented by use of color and color infrared aircraft photo- 
graphs and, in addition, photogeo logic maps have been prepared 
with aircraft products in two key areas (Appendix E and H) . 

Topical Mapping : The regional or synoptic view of ERTS 

has made it possible to undertake a number of topical mapping 
projects. These include mapping of glacial features and 
sedimentary deposits of glacial origin, mapping of geomorphic 
landforms and Quaternary deposits, mapping of dune fields, 
mapping of playa deposits (modern and fossil), and study of 
regional facies changes. 

The detailed work (Appendix I) with glacial features and 
sedimentary deposits of glacial origin in northwest Wyoming 
shows that the regional extent of glaciation can be determined 
through mapping of such features as flow directions, large 
moraines, outwash fans, terraces, cirques, large ••grooves,” 

»'U” - shaped valleys, and glacial lakes. In some areas, 
such as- the Yellowstone plateau, remarkable detail of glacial 
structure can be seen. It was also determined that ERTS band 7 
is expecially well suited to mapping the modern ice field, 
monitoring glaciers, and mapping ablating areas (Appendix J) . 

General geologic mapping programs accomplished using 
ERTS imagery also demonstrated that large-scale landforms can 
be readily mapped with ERTS images. Pediments may be readily 
identified because the deposits on pediment surfaces are often 
characterized by vegetation communities that differ from the 
vegetation of surrounding areas. Mapping of pediments can 
help reconstruct the history of paleo-rivers and streams. 

Other major geomorphic features such as large landslides, 
surfaces, and terrace deposits can also be recognized. ERTS 
images were utilized to establish a geomorphic province 
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classification for Wyoming - a first in this area (Appen- 

* ~ J3 >} 

dix C) . 

In windy, semi-arid regions such as Wyoming, sand dunes 
cover large areas. Wyoming dune fields were partially known 
from mapping of small areas, but the regional extent of stabil- 
ized and active dunes had never been determined. The ERTS-1 
imagery provided an ideal base for mapping the large dune 
belts and many of the smaller dune fields in Wyoming (Appen- 
dix K). In addition, basic new scientific information of 
regional mapping of dune deposits using ERTS has a number of 
direct economic applications. For example, ancient dunes 
form shallow water aquifers in much of Wyoming, they affect 
land use because they support a characteristic vegetation and 
active dune fields are not stable areas. Active and stabil- 
ized dunes also yield information on long term wind patterns. 

Evaporite minerals (gypsum, anhydrite, halite, trona, 
etc.) and rocks rich in these minerals are highly reflective. 
Thus, it is relatively easy to locate and map playa lakes 
that are entirely or partially dry and that have partially-to 
fully developed salt pans. In Wyoming a number of these 
lakes have been or are being mined for magnesium and sodium 
salts so an inventory of these lakes is of economic, signifi- 
cance. These lakes are seasonal and are best mapped in late 
fall illustrating the advantage of a repetitive mapping system 
in geologic studies . 

One of the surprises of our initial ERTS study was the 
determination that fossil playa lakes could be mapped if 
exposed in a sedimentary succession. We have been able to 
map the trona-bearing playa lake Wilkins Peak member of the 
Green River Formation in the Green River and Washakie Basins . 

In one area the interpreter was familiar with the outcrop; 
in the other he was not. In both areas interpretation corres- 
ponds closely to a field map of the Wilkins Peak as recently 
mapped by the U.S. Geological Survey. 

One of the investigation goals was to evaluate ERTS 
images as an aid in mapping regional facies changes. Mapping 
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of the Wilkins Peak Member of the Green River Formation dis- 
cussed above is one example of this. The playa lake facies 
(Wilkins Peak Member) was distinguished from other lake beds 
of the Green River Formation mostly on the basis of tonal 
contrast. Topographic changes along strike of a formation 
can also be an aid in recognizing regional facies changes. 

Such changes may be recognized by using ERTS images with 
side lap or using a summer image that shows a formation by 
tone in conjunction with a snow-image that has a pseudo- 
topographic expression. An example of change in facies 
expressed topographically is the change from resistant lime- 
stones to non-resistant siltstones and shales in the Casper 
Formation on the west flank of the Laramie Mountains (ERTS-1 
image 1514-17104) . The resistant limestone is marked by a 
prominent hogback that gradually becomes less well developed 
as shale and siltstone replacing limestone northward (Pig. 2). 

Some lithologic units have characteristic vegetation 
which is a distinct aid in mapping, especially utilizing 
color infrared composites. Where these units are continuously 
exposed over large areas, vegetative change along strike may 
be related to lithologic or facies changes. For example, in 
preparation of a regional geologic map from ERTS images in 
the Arminto area (Appendix E), it was apparent that certain 
sandstones and siltstones of Cretaceous age developed hog- 
backs with a characteristic heavy sagebrush growth on the 
dip slope of the formation. Some of these sandstone units 
were mapped. A facies change from sandstone to siltstone 
and to shale was located by use of this characteristic 
vegetative control . 

Red-bed facies of the continental Tertiary of Wyoming 
can be identified on ERTS images. These units have a distinc- 
tive yellowish hue on the standard false-color infrared com- 
posites, but they exhibit somewhat higher reflectance than 
the Permo-Triassic red beds that also appear yellow on the 
false-color infrared. Mapping in the Arminto area demon- 
strated that these continental red bed facies could be dis- 
tinguished from other Tertiary sedimentary rocks. 



Figure 2. Snow scene of the Laramie Range (Image 1514- 
17104). Changes in topographic expression of 
sedimentary units along the west flank of the 
Laramie Range (arrows) correspond to facies 
changes . 
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mineral Emigration : Mineral exploration using spacecraft 

images can be either direct or indirect. Direct exploration 
might be accomplished through the photointerpretive recognition 
of some characteristics of a mineral deposit that identifies 
a specific type of deposit. Indirect exploration techniques 
might employ a space image to target features such as structure 
or rock type that may then be studied in detail by use of 
aircraft images, geophysical techniques, or surface geologic 
mapping to further isolate economic mineral deposits. This 
latter technique might also be referred to as multi-level 

sensing for mineral exploration. 

Alteration halos that develop in conjunction with deposi- 
tion of some types of minerals or through oxidation and hydra- 
tion of mineral deposits exposed at the surface present a 
much larger target for exploration than the actual deposit. 
These alteration halos may be many times the size of the 
mineral deposit and they may form a "color anomaly." Two 
types of alteration were studied in Wyoming — that associated 
with disseminated copper and molybdenum deposits in igneous 
intrusive bodies in the Absaroka Mountains, and that associated 
with uranium deposited in sedimentary rocks of Tertiary age 
in the Wyoming basins. In the Absaroka Mountains, 35 of the 
49 known intrusive centers could be recognized using color 
infrared ERTS images. These 35 intrusives were studied by 
several interpreters to determine if altered intrusives 
(mineralized) could be distinguished. Results were mixed 
some altered intrusives with pronounced alteration were recog- 
recognized by all interpreters whereas other altered intru- 
sives were variously interpreted and a number of areas were 
selected by interpreters where altered bodies were not known 
(Appendix F) . It is clear that visual recognition of altered 
areas is not statistically reliable, but it is also clear 
that areas of strong alteration can, indeed, be recognized. 
Therefore, in a given area where the presence of altered 
intrusions is suspected, recognition of a "color anomaly" 
by several interpreters may be considered a reliable indica- 
tion of an exploration target. 
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In the Wyoming sedimentary basins oxidizing groundwater 
has moved through and altered arkosic sandstones of Tertiary 
age and uranium has been concentrated at the interface between 
altered and unaltered sandstone. Some of these sandstones 
have been uplifted and exposed at the surface so that now the 
alteration can. be used as a key to exploration for uranium. 

A study area was selected in the southern part of the Powder 
River Basin of Wyoming, and several interpreters were asked 
to examine the imagery of the area to see if altered sand- 
stone could be mapped from space. This attempt met with 
limited success • The altered areas interpreted from the 
imagery were poorly defined and did not correspond to altera- 
tion zones mapped on the surface. This was essentially the 
same approach employed in the previous study of alteration 
resulting from oxidation of sulphide minerals in a desseminated 
copper deposit. In the areas of copper mineralization the 
iron oxide and hydrous iron oxides hematite and goethite are 
common alteration products, and the color anomaly is often 
strong. The color anomalies associated with the uranium 
properties are generally faint. This may account for the 
relatively poor success in using ERTS color infrared images 
to identify the color anomalies in the areas of uranium 
mineralization . 

Attempts were also made to use the ERTS-1 imagery to 
identify fossil black sandstone deposits that are present in 
sandstone of Late Cretaceous age. Color infrared composite 
images and band-7 images were used to try to identify these 
deposits because it was expected that the iron absorption 
band in the near infrared would enhance the contrasts between 
the iron-rich black sandstone and associated iron-poor 
quartz sandstone. Unfortunately, the black sand deposits 
are so small that they are below ERTS resolution limit and 
could not be satisfactorily identified by visual methods. 
Larger deposits in northern Montana and large modern deposits 
on beaches might be recognized by use of ERTS images » 
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An example of the indirect application of ERTS to mineral 
exploration is the use of ERTS images in mapping greenstone 
belts in rocks of Precambrian age. Interpretations of aerial 
photographs of the greenstone terrane identified several 
suspicious outcrops that proved to be iron formations. The 
iron-formation outcrops could not be dectected on the ERTS-1 
imagery, but the broad outcrops of greenstone contain in the 
iron formation are mappable, and can serve as a guide to 
exploration. This multi-level sensing approach is discussed 
in detail in Appendix D. 

Multi-level sensing techniques can also be used to supple- 
ment direct exploration methods. For example, verification 
of color anomalies by use of aerial photography may be less 
expensive than field checking. For example, in the rugged 
Absaroka Mountains of Wyoming an aerial survey would certainly 
be faster and less expensive t ban field checking the anomalies 

located from the ERTS-1 imagery ■ 

To date, structural studies and regional geologic mapping 
with ERTS and Skylab images has not led to the isolation of an 
economic prospect in Wyoming, but the regional structural 
studies already completed or in progress may eventually lead 
to such discoveries. Oil and mining companies are now using 
the ERTS imagery in their exploration programs. Some are 
studying the linear features as an indirect guide in mineral 
exploration. Interpretation of ERTS-1 imagery of Precambrian- 
cored uplifts in Wyoming reveals thousands of linear features 
(Appendix L) . Some are faults or fractures. Some can be 
traced into the basin areas . These may represent basement 
fractures that have been reactivated and now affect the over- 
lying sedimentary rocks . Such features may control petroleum 
migration and accumulation, but prospects located by use of 
these indirect exploration methods may take years to develop. 
Attempts were made to identify some of the large linear 
features recognized in the Wyoming basins, but the field 
studies failed to identify or locate many of the major linear 
features. Nevertheless, investigators suspect that these may 
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be faint surface expressions of basement fractures • Their 
existance should eventually be verified or disproved by 
explorationists who have access to seismic data or other 
geophysical information. 

It was obvious from a first look that ERTS images could 
be used to map closed structures such as anticlines, but, in 
Wyoming, no new closed structures were detected, and some 
features that appeared to be new closed structures proved to 
be topographic pseudostructures in flat or homoclinal sedi- 
mentary rocks (Appendix M) . stereoscopic coverage would 
undoubtedly be of considerable value in segregating these 
pseudostructures from the true closed anticlines and synclines. 

Structural Geology : Most structural studies using ERTS 

imagery have delt largely with the linear features. Lineation 
study is, indeed, a major aspect of all ERTS structural 
studies, but these types of studies should be combined with 
regional geologic mapping and used in conjunction with regional 
geophysical data to be most effective. In all of these studies 
the essential field check must be included if the results are 
to be considered valid. 

More linear features can be mapped in a day than can be 
properly evaluated in several years. We have, in effect, 
been overwhelmed by linear features. We have determined that 
linear features seen on ERTS range from joint systems to mega- 
fractures. In the Precambrian-cored uplifts, abundant linear 
features include fractures and dike systems ranging in age 
from Early Pre Cambrian to Recent. In the basins the more 

recent features are more apparent. 

All three major uplifts in Wyoming have been studied (Wind 

River, Bighorn, and Laramie Mountains). R. B. Parker (Appen- 
dix N) has determined that most of the Wind River linear 
features recognized on ERTS can also be delineated on aerial 
photos and verified by field check. Many of the larger linear 
features represent zones of dislocation, cataclases, and meta- 
somatism or alteration. The fractures may or may not be 
filled by dikes. As yet, many of the major linear features 
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of the Wind River Mountains have not been verified and are 
not shown on any published geologic maps. 

One approach to analysis of linear features detectable 
on the ERTS imagery is that of Blackstone (1972; 1973a, 1973b, 
1973c) in which the trends of the linear elements in uplifts 
are compared to trends observed in structures of the adjacent 
basins to determine possible interrelationships 0 A study of 
linear elements in Bighorn and Pryor mountains determined that, 
in the Bighorns, the trends of 52 linear elements do not 
correspond to the axial trends of 83 well-defined anticlines 
and monoclinal features of the flanking basins (Appendix 0). 
The major trend of the fold axes in the basins differs in 
azimuth from the principal trend of linear elements in the 
mountain core. This may mean that fractures in the Precam- 
brian uplift did not exert a strong control on the basin 
structure. 

Another study by Blackstone and others of 616 linear 
elements of the Laramie Mountains of Southeast Wyoming does 
show a relationship between linear elements of the uplift 
and folds in the basement (Appendix P) . Northeast-trending 
linear elements in the Laramie range correspond well with the 
dominant orientation of numerous Laramide folds in areas 
marginal to the uplift. Blackstone also mapped domains of 
linear elements within areas of different major lithologies 
and found that the different lithologies within the Precam- 
brian terrane of the Laramie Range do not influence the 
orientation of linears . 

Subsequent work comparing maps of linear features inter- 
preted from ERTS imagery with those interpreted from Sky lab 
and aerial photography indicates a very strong correlation 
between dominant trends in linear elements identified from 
imagery and the sun’s azimuth at the time of the overpass 
(Appendix S). This affect is particularly pronounced in 

high relief areas. 

This complicating factor casts some doubt on the direct 
correlations drawn by Blackstone in his earlier analyses of 
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the Bighorn and Laramie Mountains because the ERTS scenes of 
Wyoming are consistantly illuminated from the southwest — 
thus, creating a natural bias toward northeast-trending 
linear elements. Enhancement of northeast-trending linear 
features would adversely affect straightforward correlations 
with northwest-trending features (such as the fold axes of 
Bighorn Basin structures) and strengthen correlations with 
northeast-trending features (such as the structures surround- 
ing the Laramie range). Consequently, the possible inter- 
relationships between linear elements and fold axes in these 
areas should be reconsidered. If possible, the ERTS data 
should be supplemented with other data obtained under differ- 
ent illumination. A composite map of linear elements inter- 
preted from images having different illumination directions 
should help to eliminate the bias created by directional 
illumination . 

Linear elements in the basins of southwest Wyoming have 
also been examined by Houston and Short and by Blackstone 
(Appendices Q and R). Prominent sets of linears striking 
north-northeast and northwest are typical of this area. 
Blackstone has described two major linears which he refers to 
as the Farson and Red Desert linear elements. He has tenti- 
tively interpreted these linear elements as the surface trace 
of deep-seated structural elements, probably faults, which 
increase in magnitude of displacement with depth. These 
linear elements are thought to control the location of vol- 
canic fields and the geometry of major Laramide folds. Black- 
stone further states that the existence of these and other 
major subparallel linear elements, extending for tens of 
miles ^ suggests a compressional stress field over large areas, 
present at least in the early stages of regional deformation. 

Field checks have shown that some linear features in the 
Great Divide Basin of western Wyoming represent linear strips 
of stabilized sand dunes separated by blowouts. These dunes 
are marked by a distinctive vegetative expression. Other 
major linear features in this area are probably real but 
could not be identified in the field, although faults of 
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minor displacement that parallel these linears have been 
mapped in this area. 

One conclusion drawn from this initial study is that 
the wind created phenomena are far more pervasive than 
previously considered in this particular region. We thus 
have a greater appreciation of the geomorphology. 

A more comprehensive evaluation of Wind River Moun 
tain linear features has been made by Tomes and others. 
(Appendix S) using ERTS, Skylab and aircraft data. This 
comparative study shows the value of both the synoptic 
view and increased resolution in lineation study and also 
demonstrates clearly that sun azimuth and sun angle play 
a major role in determining which sets of linear features 

will be emphasized in mapping. 

A regional tectonic analysis of Wyoming is currently 

in progress (Appendix M) . This analysis will consider 
all major linear features that can be recognized using 
both ERTS and Skylab images . It will not be possible to 
field check the Wyoming linear features map, but the linear 
features will be checked against known faults and fractures. 
Field evaluations and correlations with geophysical data 

will be made in key areas » 

In summary, we find that spacecraft images can be 
used to find unmapped fractures and extend known fractures 
beyond mapped limits. Most fractures recognized in general 
studies are seen in two-dimensions and are difficult 
(if not impossible) to date. Therefore, an extensive field 
study (on a regional basis) is needed before these frac- 
tures can be geneticly interpreted or used in models for 

regional structural analyses. 

In areas where rocKs are well-exposed or strongly 
expressed by vegetation, the ERTS image serves as a 
physiographic map and displays much of the structural 
detail. The image can either be annotated to show structure 
or an analysis can be made by comparison with an overlay (Fig 
3). Analyses of the information gained from image study 



ERTS-1 mosaic of Sweetwater County, Wyoming compared 
to a structure contour map representing the top of 
the Cretaceous Dakota formation. Mosaic compiled 
by U.S. Soil Conservation Service, Regional Tech- 
nical Service Center, Portland, Oregon. Structure 
contour map courtesy of Petroleum Ownership Map 
Company, Casper, Wyoming. 
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must include supplementary use of other available data, 
but the image serves as an excellent base for tectonic map 
preparation. In Wyoming studies, the imagery has been used 
as a mapping base while providing the necessary new informa- 
tion to fill gaps in mapping. In areas that are entirely 
unmapped the ERTS-1 imagery has even greater potential as a 

mapping tool . 

Geothermal 

ERTS images are generally useful for regional studies 
of the earth's surface temperature. An attempt to use ERTS 
for geothermal studies involved examination of snow scenes 
in areas of known thermal anomalies to see if melting of 
light, moderate, or heavy snow-cover might be used as an 
indicator of thermal anomalies. The two major thermal areas 
of Wyoming, Yellowstone Park am. the Thermopolis hot springs 
areas, were used as test sites. The light snow-cover was 
so non-uniform that areas of melting could not be confidently 


defined from the ERTS-1 imagery. The Yellowstone thermal 
area was consistently open but the Thermopolis area snow 
pattern could not be directly related to areas where the 
surface is known to be anomalously warm. The Yellowstone 
hot springs and fteyser areas remained open even in moderate 
and heavy snow periods (Pig. 4) suggesting that areas of 
intense thermal activity could be detected by monitoring 
snow cover . 

These studies suggest that any areas that are consistently 
open in ERTS snow scenes should be examined as possible thermal 
areas, .and that ERTS snow scenes might therefore be used m the 
search, for intense thermal anomalies in remote areas. 
Agriculture and Forestry 

Studies in agriculture and forestry have been carried 
out by botanists, plant scientists, and geologists m a 
combined attempt to define interrelationships between rock, 
vegetation, and soil. Early in these studies it was found 
that there were essentially no vegetation maps of Wyoming, 
at any scale. Consequently emphasis was placed on vegetation 
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Figure 4. Snow scene of Yellowstone National Park (portion 
of ERTS-1 image 1213-17475-7). Numbers indicate 
thermal areas detectable as zones of melted snow 
or ice. 
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mapping. A program was outlined which would use ERTS as a 
major tool to produce vegetation maps for use in land planning 
and in preparation of environmental impact statements that are 

essential background for energy development. 

Studies undertaken by botanists include the preparation 
of vegetation maps in basin and mountain areas of southeastern 
Wyoming and crop calendar studies in irrigated farm districts 
of eastern Wyoming (Appendices T and U) . These studies led 
ultimately to the preparation of the first state-wide vegeta- 
tion maps of Wyoming (Appendix V) . 

The Laramie Basin was selected for pilot vegetation studies 

because both aircraft and ERTS data are available and because 
the initial studies could be completed and checked quickly and 
inexpensively because of the proximity to Laramie. Vegetation 
in the Laramie Basin test areas was successfully mapped by use 
of ERTS. The maps were checked against aircraft data and field 
studies. The original interpretations in this area were done 
by Mr. R. F. Redfem and the fieldcheck was completed later by 
Mr. M. A. Evans. Consequently, we believe that the map evalu- 
ation of the ERTS-1 data in this application is unbiased. 

Field checks of the Laramie Basin vegetation maps showed 
the interpreter was able to delineate six major vegetation, 
communities in the basin and eight vegetation communities in 
the mountains (Table l) . 

These studies show that ERTS imagery is adequate for 
regional vegetation mapping and that density slicing is use- 
ful for refinement of classification. The results also sug- 
gest that the best time for mapping vegetation in this area 
is early fall and late summer, and that basin and mountain 

areas should be mapped individually. 

Studies by Evans and Marrs of the Wheatland irrigated 
farm area (Appendix U) show that crop inventories can be 
successfully accomplished with ERTS provided seasonal data 
and interpretation aids such as density slicing are utilized. 

An initial attempt was also made to use the Lars Programs 
and thus accomplish vegetation mapping automatically and with 
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Table 1 : Subdivisions of Vegetation Discriminable 

from ERTS-1 Imagery 

Basin Subdivisions 


Classification No. 

Species 

Major Types 

I. 

Alkaline areas 

6 

greasewood 

II. 

Light grassland 
(Moderate alkaline) 

10 

western wheatgrass 

III. 

Grassland of light 

13 

blue gramma 


cover 



IV. 

Grassland of light 

20 

blue gramma plus 

cover with greater 
moisture 


sedges and rushes 


V., 

VI., VII. Flood plains, 
no separation in early 
summer when field check 
was made — may be pos- 
sible in fall when gray 
scale or tonal map was 
made 

20+ 

lush vegetation 

VIII . 

Water 




IX o, X. Conifers and vegeta- 
tion on Casper formation — 
further subdivisions _ in 
mountain classification 


Mountain Areas 


Classification No. 

Species 

Major Types 

I. 

Limestone outcropping 
and bare areas 

23 

species of grass 

II. 

Scattered brush within 
grassland of light cover 

13 

mountain mahogany 

Ill . 

Moderate to dense brush 

13 

big sage brush 

IV. 

Dense brush 

24 

thick sage brush 

V. 

Broadleaf vegetation 
(moist area) 

29 

quaking aspen 

VI. 

Wet grasslands 

21 

grasses, sedges, 
rushes 

VII 0 

Scattered Conifers 

32 

Limber pine plus 
zone 5 vegetation 

VIII . 

Mixed pine and aspen (can 
be separated by isodensi- 
tracing 




J 


3 


1 ,a„, 


31 


greater accuracy. A number of tapes -were processed at God- 
dard and print-outs were studied by the Wyoming team. In 
this first attempt the classification proved difficult to 
transfer to Wyoming basins and was not as satisfactory as the 
visual analyses. 

The final botanical project was preparation of four 

state vegetation maps including: 

1. Map defining the agricultural areas of Wyoming, 
three divisions — — native hay meadow, irrigated 
crops, and dryland crops. 

2. Vegetation map of Wyoming showing major river 
bottoms and irrigated croplands. 

3. Vegetation map of Wyoming showing foothills scrub, 
grasslands and desert and basin. 

4. Vegetation map of Wyoming showing alpine vegetation 
areas and timbered mountain slopes. 

These general maps were prepared by use of ERTS infrared 
mosaic of Wyoming. More detailed vegetation maps were pre- 
pared in key areas within each subdivision (Appendix V) . 

Range studies have been underway in the Baggs area of 
south-central Wyoming for the past five years . The Baggs 
test area was, therefore, chosen as a control test site for 
application of ERTS-1 data to range mapping and biomass estima- 
tions. The area supports both domestic and game animals. 

Vegetation maps were prepared for this area at a scale 
of 1:24,000 using ERTS imagery, Skylab photography, and 
aerial photographs (Appendix W) . The maps compiled from 
interpretations of the aerial photography were field checked 
and subsequently used as a control base for comparison with 
interpretations of satellite imagery. 

The area was monitored during five ERTS passes and two 
EREP passes. Field spectral measurements and clippings were 
taken to determine the reflectance characteristics and above 
ground biomass (available for forage) in 25 randomly located 
plots (5' x 4*) o These reference values were to be used as a 
guide for estimation of above ground green biomass from ERTS 
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scene brightness values, but the field reflectance measurements 
could not be reduced to real values. Consequently the biomass 
reflectance dependence was determined by plotting the biomass 
estimations directly against the ERTS scene brightness values 
for each of the ERTS data sets. The initial comparisons in- 
dicate a direct relationship between biomass and reflectance 
in areas of high biomass. In areas of sparse vegetation, 
other factors such as soils and soil moisture influence the 
scene brightness such that the biomass cannot be estimated 
directly from the scene brightness. Nevertheless, biomass 
estimations may be derived for these areas by extrapolating 
the biomass changes estimated for the more heavily vegetated 
areas into the sparsely vegetated areas. 

Hydrology 

Mapping of surface water bodies and drainage networks is 
straightforward using ERTS Band-7. Mapping techniques vary 
from visual interpretation of band-7 film positives, density 
slicing using the isodensitracer , to the generation of gray- 
level maps by using digital tapes. Small water bodies (> one 
acre) can be mapped with ERTS Band-7 images; so it is possible 
to inventory water in small impoundments and lakes. In areas 
where water is critical to development one may even detect 
small man-made ponds where small streams have been impounded 
and water used for local irrigation or stock - watering . 

1972, 1973, and 1974 were marked by quite different 
rainfall in Wyoming. 1972 was a year of normal rainfall, 
snowpack, and runoff, 1973 an unusually wet year, and 1974 
a drought year. The ERTS-1 imagery records the changes in 
vegetation and water supply that occured during this period. 
One part of a continuing ERTS program (beyond the present 
contract) is a plan to study the Powder River basin and ad- 
jacent mountain uplifts comparing snow pack in the mountain 
areas, run-off in streams, extent of irrigation, duration of 
irrigation seasons, and water impoundment levels at different 
times of the growing season. This hydrologic picture is 
essential baseline information for regional environmental 
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impact studies which must be completed prior to developing 
the mineral resources of the Powder River Basin. 

Seasonal water impoundment studies were begun in the 
Laramie Basin where many of the small and large lakes are 
dependent on perched water tables or run-off from ephemeral 
streams fed by melt water from mountain snowpack. A color- 
additive comparison of these basin lakes for the 1972 (average) 
and 1973 (wet) seasons shows the striking contrast (Fig. 5). 

The clear relationship between the water levels in these 
lakes and climate is graphically shown. Because lakes of 
this type are a direct reflection the groundwater levels, their 
aerial extent can be used as an indicator of the groundwater 
level o Monitoring such lakes may also help determine whether 
or not water tables, are lowered by excess pumping of near 
surface aquifers for irrigation and/or industrial use. 

Shallow lakes can often be distinguished from deep lakes 
in the Laramie Basin by direct interpretation of green-band 
images. Shallow lakes (less than 6* deep) usually have a 
high suspended sediment level that results from wind action 
disturbing bottom sediments. Figure 6 shows that these lakes 
can be distinguished simply by study of band-5 images that 
show this suspended load. The contrast can be further en- 
hanced by use of band-5 and 7 combinations (Fig. 7). 

Plant growth in basin lakes can also be monitored from 
ERTS . Algal plumes that normally develop in late summer in 
the basin lakes can be detected. This information is useful 
to algalogists and limnologists , and also to personnel of the 
Game and Fish Commission who are interested in the deteriora- 
tion of sport fishing that parallels the algal growth. 

Continued use of the ERTS imagery in conjunction with 
Skylab research has demonstrated a practical application of 
ERTS, Skylab, and aircraft imagery in prospecting for ground- 
water. The major aquifer of the Powder River Basin is the 
Madison Limestone that crops out in four uplift areas. Black 
Hills, Big Horn Mountains, and Hartville Uplift, that surround 
the Powder River Basin on three sides. During the Paleozoic eon 
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Figure 5. Color-additive comparison of ERTS-1 imagery obtained 
in autumn, 1972 (1047-17184-7) and in summer, 1973 
(1353-17190-7) o The striking contrast in surface 
water distribution is readily apparent on this 
duotone composite «, 


TOTG INAL PAOB 
££ POOR QUALTIY 




A red-band ERTS -1 image of the Laramie Basin 
( 1353 - 17190 - 5 ) showing contrasts in basin lakes 
related to water depth and turbidity. 
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A duo-tone image composite of the Laramie Basin 
area (ERTS-1 image 1353-17190 bands 5 and 7). 

This composite shows surface water in contrasting 
colors which can be interpreted in terms of water 
depth and turbidity. 
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this limestone was exposed to erosion and an extensive karst 
topography was developed on the upper surface of the Madison 
formation. The karst was subsequently buried by later sedi- 
mentary deposits. This paleo-karst formerly was thought to 
be an important zone of permeability and porosity for the 
collection and transmission of ground water. But, recent 
surface studies of the Madison have shown that this paleo- 
karst is sealed by relatively impermeable sediment and is not 
a good source of groundwater. Instead, a recent karst has 
developed on the presently exposed Madison surface along the 
basin flanks. Initial studies suggest that water in this 
less extensive karst system is funnelled towards major fault 
systems and carried out into the basin. ERTS, Skylab and 
aircraft images have been employed in the study of these 
fault systems and associated drainage patterns. Studies of 
the fault systems in the uplift areas and correlation of 
these with basin 1 inears should provide better background 
upon which to base decisions in the search for groundwater 
in the Madison Formation of the basin. 

Environmental Quality 

A major problem facing the people of Wyoming today is 
that of the preservation of environmental quality in the 
face of major energy developments. Active and proposed 
developments include coal and uranium strip mines, oil shale 
developments, and accompanying development of power plants, 
coal-gasification plants, coal-liquif ication plants, and new 

rail and slurry pipeline transportation facilities. 

At the time of the initial ERTS proposal submitted to 
Goddard in 1971, projections indicated limited development of 
these resources in the 1970 ’s and more extensive development 
in the 1980’s with a large portion of production supporting 
coal— fired power plants. Consequently, environmental quality 
was considered a secondary aspect of the ERTS studies when 
that proposal was written. Interest in the development of 
the coal resources of Wyoming accelerated enormously as a 
result of the energy shortages of 1973 and 1974 until 
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predictions for new strip coal mine developments in the Wyo- 
ming Powder River Basin alone range from 10 (announced or 
committed 1974-1980) to 20, and these new developments will 
range in size from mines with an annual production of 3 
million tons to 20 million tons. The largest active strip 
coal mine in the United States presently produces less than 
8 million tons per year. 

We now recognize that it is essential that we establish 
baseline environmental information for areas such as the 
Powder River Basin as soon as possible. The University of 
Wyoming remote sensing laboratory has therefore undertaken a 
series of projects designed to demonstrate the potential of 
repetitive satellite imagery for the preparation of regional 
baseline maps pertinent to environmental impact and for 
future monitoring and change detection. 

The following pilot studies have been completed: 

Land use maps have been prepared for several test areas 
in Wyoming using ERTS, Skylab, and aircraft products. These 
maps have been prepared using zoom transfer scopes; or with 
stereo data, using a Kern FG-2 plotter to transfer informa- 
tion on film positives to topographic and planimetric base 
maps. The Kern Plotter is ideal for Skylab and aircraft 
photographs is not effective for ERTS images because of lack 
of stereocoverage. 

An evaluation was also made of the LARS system through 
a cooperative effort with N.M. Short of Goddard. Print-outs 
were prepared of selected areas in Wyoming at Goddard and were 
analyzed by Wyoming investigators. This study suggested that 
some modification of programs might be required before the 
LARS programs could be used for land use analyses in Wyoming. 

Examples of general land \ise maps prepared from ERTS 
imagery are shown in Fig. 8 and Appendix A. It is clearly 
possible to distinguish major categories such as urban and 
built-up land, agricultural land, rangeland, forest land, 
water, and barren land, and, in addition, some 13 subdivisions 
can be made within these categories • Maps of this type at 
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scales of Is 250, 000 and 1:500,000 are excellent for regional 
planning. The Wyoming remote sensing laboratory is currently 
negotiating with the II S. G. S. to develop a series of 1:250,000 
scale maps for use of regional planners in the Powder River 
Basin of Wyoming and Montana. 

ERTS imagery has also proven useful for topical studies 
such as mapping of flood plains (with distinction between 
natural vegetation and irrigated flood plain), regional vege- 
tation maps (previously discussed) , rangeland mapping (pre- 
viously discussed), mapping of major soil types, slope estima- 
tion, mapping of strip mining operations, mapping of water 
impoundments and drainage networks (previously discussed), 
mapping of potential areas of mineral development, and mapping 
of small urban areas. Soil maps compiled from ERTS image 
interpretations closely match regional maps prepared using 
genetic soil classifications, but these maps do not fit 
current soil classification schemes. Still these maps can 
be used for regional planning, especially in areas where soil 
types are rather related to bedrock lithology. 

A method of slope estimation was developed using drainage 
networks as a measure of relief (Appendix A) . This method 
proved quite rapid and compared favorably with similar deter- 
minations made from regional topographic maps, but the method 
is of little value where good regional topographic maps are 
available. 

ERTS imagery proved to be an excellent tool for mapping 
of large strip mining operations and associated spoil piles . 
The imagery should prove equally valuable in monitoring vege- 
tation recovery in these mining areas . Density contouring 
was used successfully for extracting maximum detail in some 
of the mined areas (Appendix A)* 

Imagery obtained with light snow-cover on the ground 
is very effective in some land use application. Transporta- 
tion networks are very apparent on this type of imagery and 
areas of rural activity, such as farm yards are readily de- 
tected. Most urban areas in Wyoming are small and difficult 
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to detect with standard or color-infrared imagery, but these 
small urban areas were outlined clearly on ERTS snow scenes. 
Such images, if taken at the proper time, show the level of 
activity in a small urban area. The business district where 
human activity is most highly concentrated will appear darker 
than the residential areas which, in turn, are darker than 
the rural areas. The snow melts faster where the activity 
is greatest. Density contouring can be of considerable 
assistance in defining boundaries of these "activity zones 
and sometimes allows further land use subdivisions to be 
made within some zones ( Fig . 9 and 10 ) . 

Fortunately emission plumes from industrial activities 
in Wyoming are seldom large enough to be detected on ERTS 
imagery. However, a major oil well fire occurred m the 
southern Powder River Basin in the fall of 1973. The fire 
burned for several months before it was contained. The smoke 
plume from this fire can easily be seen on the ERTS-1 imagery 
Fig. 11a and lib) • In some areas, such as the Green River 
Basin, industrial emission plumes are visible on the snow- 
covered imagery (Fig. 12). Study of ERTS imagery confirms 
that fly ash is being reasonably well controlled at present 
industrial sites in Wyoming, and the 1972, 1973, and 1974 
imagery can be used as a base line for comparison of new 
plant sites anticipated for this decade. 

Further studies of this problem are being under takexi 
using the higher resolution Skylab photographs. 

Ecology 

Most of the ecology studies suggested in proposal 434 
were eliminated because reviewers believed that the proposed 
programs could not be effectively completed with low-resolu- 
tion imagery. For example, it was suggested that selenium 
rich areas might be mapped either through recognition of a 
characteristic flora or selenium-bearing rock type or both. 
The selenium areas present a severe ecological problem be- 
cause selenium converter plants that grow in these areas are 
toxic to animals. Preliminary study of Skylab and aircraft 
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• Snow scene of the Casper area, central Wyoming (ERTS-l 1101-1 / 1'^ 6 ) • 
The area outlined in black corresponds to the area of the isodensity 
contour map. 


Figure 9 
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Isodensity contour map of the Casper area. This map represents a 50X enlarge 
ment from the ERTS 70mm transparency. The original map was produced in four 
colors using a Joyce Loebl/ Tech Ops isodensitracer . 
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Figure 11A. A portion of ERTS-1 image 1515-17160-4 showing 
a smoke plume from a fire at the American 
Quasar #1 Patterson well north of Douglas, 
Wyoming. The smoke plume is 60 miles long 
and the fire is large enough to be seen as 
a bright dot on the ERTS image . 



Figure lib 


A ground-level view of the fire at the American 
Quasar #1 Patterson well. (Photo courtesy of 
Wyoming Oil & Gas Commission, Casper, Wyoming). 
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Figure 12. Emission plumes from trona mills and other 
industrial sites in Wyoming can be seen on 
this winter ERTS scene (1140-17365-7). 
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images in a test area in the Powder River Basin suggests that 
these areas can indeed be mapped with the higher resolution 
imagery (Appendix X) , but not with ERTS. 

Despite resolution limitations, ERTS images have proved 
to be useful for ecological studies. For example, a study 
of mule deer in the Medicine Bow Mountains of southeastern 
Wyoming determined that snow distribution during the winter 
controls the available winter range and, indirectly, the 
deer population. Field observations indicated certain areas 
are consistently used by mule deer as winter range. The deer 
concentrate in these areas and stay until spring when they 
can return to the high country. The amount of forage avail- 
able in these winter-ranges must be sufficient to support 
the local herd or many of the deer will die before spring. 

The amount of forage available depends both upon the condi- 
tion of the range and upon the size of the area that remains 
snow-free during the winter months. The size and distribution 
of the snow-free areas where deer concentrate were determined 
from snow distribution maps compiled using ERTS imagery. 

Images taken at various dates during the winter of 1972-1973 
were examined to determine the progression of snow cover. 

In addition, maximum snow-cover during this winter and the 
following winter of 1973-1974 was determined. Images in Band-7 
proved the best for mapping snow distribution. An outline 
of snow cover was transferred to base maps with the aid of a 
Bausch and Lomb zoom transfer scope. These maps were supple- 
mented with data on snow depth taken at several monitoring 
stations extablished throughout the area. 

The snow distribution maps (Fig. 13) indicated a maximum 
snow cover around the end of January during both winters. In 
January 30, 1973 approximately 64.1 square kilometers were 
without heavy snow cover. In 1974, at about the same time of 
year, February 12, approximately 285.5 square kilometers were 
without heavy snow cover. This was an increase of 224.5 per 
cent in available winter range in 1974. The heavy snow cover 
in 1972-1973, as can be seen from the map, forced deer from 




ORIGINAL PAGE IS 
OF POOR QUALITY 


00 



0 s 10 kiMMntMn 


Figure 13b. 


January , 1973, ( figure xo t .. •T' rpnrpsented by the January 

^iorslTS^rSSS-^i JSs o? approximately 80 * of the 
resident deer population. 



i 


J 


1 


49 


their normal winter range into the Baggott Rocks, Bennett 
Peak and Beaver Hills areas to winter* This extreme concen- 
tration of deer quickly consumed the available forage and 
resulted in an extreme winter die-off estimated at 80 percent 
of the herd. 

The ERTS photos will allow investigators to determine 
accurately the minimum critical mule deer winter range in the 
study area. Similar maps for subseguent winters will provide 
a measure of the normal fluctuation in snow cover and allow 
continue correlation with changes in the deer population. 
Continued monitoring with ERTS will allow game managers to 
act quickly and provide supplementary feed for the deer 
herds as the winter range shrinks below a critical size. 

Other Applications : 

Sand dune maps prepared from ERTS imagery are being used 
by archeologists in Wyoming as a source of information about 
dunal trends and possible archeological site locations. The 
more promising sites may then be excavated in order to find 
and interpret ancient artifact and bone matter. In addition, 
landforms , badlands, and aeolian erosional features can be 
related to specific types of prehistoric food procurement 
systems, where they were used as trap elements. 

Paleontologists note shifts in wind patterns which are 
recorded in surface or paleo landforms during the Holocene 
environments . 

Similarly, climatic cycles recorded in glacial landforms 
and in patterns of fluvial deposition or erosion can be inter- 
preted in terms of their effect on ancient man or in terms of 
paleo environments. 

USE OF AIRCRAFT IMAGERY 

As a result of the ERTS and Skylab programs plus studies 
done by the U.S.G.S. for NASA in Wyoming in recent years, some 
65 percent of the state has been covered by up-to-date aerial 
photography. The primary use of this aircraft support has 
been in checking interpretations of satellite imagery and in 
adding essential detail in these studies. But the potential 
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for long term use of aircraft imagery for a variety of studies 
and the critical environmental baseline information is even 
greater o Partially because of the potential of the aerial 
photography and satellite imagery in other applications, the 
University of Wyoming remote sensing laboratory has developed 
a color photographic laboratory to process the photographs for 
research support and for local users* In addition, the labora 
tory is now a repository for other aircraft photography ob- 
tained by State and Federal agencies. 

There is nothing really new in the use of aircraft photo- 
graphs for geologic mapping, but the techniques used today 
are vastly improved (i.e. Kern PG-2) and the color imagery is 
far superior to the standard panchromatic prints. Consequently, 
it is reasonable to assume that geologists will maXe extensive 
use of these new data and techniques as a supplement to stan- 
dard field mapping methods. We anticipate that most students 
will make an analysis of all available imagery and prepare 
photogeologic maps prior to beginning their field studies. 

This will allow each student to study all areas that are well- 
exposed and to map units readily recognized on the photograph. 
They will also be able to select key areas where units have 
to be mapped on the ground and where anomalies require careful 
field check. The student may then proceed to field studies 
and finally to a compilation of all data on a topographic 
base. Completed studies (Appendix H) suggest that geologic 
maps can often be completed more accurately and in less time 
using these techniques. 

Aerial photography is presently being used in uranium 
exploration, petroleum exploration, for geologic mapping and 
for the preparation of environmental impact statements. Most 
of the several hundred visitors to the Wyoming remote sensing 
laboratory have shown a major interest in the aircraft imagery. 

DEVELOPMENT OF COMPUTER PROCESSING TECHNIQUES 

Several approaches have been used to determine the value 
of computer processing of ERTS magnetic tapes . The most 
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straight-forward approaches have been ground interpretations 
of print-outs of LARS programs. In these tesr. > several 
Wyoming tapes were processed at Goddard. We nave also 
selected certain key mineralized areas to be processed by 
the GE Image 100 system at the EROS data center and other 
advanced computer systems, but these analyses are still 
in progress. 

A major effort has been to establish a capability 
at Wyom in g for computer processing of ERTS tapes o Four 
basic computer processing techniques have been developed 
and investigated ■under the sponsorship of the NASA ERTS 
and EREP projects. These are (l) mini-computer output 
possibilities, (2) sequented Q-mode factor pattern recogni- 
tion techniques, (3) frequency space cluster search and 
alpha mapping procedures, and (4) map feature enhancement 
with a type of multivariate ratio. Each of these topics 
will be outlined under a separate heading. 

Investigations in the four areas were carried forward 
until their feasibility was verified and preliminary appli- 
cations in geology could be made. The conclusions so far 
are tentative since most of the project time was used in 
technique development. More complete application testing 
of the procedures will be carried out in the continuing 
research of the Remote Sensing Laboratory. 

Mini-Computer Output Possibilities 

One of the major difficulties of the use of digital 
computers in geological applications is the necessity of uti- 
lizing alphanumeric map plots output from the on-line printer. 
Maps of this type typically involve north-south scaling distor- 
tion as compared with east-west scaling. In addition, site by 
site mapping from ERTS magnetic tape data of regions large 
enough to treat typical geological problems involve patching 
together many line-printer output paper strips and produce 
maps which are ten feet by ten feet typically. These 
might be reduced photographically but even so they are 
difficult and inconvenient to work with. 
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One alternative to the above output processes is the 
use of mini -computer-typewriter couplings which communi- 
cate with the large digital computer used in the bulk 
processing of the data. Such devices have the advantage 
of relatively low price (less than $15,000), greatly- 
improved resolution, avoidance of directional distortion, and 
full-time dedication to the project studies. The devices 
have the disadvantage, typically, of being relatively slow. 

The Wang 700B mini-computer coupled with a IBM 
selectric typewriter was used to test the feasibility of 
preparing high resolution maps. The Wang 700B computer 
has a 859-step program memory and 120 storage registers. 

It receives data and programs by cassette tape. Other 
input options, such as paper tape, disk storage, and punched 
cards can be used with the Wang 700B although these are 
not owned by the department. 

The IBM selectric typewriter coupled to the Wang 
700B has an up-down, left-right carriage motor capable of 
plotting to 0 o 01 inch vertically and horizontally. In 
addition, colored ribbon cassettes can be snapped in or 
out easily, thus guickly changing colors for colored map 
preparation. 

Unfortunately, the Wang 700B cannot communicate 
directly with the large campus -wide digital computer. 

A new slightly more powerful mini -computer , the Wang 
Model 2200B-1, however, has this capability. It is coupled 
to essentially the same typewriter-plotter and, in addition, 
serves as a teletype connection with the large computer. 

As a means of evaluating the usefulness of mini- 
computer output in remote sensing investigations and 
determining whether it would be desirable co upgrade 
the Wang equipment to the 2200 level to support the NASA 
investigations, various mapping feasibility studies were 
made with the Wang 700B equipment. 

A program to prepare a map was written in the Wang 
language o The map input was coded for storage on a cassette 
tape which was automatically read line-by-line as needed. 


The procedure produced an adequate grey-scale map by plotting 
asterisks at various spacings. The shades varied from 
essentially a solid color to a widely-spaced spotted pattern. 

The use of colored typewriter tapes was not as effec- 
tive because the colors available commercially were dark 
shades of brown, blue, green, and red. Of these, red, 
blue, and black presented the best map contrast. 

A reasonable sized-map (ten inches by ten inches) 
required several hours output time. Less time was needed 
if wide grey-scale spacings were permitted for many of the 
outcrop shadings. More time was required if essentially 
solid color maps were made. The time was not a big factor, 
however, since the computer was available for almost full- 
time use by the project and night runs could easily 
handle longer jobs. 

The study was judged to have demonstrated the feasi- 
bility and desirability of mini -computers in remote sensing 
output mapping and, accordingly, funds for acquiring a 
Wang model 2200B-1 were included in the project extension 
proposal submitted to NASA 0 However, the extension 
was not funded and the mini -computer investigation was 
deferred until other funding sources could be located. 
Sequential Q-Mode Factor Pattern Recognition Techniques 

The use of factor or principle component analysis in 
geological investigations is becoming more and more widely 
accepted (Imbrie and Van Andel, 1964; Imbrie and Purdy, 

1962; Griffiths and Ondrich, 1969; Klovan, 1966; Davis, 

1970) o This is a technique in pattern recognition 
(Andrews, 1972, pp. 24-32) which searches for new linear 
combinations of the variables which are more sensitive 
for classification of sites and for geologic source studies. 

Let x. . be the (i, j) element in a data matrix where 
within each column of data are the various measurements 
at a site. For the ERTS data, the column would consist 
of four values which are the measured green, red and two 
infrared reflectance values for each site. If there are 
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N sites in a region of interest, the data matrix in this case 

would be four rows by N columns. 

Let X. be the arithmetic average over all the sites in 

the i-th row of the data matrix and let <r^ be the standard 
deviation of the data values in the i-th row. The row data 
values can be normed by subtracting the row mean and dividing 
by the row standard deviation to yield a normed data matrix 
Z with element Z^j where: 

z ij = (x ij - /<T i 

The usual statistical correlation matrix, R (for ERTS data), 
for the data would be the 4x4 matrix given bys 

R = Z Z B / N 

The eigenvectors for R represent the principle components 
and the corresponding eigenvalues when divided by 4 yield 
the. fractions of the data variance "explained" by each com- 
ponent o Ordinarily only several of the four eigenvectors 
are needed to "explain" most of the variance and site-by- 
site alpha maps for the component of the site data values 
parallel to a particular principle component yield the 
overall picture of the site loadings on that component. 

The map giving loadings for the component with the largest 
eigenvalue will "stretch out" the data to its maximum range 
of variation. 

The above procedure is called R-mode factor analysis in 
geological statistics because the correlation matrix, R, 
was the basis for the eigenvectors and eigenvalues. Essen- 
tially the principle components were determined by the over- 
all average correlation of one reflectance value with the 
reflectance value for another color.. The computation is 
trivial , requiring only the storage and manipulation of the 
R matrix (4x4 for the ERTS data) . 
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Geologists have found another type of factor analysis 
more useful in many applications. This is called Q-mode 
because the Q or cosine matrix is the basis for the eigen- 
vector computations. Each column of data matrix is regarded 
as a vector with lengths 

L j = V^i X ij 

The cosine of the multidimensional angle between the data 
vectors for sites j and j* is given by: 

Qj = cos (site j, site j*) = (Z^ Xyj X ij*^ ^ 


The cosine matrix Q with elements Q j j* is a N x N matrix 
where N is the number of sites. The eigenvectors and eigen- 
values of Q form the basis for Q-mode factor analysis. 

Two points need to be madeo Sites with very similar 
reflectances in the several bands or which keep about the 
same ratio relations between reflectance components will 
give cosines close to 1.0. This is because all the site data 
vectors are normed to unit length initially (by division by 
L.) and so intensity is essentially removed. The cosine 
matrix, thus, forms a natural basis for cluster search. 

Similar sites have nearly unit cosines with one another. 

The second point is that in the context of remote sensing, 
where N, the number of sites, is extremely large, the N x N 
matrix is so large it is impossible to store it in even the 
largest computer. Hence, the advantages of Q— mode are diffi 
cult to attain except in small problems where the N x N 
cosine matrix is not too big for computer manipulations. 

Statistical research was carried out under the ERTS 
contract in the search for a sequential Q-mode factor pro- 
cedure which would allow the analysis of many sites by Q— 

Under the direction of the ERTS project 
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statistician, a master's thesis entitled "Sequential Q-mode 
Factor Analysis," by Mr, Albert J, Wolff was completed in 
May, 1974, This research established the mathematical founda- 
tions for the sequential procedure. 

The procedure groups the sites into groups of a hundred 
or so and completes the GsrlOO x 100} cosine matrix for each 
group o The eigenvalues and eigenvectors are computed for 
each matrix and processed together with the data to find a 
parsimonious set of data vectors which span the data vector 
space o The (typically) small sets of spanning vectors from 
each cosine matrix is combined to form a new "data" matrix 
which is then Q-mode factor analyzed to reduce these vectors 
to a set of overall spanning vectors. This recombination 
and summarization for larger and larger groups of sites con- 
tinues until the complete map area is covered o 

In addition to the development of the basic method for 
sequential Q-mode factor analysis, the master's thesis also 
included a study of percentage of the original data informa- 
tion retained at each recombination step. Various measures 
of the information retention were studied and compared. For- 
tunately, a relatively simple measure based on the trace of 
one of the analysis matrices was found to be a good index of 
information retention in the sequential analysis „ 

During the summer of 1974, a complete set of computer 
programs were developed and tested to achieve practical use 
of the Q-mode factor analysis . These programs are listed in 
Appendix Y. Various applications in map analysis and other 
geological problems are now underway and the programs are 
subject to continuing modification to make them more useful 
in geological studies 0 
Frequency Space Cluster Search 

As an empirical check against the factor analysis, com- 
puter subroutines were prepared for making one-, two-, three-, 
and four- dimensional frequency histograms for the various 
reflectance intensities versus each other. These were dis- 
played in alpha-numeric patterns for the two or more 
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dimensional cases. Slice patterns were displayed for the 
three- and four- dimensional histograms. 

The two— dimensional histogram, particularly that for red 
versus the second infra-red band proved to be the most useful 
in delineating groupings or clusters of sites into mapping 
units o Alpha geographic maps based on the Red-IR2 tuo -dimen - 
sional frequency histograms were prepared for an area near 
Greybull, Wyoming. The maps were compared with area photo- 
graphs and most features shown in the alpha map were identi- 
fied o It appeared that many of the features were delineated 
more precisely on the alpha map than could be outlined by 
eye on the photographic representations. However this is a 
subjective judgement which is difficult to validate numerically. 

Comparisons of the 2-D histogram with the results of R- 
mode factor analysis yielded an interesting result for test 
area used in the study. The direction of the principle 
eigenvector was aligned along a very characteristic streak 
of frequencies which dominated the histogram 0 The streak 
stretched along a line extending from the origin out along a 
fixed ratio of one band’s intensity to the other . The second 
eigenvector then represented the mext main direction of scatter 
perpendicular to the principle direction or frequency streak. 
Thus, mapping data components site-by-site on the first 
eigenvector direction is equivalent to mapping intensities 
along a fixed ratio „ Similarly, mapping data components on 
the second eigenvector direction is equivalent, at least in 
the best area, to mapping by ratio. It follows that the em- 
pirical histogram cluster method and the R-mode factor approach 
both lead to mathematical extensions of a type of ratio mapping. 

Similar results are anticipated where Q-mode factor methods 
have been similarly compared with the frequency histograms. 
However these studies are still in progress and will have to 
be reported later. 

Multivariate Ratio Techniques 

The success of preliminary applications of the above 
suggested that a study of various multidimensional versions 
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of ratio procedures deserved study. Consider for example, the 
ratio of : 

X — red reflectance intensity 

Y s infra-red reflectance intensity 

Ratio mapping could be based on values of Y/x or, alternately, 
on values of 0= arc tan (Y/x) . In the two-dimensional of the 
X and V values, this turns out to be the same as a conversion 
to polar coordinates where : 

X = r cos 6 

Y = r sin 9 

and 

r = Vx 2 + Y 2 , 0= arc tan (Y/X) 

The concept can be extended to four-dimensions by con- 
version to generalized spherical coordinates o Let U, V, X, 
and Y be the intensities of the four reflectance bands. Then 
generalized spherical transformation relates these to four 
new variables 0, V 2 > and the relations: 

U = r cos y ± 

Y =r r sin V ^ cos y 2 

X s= r sint'-L siny 2 cos 

Y s r siny i siny 2 sin 

The reverse relations are; 

7777 

r = tr + + X + x 

0 = arc tan (Y/x) 
y = arc tan (X/V cox0) 
y = arc tan (V/U cosf 2 ) 

Thus, mapping on the vector (0, Y 2 , represents a multi- 

variate generalization of ratio mapping. 
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The choice is not unique. The ordering of the reflectance 
bands smong U, V, X, and Y produce apparent differences m 6, 
y , and y . In a spherical sense these are equivalent and 
merely represent different four-dimensional choices corres- 
ponding to the selection poles, equator and longitude origin 
in three-dimensional space. There are, however, more subtle 
choices in that the listed equations are not the only general- 
ized spherical transformation possible (Hiller, 1964, p. 9 - 
12 j Shelupsky, 1962). 

A test area near Pumpkin Buttes, which is geologically 
interesting in that surface coloration indicates uranium de- 
posits, was selected for study relative to map enhancement by 
generalized spherical transformations „ The techniques of 
factor analysis and frequency histogram cluster analysis were 
also applied to this same area so that the relative advantages 
and disadvantages of the various techniques could be compared a 
This study is still in progress. 

COST BENEFIT ESTIMATES 

Accurate cost/benefit analyses cannot be made for most 
applications of the ERTS data to natural resource studies be- 
cause most of the benefits are not immediately realized., 

General geologic mapping, tectonic studies, and mineral ex- 
ploration may eventually lead to the discovery of mineral 
deposits o Likewise, sound planning for resource development 
may save millions of dollars worth of resources that might 
otherwise be wasted or minimize damage to the environment. 

Such benefits will only be realized over a period of many 
years and are largely dependent upon acceptance of the ERTS 
data as a standard tool by industry, governmental agencies, 
and the public „ Such acceptance will probably be very gradual 
because the use of the ERTS data often necessitates considerable 

change in standard procedures o 

Nevertheless, some rough cos*t benefit estimates can be 
made for general mapping projects by estimating the improve- 
ment in mapping efficiency resulting from the use of the 
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ERTS-1 imagery. Unfortunately, even these estimates consider 
only the immediate benefit resulting from increased speed 
and accuracy of the mapping program. The estimates do not 
provide a measure of the long-term benefits or those bene- 
fits resulting from the applications that could not be accom- 
plished "without ERTS. 

An estimate was made of the cost effectiveness of ERTS 
in the regional land-use mapping program conducted in the 
Rawlins area (Fig. 8). Comparing the initial cost of the 
ERTS imagery to the cost of aerial coverage of the same area 
a 10-fold (or greater) savings was realized. This estimate 
does not consider the initial cost of the ERTS launch and 
its subsequent operation nor does it take into the account 
the actual cost of flying and aircraft mission. In terms of 
time saved in compiling the 1 : 250, 000-scale land-use map, a 
4-fold reduction in interpretation time was estimated using 
the ERTS imagery in lieu of high-altitude aerial photography. 

Some detail that might have been derived from interpretation 
of the aerial photography was not available from ERTS, but 
this loss is somewhat compensated by the additional under- 
standing of regional relationships gained from the synoptic 
viewo 

A similar time/cost saving can be achieved by using 
ERTS imagery for regional geologic mapping or land-form 
studies. Typically, a compilation /mapping program covering 
an area of some 7000 square miles (the area of a 2° AMS map) 
would require 9 to 15 months and an investment of $ 20,000 
to 30, 000 o Using the ERTS imagery as a mapping base and 
using published information and field checks to obtain 
auxilliary control data, the same area might be mapped in 2-3 
months and cost less than $ 10,000. The 3-fold savings in 
cost is only a partial estimate of the benefit because addi- 
tional benefit is provided by the savings in time and man- 
power. In some applications the time factor may be crucial. 
Certainly, the timelyness of ERTS-assisted mapping programs 
enhances the value of the products and adds to the indirect 
benefits by making the information available earlier. 
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SUMMARY 

The ERTS-1 satellite data have proven extremely effective 
in a number of resource applications • Regional repetitive 
coverage is a major asset in general geologic mapping, vege- 
tation mapping, land-use studies, and tectonic analyses. The 
spectral information provided by ERTS is particularly valuable 
in vegetation studies, mineral exploration, and other 
applications involving surface geochemistry or soil moisture. 
The ERTS data were of very limited value in geothermal 
studies , because only the large and relatively intense 
thermal anomalies are detectable with ERTS, and even this 
capability is dependent upon having the proper surface 
conditions during the overpass. 

The most severe limitations of the ERTS-1 data were 
found to be the lach of stereoscopic coverage (an in- 
herant limitation of the Scanner) and the relatively low 
ground— resolution of the system. Comparisons of ERTS 
image interpretations and interpretations of Skylab 
photographic data show that the additional resolution 
provided by the SKylab data is valuable in many resource 
studies. Yet, the greater resolution is of little ad- 
vantage if complete and repetitive coverage must be 
sacrificed to achieve the increase m resolution. 

The computer compatibility of the ERTS multispectral 
scanner data lends greater flexibility to the ERTS-1 
data by providing a means of performing complex spectral 
or spatial pattern recognition or by allowing rapid 
classification of some features. However, most resource 
studies require that the interpreter emphasize certain 
aspects of the data while ignoring other features which 
he considers "noise." The interpreter relies upon complex 
reasoning, correlation, and experience to make these de- 
cisions and the result is somewhat subjective. It is 
difficult or impossible to program these complex decision 
criteria into a computer , so in most resource applications 
the roll of the computer is limited to providing the data 
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user with a more interpretable product or to help him 
in making rapid comparisons „ 

APPENDICES 

(see following pages) 
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Appendix A : Summary of Applications of ERTS-1 data to 

regional land-use studies in Wyoming. 


from a report compiled for the Wyoming Department of 
Economic Planning and Development: 

Breckenridge, R. Mo, Marrs, R. W., and Murphy, 

D. J., 1973, The application of remote-sensing data to 
planning in Wyoming — A demonstration project. 

Results of this study were also reported at the March, 

1973 Symposium on Significant Results Obtained from 
ERTS-1: 

Breckenridge, R. M. , Marrs, R. W., and Murphy, 

D. u . , 1973, Remote sensing applied to land-use studies 
in Wyoming, in Proceedings of the Symposium on Signi- 
ficant Results Obtained from ERTS-1; NASA/Scientific 
and Technical Information Office, Washington, D. C., 
v. 1, p. 981-989. 

Land-use considerations have become critically important 
in Wyoming where greatly increased development of natural 
resources is anticipated. Of prime concern at the present is 
the development of coal resources in the Green River and Powder 
Basins as fuel for huge generators. The direct impact of 
mining, the diversion of water, and the effluent from coal- 
fired electrical generators will be considerable, but the asso- 
ciated increase in population and the necessary urbanization 
will have far-reaching effects on the environment in Wyoming. 

The construction and operation of the power plants and supporting 
mines will result in a substantial increase in the state *s 
population (perhaps as much as 50%) . These population increases 
will, in turn, affect almost every aspect of land-use. 

Governmental agencies have recognized the need for 
thorough land-use planning in order to properly control develop- 
ment and prevent costly or unnecessary damage to the land and 
its inhabitants but there is very little land-use data 
available for use as a basis for such planning. As is commonly 
the case, the problem requires an immediate solution. It 
is totally impractical to begin a program for gathering 
detailed land-use data by conventional means because 1) such 
a program would be too expensive, 2) insufficient time is 
available, 3) much of the data would be obsolete by the time 
the program is complete. Consequently, planners must find 
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a rapid and efficient means of obtaining the necessary land- 
use information and a -way to continually up-date that informa- 
tion. 

As soon as it became apparent that the ERTS satellite 
program offered considerable potential for obtaining much 
of the needed land-use data, the state planning office of the 
Department of Economic Planning and Development contracted 
the University of Wyoming Remote Sensing Labor atory for a 
pilot study to determine what land- use data might be gained 
from the ERTS-1 imagery. The pilot program was to be a quick- 
look study designed to assess the ERTS-1 imagery as a tool 
for acquisition of land-use data. The Powder River Basin was 
chosen as a general test area (Fig. A-l) and suitable sub- 
areas within the basin were selected for testing individual 
application. 

Because the pilot program was time-limited, a suite of 
tasks was selected as a representative group of ERTS appli- 
cations which might yield information necessary for land-use 
planning. Included in this group were 

1. Identification and mapping of water impoundments 

2 o Detection of potentially unstable areas 

3. Mapping of flood plains, irrigated land, and crops 

4. Discrimination of coniferous and deciduous trees 

5 o Range mapping (type and/or quality) 

6. Mapping of soil types 

7 . Slope estimation 

8. Mapping of mined areas 

9. Mineral resources mapping (actual and potential) 

10. Areas having potential for construction materials 

11. Mapping of areas of urban and rural development 

Application of standard photointerpretation, color- 
additive viewing, and/or isodensitometry to these tasks 
provides results which suggest that many of these tasks can 
be done routinely. Figure A-2 is a land-use classification 
map for the Powder River Basin test area. Constructed from 
interpretations of ERTS imagery. In addition to the broad 
land-use categories represented on this map, there is a great 
deal , 'of other significant land-use information available from 
the ERTS data. 

Water impoundments can be readily mapped from the infrared 
(band 7) ERTS imagery. Impoundments covering 5 acres or more 















Figure A-2. Land-use classification map for the Povder River 

Basin, Wyoming constructed from ERTS-1 image inter- 
pretations. 
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can be mapped reliable and most impoundments larger than 
1 acre can be located 0 However, the actual size of these 
smaller impoundments cannot be accurately estimated because 
they are represented by too few scanner resolution elements. 

Although in some parts of Wyoming, local unstable areas 
(landslides, slumps, etc.) have been identified from the 
ERTS imagery none could be identified in the Powder River 
Basin test area. Potentially unstable areas were identified 
on the basis of their regional associations. The potentially 
unstable areas in the Powder River basin mappable from ERTS 
include 1) dune fields, 2) flood plains, 3) areas of potential 
mineral development (coal, oil, uranium), and 4) areas likely 
to be underlain by clay-bearing soils. Identification of these 
areas from the ERTS imagery required integrated interpretation 
of all ERTS bands. Color-additive viewing provided some 
enhancement of the flood plains and aided in soil and litho- 
logic distinctions . 

Color-additive viewing was of considerable utility in 
mapping of natural flood plain vegetation and irrigated farm 
land (Fig„ A-3). The red and infrared bands (5 and 7) have 
the greatest utility in distinguishing contrasts within the 
vegetation. Interpretation of the ERTS-1 imagery provided 
adequate detail for mapping broad-scale agricultural patterns. 
Some detail features were recognized in the agricultural 
pattern 0 These sometimes allow crop identification if supple- 
mented with other image interpretations showing seasonal 
variations 0 

Maps of distribution of conifers and deciduous vegetation 
were compiled from ERTS-1 interpretations for the Bighorn 
Mountains (Fig. A-4). These maps are necessarily generalized 
due to the very small scale of the ERTS images, but they 
yield a fairly reliable estimate of the distribution of these 
major forest subdivisions. 

Rangeland is readily distinguishable from other agri- 
cultural lands and can be subdivided on the basis of tone 
and texture (Fig. A-5). The resulting map can be correlated 
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streams in the Powder River 
ERTS-1 imagery. 







Figure A-4. Map of forest and rangeland in a portion 
of the Bighorn Mountains (interpreted 
from ERTS-1 image 1013-17291). 
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Figure A-5. Range map of part of the powder River Basin test site* 
(interpreted from ERTS-1 images 1030-17233 and 1030- 
17235) . 
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to standing biomass at the time of the satellite overpass, but 
range type and quality (a time -dependent factor) cannot be 
established "without control data and complementary interpre- 
tations from repetitive coverage. 

Broad soil classifications can be interpreted from the 
ERTS imagery on the basis of tone, texture, and regional associ- 
ation (Fig. A-6). The resulting maps show excellent correspon- 
dence with available soil maps (Dunnewald, 1957 5 Thorpe, 1939) 
which classify soils on the basis of composition (which 
relates to parent material) 0 Very little correlation is 
evident between the soil interpretations derived from ERTS 
and the more detailed soil maps based on the newer Seventh 
Approximation classification, but similarities may exist 
between the larger soil groups and the ERTS interpretation. 

Slope estimations were based on stream pattern analysis. 

The frequency and patterns of drainages were used as criteria 
for subdividing the Powder River Basin into zones of equivalent 
slope (Fig. A-7). The actual slope values for the various 
areas were then determined from the topographic maps. This 
application was particularly successful in the Powder River 
Basin which has a very intricate system of small drainages 
that reflect even subtle changes in regional slope. 

Abandoned, producing, and developing open-cut coal mines 
are visible on the ERTS imagery, but small open-cut uranium, 
bentonite, and gypsum mines and surface evidence of under- 
ground mines in the Powder River Basin could not be seen on 
the ERTS imagery. However the large uranium mines in the Gas 
Hills district of central Wyoming can be seen on the imagery. 

A map of known and potential coal resources was constructed 
from the combined mine locations and geologic interpretations 
of the Powder River Basin imagery. This generalized map 
(Fig. A-8 ) is useful for defining areas of probable or possible 
coal mine development for the proposed electrical power 
installations or other energy needs. Maps of potential areas 
for the development of other mineral resources were compiled 
with the favorable locations defined on the basis of geologic 



Figure A-6. Soil map of the west-central portion of the Powder 
River Basin. (interpreted from ERTS-1 images 1030- 
17233 and 1030-17235), 
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Figure A-8. Coal reserves of the central Powder River Basin. 

Coal outcrop patterns were interpreted from 
ERTS-1 data. 
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interpretation (Fig. A-9). For example, dunes and terraces 
■were designated as potential sources of sand and gravel, 
areas of limestone outcrop were mapped as possible sources 
of lime, and areas of structurally favorable geologic patterns 
i were designated as areas of potential petroleum resources 

(although virtually all structures visible on the surface 
have already been tested) . 

Perhaps the most difficult task attempted was the mapping 
of urban areas. Cities and towns of Wyoming are relatively 
small (population 60,000 or less) 0 Early attempts to outline 
the areas of urban development by standard photo interpretation 
techniques were unsuccessful. Although the urban areas could 
generally be located on the imagery, there was often insufficient 
contrast with the surrounding rural development to allow 
definition of the urban area. Color-additive viewing was 
useful in some instances where subtle contrasts on two or 
more bands could be enhanced by composite display , but it 
was soon discovered that the ability to discern urban areas 
was largely dependent on the overall contrast with the sur- 
rounding area, and that the urban areas were generally too 
small to exhibit a characteristic texture or pattern. Even 
under magnification the urban areas generally appear as fuzzy , 
grey areas with subtle variations in internal tone and texture □ 
The patterns are too intricate and tonal variations too subtle 
for visual interpretation,. 

This problem was circumvented by employing the isodensi- 
tracer as a tool for enlarging the images for the urban areas 
(up to 50x) and enhancing subtle density contrasts. Inter- 
pretation of the resulting isodensitj' contour maps usually 
allows separation of the urban areas from the surrounding 
rangeland and natural vegetation but still provides little 
distinction between urban and agricultural land. A more 
powerful technique was devised to segregate urban areas from 
agricultural land. In this technique, four color isodensi - 
tracings are prepared for each of the multispectral images and 
l the isodensitracings are then composited to form a color 

J composite map. V 

a 
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Figure A-9. Areas of potential! mineral resource development along 
the west fland of the Bighorn Mountains. Areas were 
identified by interpretation of ERTS-1 imagery. 
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Both isodensitracing techniques reveal patterns within 
the urban areas that can often be related to type or intensity 
of development. In some instances the industrial and/or 
business districts can be segregated from residential areas. 

New residential areas are commonly distinct from older resi- 
dential areas which have mature trees. Parks or recreational 
areas are sometimes identifiable. 

Interpretations of the isodensity contour maps from 
winter ERTS images compare very favorably with interpretations 
of comparable aircraft photography and provide considerable 
data for updating most city maps (Fig. A-10) . 

This study demonstrates, that even though photointerpretative 
techniques using ERTS data cannot provide detail that many 
land-use applications require, it is capable of providing a 
great deal of broad-scale, basic land-use data on a continuing 
basis. A classification scheme such as that proposed by 
Anderson and others (1972) may serve to convert the image data 
to useful land-use maps. This 'information can be derived 
from the ERTS data far more quickly and efficiently than by 
conventional techniques. Consequently, the ERTS program can 
supply information for land-use palnning as long as there is 
a need for constructing or revising the broad-scale land-use 
data base. 
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Appendix B s Application of remote-sensing data to geologic 
problems encountered by the Wyoming Geological 
Survey. 

summarized from a Wyoming Geological Survey Reports 

Breckenridge, R. M. , 1973, Geology of the proposed 

Deer Creek reservoir site, Wyoming Geol. Survey, 

Laramie , Wyoming . 

The proposed dam site is located in the SE% of Sec. 11, 

T. 31 N., R. 77 W. An alternative site is located in the NE% 
of Sec. 13. 

The area is located at the northern end of the Laramie 
Range which was uplifted during the late Cretaceous. Sub- 
sequent erosion has removed thousands of feet of Paleozoic 
and Mesozoic rocks thus exposing the Precambrian core. The 
northern margin of the uplift is extensively folded and faulted 
as evidenced by the disturbance of the flanking sediments. 

No detailed maps of the area are published. A part of 
the area was mapped by Sears (1948) before topographic control 
was available and does not agree with the present-day land 
survey . 

A generalized reconnaissance map of the proposed site 
was compiled using color aerial photographs of the area 
(NASA FLT 213, Sept., 1972) from the University of Wyoming, 

Remote Sensing Laboratory. The immediate area of the dam site 
was field checked. The proposed location is in a steep canyon 
cut into Precambrian rocks. Both abutments of the structure 
would be in bedrock of granite with some mafic dikes and 
occasional metasediments. The north canyon wall is exposed 
bedrock while the south side has a mantle of soil and 
colluvial material. 

Major joint sets trend approximately N 50° E in this area. 

Leakage along joints and fractures could result in percolation 
through the "gooseneck" or meander. 

The alternate location is at the entrance to the canyon 
near the confluence of Duck and Deer Creeks. The site is also 
in Precambrian rocks. The west bank exhibits several cut 
terraces and is probably covered with a mantle of unconsoli- 
dated alluvium. 
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In the case of either dam site, the water impounded would 
cover an area of low relief called Duck Creek Flats. Several 
types of geologic materials are present in the area to be 
covered by the reservoir. Prefcambrian granites, gneisses 
and schists form the bedrock and are commonly exposed at the 
surface. In some places, however, the bedrock is covered by 
weathered colluvium in addition to Tertiary sediments. 

Although none of the existing maps of the immediate area show 
Tertiary rocks. White River or Arikaree sediments have been 
mapped in the adjacent area of the west fork and probably 
extend into the reservoir area. The low valley drained by 
Duck Creek and West Deer Creek is a likely place to find more 
isolated patches of these Tertiary sediments. A number of 
springs and seeps are present in the reservoir area which 
are controlled by factures in the bedrock. 

Conclusions 

1. ) Both the proposed and alternate sites would provide 
a Precambrian bedrock foundation with minimal excavation of 
alluvium. 

2. ) A number of well-developed joint systems are present 
in the Precambrian rocks. It is recommended that the joints 
be further evaluated with respect to possible leakage. 

3. ) The reservoir would cover areas of Precambrian, Ter- 
tiary (?), and Quaternary rocks. Some springs are present 
along fractures in the Precambrian. The Tertiary White River 
sediments characteristically contain clays which exhibit 
swelling properties. 

4. ) No active faults were recognized in the area. Although 
the area has been tectonically active in the geologic past, 

it is structurally stable at present. 

5. ) No other obvious geologic hazards, aside from the 
joints, were observed in this brief study. 
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Appendix C : Summary of a study of regional geomorphologic 

relationships in the Wyoming Basin province. 


from a special report by R. M. Breckenridge, Brainard 
Hears , Jr., and others, Geomorphic Sections of the Wyo- 
ming Basin Province (in progress). 

The Wyoming Basin geomorphic province (Fenneman, 1931) 
has never been satisfactorily differentiated into its consti- 
tuant divisions using geomorphology as the principle criterion. 
This study attempts to accomplish this subdivision using 
interpretations of ERTS?-1 imagery, aerial photography and 
field data. 

The province was divided into ten basins and associated 

uplifted regions (Fig. C-l). This delineation was based on 

distinctions of geomorphic entities apparent from the ERTS-1 

£ 

color— composite mosaic of Wyoming and upon available geologic 
and topographic data. Topography, geologic structure, and 
stratigraphy were considered along with the geomorphic criteria 
in defining subdivision boundaries. A description of each of 
the basin areas was then compiled setting forth its major 
characteristics and discussing the decision criteria used 
in selecting each basin boundary. Included in each descrip- 
tion is a summary of the geography, size, boundaries, drainage, 
topography, and general geology of the area. The following 
description of the Laramie Basin (one of the ten basin sub- 
divisions) presents a typical set of physiographic and boundary 
conditions upon which the basin subdivisions were made: 

The Laramie Basin of southeastern Wyoming encom- 
passes a north-south trending area 90 miles long 
and 30 miles wide. It lies between the Laramie Range 
on the south and east, the Medicine Bow Mountains 
on the west and the Shirley and Hanna Basins on 
the north. The northeast-flowing Laramie River is 
the major drainage for the area, and the basin con- 
tains extensive alluvial deposits along the river 
and its tributaries. Elevations on the basin floor 
are usually from 7000 to 7500 feet. 

*The 1 :1, 000,000-scale mosaic used in this study was completed 
by the General Electric Photographic Engineering Laboratories, 
Baltsville, Maryland. 
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Structurally the Laramie Basin is a broad north - 
plunging syncline containing about 8000 feet of late 
Paleozoic and Mesozoic rocks with some Tertiary rocks 
in the northwest part of the area. Alluvial deposits 
and terrace gravels of Pleistocene and later age are 
also quite abundant. Several smaller anticlines are 
found on the west and north sides of the basin, but 
are not prominent topographically as are those of the 
basins to the north. 

The eastern boundary of the Laramie Basin is the 
top of the Casper Formation and the base of the prominent 
dip slope developed on the resistant limestones and sand- 
stones of this formation o The southeast boundary, west 
of Tie Siding, Wyoming, corresponds to the Casper-Pre- 
cambrian contact and the mountain front of Boulder 
Ridge. In this area, the base of the Casper dip slope 
is at a relatively low angle and cannot be shown as a 
single line. The southern boundary is also arbitrary, 
as shown by the dashed line, but roughly corresponds to 
the unconformity between the Precambrian granite and the 
Fountain Formation. 

The western boundary is shown as the break in slope 
between the relatively flat basin and the steep mountain 
front of Precambrian from the eastern flank of Bull 
Mountain to the south, and Jelm Mountain and Sheep Mour 
tain to the north. Centennial Valley is included in the 
Laramie Basin, the boundary again being the steep moun- 
tain form of Precambrian in contact with the Fountain 
Formation on the east side of Centennial Valley . A 
fault contact defines the western and northern boundary 
of Centennial Valley to Comer Mountain. 

The north- trending fault-controlled Precambrian 
front of the Medicine Bow Mountains continues to bound 
the west side of the Laramie Basin to Arlington. This 
boundary is dashed northward from Arlington, cutting 
the center of the Cooper Cove anticline which divides 
the Laramie Basin from the Hanna Basin. The boundary 
continues to the north where it intersects the ridges 
of the Mesaverde Group south' of Medicine Bow. The north- 
ern boundary of the Laramie Basin is the top of the 
Cloverly Formation hogback as exposed on Como Bluff 
anticline. The boundary is dashed eastward from here 
to the base of the Casper Formation dip slope on the 
Laramie Range „ This northern boundary is somewhat arbit- 
rary, however, it appears best to use the most prominent 
hogback and the largest feature in the area. In addi- 
tion, the many minor anticlines and synclines in the 
northern Laramie Basin appear to be more subdued than 
those to the north. Como Ridge marks the boundary of 
this characteristic. 
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Appendix D : The use of satellite imagery and aerial photo- 

graphy for making lithologic distinctions. 


summarized from two special reports submitted under con- 
tracts NAS 5-21799 and NAS 9-13298 s 

Houston, R. S., 1973, Distinguishing major litho- 
logic types in rocks of Precambrian age in central Wyo- 
ming using multilevel sensing — with a chapter on possible 
economic significance of iron formation discovered by 
use of aircraft images in the Granite Mountains of Wyo- 
ming, NTIS , Springfield, Virginia, Report E73-10301, 34p. 

Houston, R. S., 1974, Multilevel sensing as an aid 
in mineral exploration — iron formation example; in 
Houston and others, 1974, Some illustrations of the ad- 
vantages of improved resolution in geologic _ studies ; 
National Technical Information Service, Springfield, 
Virginia, Report #74-10628, 73p. 

NASA coverage of Wyoming includes complete cloud free 
ERTS imagery (900 meters above the surface. Fig. D-l) on a 
seasonal basis (18-day cycle, about 1/3 to 1/2 cloud free); 
greater than 50 percent cloud-free Skylab (270 miles above the 
surface. Fig. D-2), and greater than 70 percent aircraft cov- 
erage (15,000-60,000 feet above the surface. Fig. D-3). Stan- 
dard products include film positives and prints of multiband 
images generated from the MSS ERTS multispectral scanner, 
film positives and prints of six bands from Skylab S190A 
cameras, high resolution black and white and color positives 
and pfints from the Skylab S190B camera and a variety of 
products from aircraft sensors including six or more bands 
from multispectral camera systems, high-resolution color and 
color infrared prints and film positives from 10-in. aerial 
cameras and various products from multispectral scanners and 
spectrophotometers . 


Multilevel Sensing 

Th^multilevel sensing approach takes advantage of this 
material to save time and money in geologic studies. Space- 
craft images can be used to target areas of geologic and/or 
economic interest and images taken at lower (aircraft) ele- 
vations can be used to study key areas in greater detail. 
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Figure D-3o Index map showing NASA aircraft coverage of Wyoming as of October, 1973 











... ..,,.1... I I I 1 1 i f 

88 


This approach saves time and money since the user will find 
the space image can be obtained and analyzed at less cost 
per unit area. 

Examples of the use of this approach are l) recognition 
of regional lineaments with spacecraft images and detailed 
study of these features with aircraft products and on the 
ground; 2) recognition of major color anomalies associated 
with mineralization by use of color and color infrared space 
images and verification or rejection of these by use of 
aircraft and ground studies; and 3) mapping of major litho- 
logic units that may be host of mineral deposits by use of 
spacecraft images and study of these units in greater detail 
by use of aircraft and ground studies 0 

Rocks of Precambrian age in Wyoming are primarily Archean 
(2o5 b.y D ) in age and consist of two major terranes — green- 
stone and granite gneiss. The granitic rocks are of economic 
interest because they may be directly or indirectly one of the 
sources of uranium now in the sandstone-type uranium deposits 
of the basins, but. most of the mineral deposits in the Archean 
rocks are in the greenstone belts that make up a relatively 
small part of the total bulk of Archean rocks. The green- 
stone belts contain all known iron formations of economic 
significance and virtually all other metallic mineral deposits 
in the Wyoming Precambrian except deposits in post-Archean 
metasedimentary rocks located in southeastern Wyoming. For 

| 

this reason, the greenstone belts are an economic target. 

The results of a preliminary study using ERTS-MSS band-5 
j for Precambrian terrane distinctions have been discussed by 

Houston and Short (1972, p. 12) greenstone belts could be 
distinguished from granite and granite gneiss terrane in the 
j southeastern Wind River Mountains and to some degree in the 

] Granite Mountains of Wyoming. 

j _ ' _ 

j In the southeastern Wind River Mountains mapping by 

j ■ ■ 

j Bayley (1965a-d) and Bell (1955) gives adequate ground truth 

| to enable us to evaluate the ERTS, MSS-5 interpretation. 

According to Bayley (1965a) the metasedimentary/metavolcanic 

j ■■ 

| 

I ■ 


rock terrane (greenstone) consists chiefly of light -brown 
to black feldspathic and micaceous graywacke, conglomerate 
and mica schist, and dark-colored mafic lava flows and bodies 
of gabbro, black graphitic schist, dark brown meta-andesite 
flows, and lesser amounts of dark quartz -magnetite iron for- 
mation and white, gray, or pale green quartzite and pale 
green mica (fuchsite) schist. The area underlain by these 
rocks is dark colored and distinctly layered. The granitic 
terrane (Bayley, 1965c) may be gray quartz diorite and grando- 
diorite or white to pink coarse-grained granite (Bell, 1955, 
p. 13). The granitic terrane for the most part is light 
colored and texturally more massive than the metasedimentary/ 
metavolcanic rock terrane 0 This tonal contrast shows quite 
well on ERTS MSS-5 imagery with the metasediment ary /metavol- 
canic terrane indicated as areas of low reflectance and the 
granitic terrane as areas of higher reflectance. 

Although similar distinctions can be made in the Granite 
Mountains the interpretation from the ERTS image is not as 
accurate as that of the southeastern Wind River Mountains . 

The accuracy of the interpretation may be judged by comparing 
Figure D-4 (the ERTS interpretation) with Figure D-5 which is 
a geologic compilation from the best available geologic and 
photogeologic mapping of this area. The major area of meta- 
sedimentary/metavolcanic rock in the northeast Granite 
Mountains can be distinguished but minor granite bodies with- 
in this area cannot. In the southwest the general area where 
metasedimentary /metavolcanic rocks are present can be recog- 
nized from the ERTS MSS-5 imagery but an area of metasedimen- 
tary rock approximately two miles wide and six miles long was 
misinterpreted as granite. 

The difficulty in making distinctions between rock types 
in the Granite Mountains is probably the result of several 
factors. The Granite Mountains are interpreted (Love, 1970) 
as a once great mountain area like the Wind River Mountains 
that collapsed to form a graben. The mountains were partially 
exhumed so that today we see isolated masses of rocks of 
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Figure D-4.Map of Precambrian terrane of the southern Wind 
River Mountains (bottom center) and Granite 
Mountains (lower right) showing greenstone belts 
pCms mappable from ERTS imagery. 
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Precambrian age surrounded by flat-lying sedimentary and vol- 
canic sedimentary rocks of late Tertiary age. Some masses 
stand well above the Tertiary floor and others are barely 
exhumed. In general, the granitic rock terrane has higher 
relief, less soil, and, as stated above, is light colored as 
compared to the metasedimentary/metavolcanic terrane. The 
higher reflectance in the red band of granitic terrane is 
probably a combination of all of these factors rather than a 
simple expression of rock color alone. In addition, some 
metasedimentary rock units in the Granite Mountains are more 
felsic in composition and lighter in color than those of the 
Wind River Mountains, and are difficult to distinguish from 
granite even with high-resolution color and color infrared 
photography. The various factors undoubtedly hinder litho- 
logic distinctions in this area. This conclusion is supported 
by limited field checks and available ground truth. 

Figure D-6 is a geologic map prepared from interpretation 
of U-2 red band and color infrared 70 mm positive: transparen- 
cies. Resolution of the U-2 photographs was adequate to dis- 
tinguish many of the large diabase dikes that cut the granitic 
rocks of this area so three Precambrian lithologies were 
mapped ; granite-granite gneiss terrane, greenstone belt ter- 
rane, and diabase dikes. The distinctions made using U-2 
images are good as indicated by a comparison with Figure D-5 . 

Direct comparison of the U-2 photography with the Mission 
184 (C-130) photography shows that only half as many dikes can 
be resolved on the U-2 photography as on the higher resolution 
C-130 photography . This is apparently a line-resolution limita 
tion rather than a reflectance difference because the wider 
dikes are clearly distinguished on the U-2 images. 

Figure D-7a and D-7b compare U-2 and ERTS images of a 
granite-metasedimentary rock contact in the southeastern 
Wind River Mountains. In this application the U-2 has no 
advantage over the ERTS images because the contact can be 
readily mapped with either image. 



Figure D-7» Images of the southeastern Wind River Mountains of Wyoming showing contact between 

Precambrian metasedimentary rock - metavolcanic rock (dark mass on left) and granite/ 
granite-gneiss; red band U-2 (A) and MSS-5 ( B ) . The contact strikes approximately 
north and can be clearly seen in the left center of the images near the fork in the 
Sweetwater River. A contact between northeast dipping Paleozoic rocks and the Pre- 
cambrian is in the northeast portion of each image. 
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Figure D-8 is a geologic map made from interpretation of the 
C-130 aircraft color and color infrared photography o If 
this map is compared with Figure D-5, the reader may note 
that the photogeologic mapping compares favorably with the 
field-checked geologic maps, but several exceptions are appar- 
ent. In T. 30 N., Rs. 92 and 93 W. a series of outcrops of 
gneissic granite, biotite gneiss, and granite gneiss are ex- 
posed in secs. 16, 17, 18, 19, 20, 21, 30, and 31, T. 30 N., 

R. 92 W., and secs 0 1, 2, 11, 12, 13, 14, and 24, T. 30 N., 

R. 93 W. According to Sheirer (1969, p. 7-22), the 
biotite gneiss is a layered rock with alternating biotite- 
rich and quartz-feldspar-rich layers and it forms brownish- 
black to dark-gray outcrops whereas the gneissic granite and 
granite gneiss (both shown as granite in Figure D-5) have less 
well developed layering are generally lighter in color than 
the biotite gneiss. These rocks are mineralogically similar, 
containing chiefly feldspar, quartz, and biotite with the 
primary difference being a greater abundance of biotite and 
better developed layering in the biotite gneiss. They are 
not easily distinguished in the field because all units 
have gradational contrasts and are inhomogeneous in mineral 
percentages. Thus, it is not surprising that these units 
are not readily distinguished by photogeologic techniques. 

The biotite gneiss which is probably sedimentary in origin 
was interpreted as either granite or foliated granite by 
photogeologic techniques and contacts between the various 
units are not correct. 

The area in the general vicinity of Barlow Gap is, perhaps, 
the best for photogeologic mapping in the Granite Mountains. 

The major rock types show strong color contrasts and lith- 
ologies show greater tonal and textural contrast on the film 
positives than in the field. As noted above, iron formation 
was discovered in the area during a field check of the photo- 
geologic map in August of 1972. This is the first report of 
iron formation in the Barlow Gap area, but a small body of 
iron formation was discovered in an area several miles 
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Figure D-5 Reconnaissance geologic map of rocks of Precambrian age in the 

Granite Mountains and adjacent areas of Wyoming compiled by use of 
U-2 and C—130 air photographs and from geologic mapping in T. 32 N. f 

Rs. 87 and 88 W by B O. Carey (1959) and in T. 30 N , Rs 92 and 93W 
by R L. Sherer (1969) Base map modified from Love (1970). 
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northeast of Barlow Gap in 1971 by Alfred Pekarek (personal 
communication, 1971), graduate student at the University of 
Wyoming . 

Three bodies of iron formation (Figure D-9) were noted 
at Barlow Gap in the one day reconnaissance survey of the 
area. The largest occurrence is in a small anticline (as in- 
dicated by primary structures in sedimentary rock noted in 
the field) located in SE 1/2, sec. 16, T. 31 N., Ro 88 W. 

(Fig. D-9). Oxide and silicate facies (James, 1954, p. 249) 
iron formation were recognized. The oxide facies iron for- 
mation consists of layers of magnetite alternating with 
layers of green silicate and brown chert. The silicate 
facies is a laminated greenish black rock probably composed 
of minnesotaite and/or iron-rich amphiboles. Thin layers of 
iron-rich carbonate are also present, but the most abundant 
iron -bearing rocks are the oxide— facies type with alternating 
magnetite and green silicate layers. Chemical analyses of 
these various rock types are m Table D— 1 . The magnetite 
layered iron formation ranges from 21 to 29.5>percent less 
than is typical for this facies as suggested by James (1954, 
p o 249) . 

The iron formation at Barlow Gap and that discovered by 
Alfred Pekarek northeast of Barlow Gap represents a series! of 
small isolated outcrops; none large enough to be of economic 
value at present and many factors; in addition to tonnage 
(Ohle, 1972), must be considered before any firm conclusions 
about the economic value of the iron formation can be made. 
However, the geologic study of this area is only reconnaissance 
and the total areas of metasedimentary rocks where iron- 
formation outcrops are known covers at least 36 square miles . 

In addition to these outcrops Love (1970, C140) reports re- 
sults of drilling which brought up numerous fragments of iron 
formation from the Tertiary Wagon Bed Formation at a depth of 
140-160 feet. These fragments of iron formation in the Ter- 
tiary beds were probably derived from iron formation in the 
Precambrian. The well in which they are found is in NW 1/2 
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Table D— 1 • 


Iron, manganese, titanium and chromium content of samples 
of iron formation from Barlow Gap, Wyoming# (Analyst -- 
J. W. Murphy, University of Wyoming, Rock Analysis Lab). 


Rock Type 


Sample No. %Fe(±0.2) %Mn(±0.003) %Ti(±0.1) %Cr(±0.02) 






Silicate 

RS72-2A 

11.2 

0 o 17 8 

Silicate 

RS72-2B 

31.1 

0.198 

Carbonate 

RS72-3A 

6.1 

0.162 

Magnitite-Silicate 

layered 

RS72-3B 

21.0 

0.213 

Magnet i te -Si 1 i cat e 
layered 

RS72-3C 

29.5 

0.268 

Magnetite Silicate 
layered 

RS72-3D 

21.0 

0.164 
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sec. 14, T. 32 No, R. 86 W. , 16 miles from the Barlow Gap 
area. These observations show that iron formation is present 
in Precambrian metasedimentary rocks in, at least, three 
townships and Fig; D-5 shows that similar metasedimentary 
rocks are present throughout much of the northern Granite 
Mountains. It is the writer* s opinion as was Love's (1970, 
p. C140 ) that this area should be studied to determine the 
extent of these iron-bearing units. 

Perhaps the most useful and least expensive tool for 
initial exploration whould be an airborne magnetometer survey. 
Such a survey might locate iron formation both in surface out- 
crops and in Precambrian metasedimentary rocks that are buried 
beneath the relatively thin cover of Tertiary sedimentary rocks. 



Appendix E : Geologic mapping from ERTS-1 imagery. 


from a special report under contract NAS 5-21799: 

R. S. Houston, 1974, Geologic mapping . using space 
images, in Houston and others, 1974, Some < illustrations 
of the advantages of improved resolution m geologic . 
studies? National Technical Information Service, Spring- 
field, Virginia, Report E74-10628, 73p. 

Photogeologic mapping using standard black and white 
products is basically an attempt to map contacts of photo 
units as expressed by tonal contrasts, texture, and geomorphic 
features, and to map faults and other structures by use of 
stereoscopic interpretation. Tonal contrasts expressed on 
black-and-white prints are useful in mapping lithology, but 
they can seldom be used to identify rock types . Laboratory 
determined reflectance spectra for common rocks show differ- 
ences in the visible and near-infrared regions due primarily 
to electronic processes in the cation or impurity ions (iron 
the most common source) and to vibrational processes in hydroxyl 
groups or molecular water, and in the far infrared (9-12 urn) 
due primarily to silica content (Hunt and Salisbury, 1971). 

The reflectance spectra also show variation in overall reflec- 
tances for rocks depending, in part, on electronic processes 
which can produce an extremely sharp cutoff in the visible or 
near-infrared spectrum, or continuous absorption throughout 
the entire range (Hunt, Salisbury, and Lenoff, 1971, p. 1-3). 
Thus rock spectra show two main characteristics : narrow ab- 

sorption bands due mainly to the presence of iron, hydroxyl 
groups, water and silica, and broad reflectance differences 
or brightness differences. A granite, for example, has a 
relatively high spectral reflectance as compared with basalt, 
basalt may show absorption bands due to iron content not 
shown by granite, and silica absorption bands in granite are 
recorded at different wavelengths than those in basalt (Fig. 
E-l). These spectral differences appear to offer a possibility 
for remote identification of rocks and a number of studies 
have been undertaken, primarily using aircraft systems, to 
determine if rocks can indeed be identified by use of spectral 
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Figure E-l. Comparisons of reflectance and emissivi 
tance curve for the mafic rock shows a 
1.0 pm. The thermal emissivi ty curves 
shorter wavelength due to the increased 


ty spectra for acidic and mafic rocks. The reflec- 
ref lectance low (absorption) at about 0.7 unr and 
show a decided shift in the emissivity low toward 
silica content of the granitic rock. 
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differences. One of the most comprehensive studies was that 
of Lyons (Lyons and Patterson, 1969, p. 527-552) at Stanford 
who has had some success in identifying rocks by taking ad- 
vantage of the shift in silica absorption band with silica 
content in the far infrared. This study made use of an air- 
craft- mounted spectrometer and radiometer . Characteristic 
spectra of common rocks were compared with those of the tar- 
get or unknown as determined by the spectrometer, and under 
optimum conditions rock types and contacts could be mapped 
along the straight-line path of the spectrometer. The method 
is subject to a number of limitations that will be discussed, 
but for geologic mapping a major limitation is that informa- 
tion is confined to the single narrow path of the spectro- 
meter and is difficult to extrapolate into other areas for 
mapping. An approach that uses the same silica shift prin- 
ciple in mapping rocks but makes use of a multispectral 
scanner instead of a spectrometer has been proposed and tested 
by Vincent (Vincent and Thompson, 1971, p. 247-252). Instead 
of using the complete rock spectra as was done by Lyons, selec- 
ted bands of the spectra (i.e. 8.2-10.9 um/9. 4-12.1 urn) are 
rationed and ratio images are produced that are useful in 
distinguishing silica-rich and silica-poor rocks. This 
method is less specific than that of Lyons, but multispectral 
scanners produce image data so that maps can be generated 
from the interpretations and identifications ■ 

The ratio approach has been extended to make use of iron- 
absorption bands by a number of investigators (Vincent, 1972, 
p. 1239-1247; Rowan, 1972, p. 60-1-60-18), and because the 
MSS scanner of ERTS includes near infrared bands the technique 
has been used to map iron-rich rocks and alteration zones by 
use of ERTS data (Vincent, 1973; Rowan and others, 1973). 

These various techniques plus others using related 
approaches (Watson, 1972; Pohn, Of field, and Watson, 1972; 
Watson, Rowan, and Of field, 1971) offer promise for eventual 
quantitization of geologic mapping, but, as is true for many 
geologic problems , the number of variables that must be 
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considered to apply these techniques to actual field problems 
is very large. Factors that affect rock spectra include 
grain size, grain orientation, degree of weathering, vegeta- 
tion cover, sun angle and azimuth, viewing geometry, slope 
angle, temperature, water saturation, atmospheric conditions 
and depth of the atmospheric column, shadowing, and such 
common factors as strike and dip of outcropping beds. No 
technique has been developed to take all of these factors 
into account although some are considered in most experiments • 
Statistical methods that essentially map significant reflec- 
tance differences in a given area or use training sets to 
map a certain tone or density that may characterize a rock in 
a given area (Anuta and others, 1971; Smedes and others, 1972; 
Watson and Rowan, 1971) are useful but, again, cannot be ap- 
plied to general mapping because of variables cited above 0 
The 1974 state-of-the-art suggests that enhancement techniques 
may be useful for mineral exploration under controlled con- 
ditions (i.e 0 where the method is tailored to a specific area) 
but general geologic mapping is still best done by the experi- 
enced photogeologist using standard photogeologic techniques. 

The best approach we have found for geologic mapping using 
space images is to try to recognize units with distinctive 
tone or texture that can be used as marker beds in a given 
region, and to use these beds along with other standard recog- 
nition parameters for mapping decisions . Figure E-2 shows 
laboratory reflectance spectra determined at Goddard Space 
Flight Center (Short and Lowman, 1973, p. 9) for a typical 
suite of Wyoming sedimentary rocks. The spectra are not 
sufficiently distinctive to be diagnostic of rock type, and, 
as noted by Short and Lowman (1973, p. 18), several individual 
rocks displayed greater variation between spectra from fresh 
and weathered surfaces an inch apart than between quite differ- 
ent rock types. Yet it is apparent that one sedimentary rock 
group, the Cretaceous shales, has exceedingly low reflectance 
as compared with other units. These beds of the upper Ther- 
mopolis Shale and lower Mowry Shale contrast strongly with 
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Figure E-2 . Spectral response curves for a group of rocks from Wyoming 
determined on a Carey 90 reflectance spectrometer (from 
Short and Lowman, 1973, Fig* 30, p. 90* 
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other units as the laboratory spectra predict and can be 
recognized through most of their outcrop area 0 Furthermore, 
with the advantage of color composite images, color can also 
be used for mapping. Throughout most of Wyoming, Permo- 
Triassic red beds make excellent marker -units. 

A geologic map of the Arminto area of central Wyoming 
was prepared from interpretation of an ERTS infrared color 
composite using a Bausch and Lomb zoom transfer microscope. 

When compared with the state geologic map, if was found that 
units mapped from tonal variations are surprisingly close 
to those chosen as formations by field geologists — particu- 
larly since most of these major units are rock -vegetation 
units . 

Figure E-3 compares the southeastern guarter of the Arminto 
AMS map with same area on the geologic map of Wyoming (Love 
and others, 1955) . Much of this area has not been mapped in 
detail and no mapping is available in published or open- 
file form. The geologic map shows nine lithologic subidivisions 
whereas the ERTS photogeologic map shows twenty^-five lithologic 
subdivisions. This area has not been fully field-checked so 
we cannot be certain that all subdivisions shown on the ERTS 
photogeologic map correspond to changes in lithology 0 The 
large area shown as Cody Shale (Kc) on the geologic map has 
been subdivided into several units ; the areas shown as Kn ( cs ) 
and Kn (s) probably correspond to calcareous shales and shales 
of the Niobrara Formation. The area shown as Kc is typical 
Cody shale. The area shown as Kstc is the Shannon Sandstone 
member of the Cody Shale. The area shown as Kc (b) is Cody 
Shale with irregular vegetative cover. Finally, the lens 
shown as Kss may be the Sussex Sandstone member of the Cody 
Shale. The reason for the irregular vegetative cover on the 
Upper Cody Shale bed has not been determined but it Kc(b) does 
coincide with the location of a major oil field. The odd-shaped 
area in T. 38 N., Rs. 81 and 82 W. is an area where the natural 
vegetation has been disrupted by extensive spraying to kill 
sagebrush so its significance is probable not lithologic. 
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Figure E-3 0 Details of the southeast 1/A of the ERTS 
Arminto map (A) compared with the^ 
same area on the state geologic map of 
Wyoming .(B, after Love and others) o 
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The Mesaverde Formation can also be subdivided into the 
Parkman Sandstone Member (kmvp) , a shale unit (Kmvs), and an 
upper unit that is probably the Teapot Sandstone member (Kmvtp). 
Distinctions were also made ■within the Fox Hills Sandstone 
(Kfh) and Lance Formation (Kl), but their significance is 
not known. A unit designated ( Tf us ) may be a lower, massive 
sandstone of the Fort Union Formation and a unit designated 
(Tuss) may be a sandstone of the Wasatch Formation. Finally, 
stabilized sand dunes (Qsc), active sand dunes (Qsa), and 
alluvium (Qal) are readily mapped. 

The southeastern quarter of the Arminto AMS sheet is 
particularly well suited for regional geologic mapping with 
ERTS color infrared because the various formations display 
good vegetation control, have low dips, and exhibit a fairly 
broad outcrop area. Other areas are not so suitable for map- 
ping because dips are steeper, structure more complicated, 
and vegetation tends to mask rock types rather than enhance 
them. 


Appendix F : Geologic applications of ERTS-1 in volcanic 

terrane . 

summarized from a special report under contract NAS 5-21799 
Breckenridge , R.M., 1973, Application of ERTS imagery 
to geologic mapping in the volcanic terrane of northwest 
Wyoming; National Technical Information Service, Spring- 
field, Virginia, Report E74-10244, 2lp. 

Two of the largest volcanic provinces in the United States 
are located in Northwestern Wyoming; the Absaroka volcanic 
field and the Yellowstone Plateau. The Absaroka Range forms 
the southern part of the Absaroka field which extends north 
into the Gallatin Range of Montana (figure F-l) . Generally 
the Absaroka Range consists of a large dissected plateau of 
layered Eocene volcanic flows, breccias and conglomerates, 
most of which are andesitic in composition. The Quaternary 
Yellowstone volcanics adjoin and cover the Absaroka field on 
the west. Unlike the Eocene eruptions, the Yellowstone vol- 
canism was not violently explosive and accumulated a lava flows 
and ash flow tuffs. Together the two provinces provide an 
excellent test area of varied volcanic features. 

Purpose and Method of Investigation 

The purpose of this study was to evaluate ERTS-1 imagery 
in interpretation of various geologic features common to this 
region. Volcanic rocks are notoriously difficult to map. They 
generally lack sufficient stratigraphic markers and are charac- 
terized by rapid facies changes over short distances. Also 
much of the section is mineralogically homogeneous. 

Although parts of the region are well-traveled and 
familiar to anyone who has visited Yellowstone National Park, 
most is wilderness area and extremely rugged. For a general 
discussion of the regional geology, the reader is referred to 
Keefer (1972); or, for a more detailed discussion, to Smedes 
and Prostka (1972) and Christiansen and Blank (1972). An ERTS 
interpretation of the glacial history of this same area was 
presented in an earlier report, Breckenridge (1973). A first 
look analysis of the area from ERTS was reported in Houston 
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Figure F-l 
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map of volcanics in northwestern Wyoming. 
Renfro and Feray ). 
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and Short (1973). Subsequent receipt of color conposites and 
application of other interpretive techniques has resulted in 
improved geologic interpretation. 

Laboratory interpretations of the ERTS imagery were made 
using standard photogeologic techniques augmented by color- 
additive viewing and isodensity contouring. Several small 
areas within the main test site were checked j,u. the field 
in order to confirm the image interpretations. Spectral 
reflectances of representative rock outcrops were obtained 
with a multi-channel filtered photometer. Each channel was 
calibrated by reading a reference grey card at the outcrop. 

Data were collected in spectral bands corresponding to ERTS 
bands 4, 5, 6, and 7. Limited time and accessibility permitted 
field checks of only a small portion of the study area, but it 
seems reasonable to assume that the interpretation has approxi- 
mately the same reliability for the entire area as was demon- 
strated in the areas that were field checked. In some areas 
aircraft imagery was available as an additional check. 

At first an effort was made to restrict the map interpre- 
tations to features which any trained photogeologist could 
recognize, but the author’s familiarity with the area pre- 
cluded an unbiased judgement. However, it is apparent that 
much extremely useful information could be derived from the 
ERTS imagery by an interpreter who is more than casually 
familiar with the area under investigation. For the purpose 
of this study, all existing knowledge of the region was used 
to benefit the conclusions. 

Structure 

An ERTS band-7 mosaic was constructed for the area to 
facilitate a uniform regional geologic interpretation of 
the area (figure F-2). Several major structural features 
and numerous smaller structures were mapped in the ERTS image 
interpretation. Most of the obvious structural features are 
in the sedimentary units along the margin of the Bighorn 
Basin. The largest structural feature noticeable on the 
imagery is the Yellowstone caldera (figure F-3) . The caldera 
is so large an ill-defined that it cannot be recognized from 
the ground and, in fact, was not known until the late I950*s 
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(Keefer, p. 34-42). Recognition of the feature on the imagery 
is possible primarily because of the series of topographically 
expressed ring fracture zones outlining the perimeter of the 
structure and hot-spring deposits which are concentrated in 
the old caldera. Also, the flooding of the caldera floor by 
late rhyolites resulted in a conspicuously flat basin. 

A number of fault zones mapped in the Yellowstone Park: 
have been described by Love (1961) as Quaternary fault zones, 
many of which are exposed as recent fault scarps. Other fault 
zones shown only on the newer map of Yellowstone (U.S.G.S., 

1972) were also mapped from ERTS. Direction of movement 
could not be determined from the ERTS imagery. In several 
places, glacial streaming and flow features were confused 
with faulting . 

Nearly all of the linear elements previously discerned 
from radar imagery (see Christianson, Pierce, Prostka, and 
Ruppel, 1966) were distinguished of the ERTS imagery. 

Faulting in the Absaroka Field is not as obvious as in 
the younger plateau volcanics. The extreme dissection of the 
mountains in a dendritic pattern produces a complex topography 
and shaded slopes which combine to mask most linear elements. 
However, some major structural features appear to exert 
control on topography and can be observed on the ERTS imagery. 

The trace of the Wood River fault (Wilson, 1964) was ex- 
tended beyond the original strike on the basis of the ERTS- 
image interpretation. The valley of the South Fork of the 
Shoshone River is distinctly linear and probably structurally 
controlled. Alluvial fans on the north wall coalesce in a 
straight line which probably represents the fault trace. Land- 
slides characterized by hummocky topography dominate the south 
side of the valley obscuring possible structural elements. 

Both the He^rt Mountain and South Fork thrusts (Pierce, 1957) 
were easily mapped from the color composite imagery. The 
thrust sheets are distinguished because of the high reflectance 
of the Bighorn Dolomite, Jefferson and Madison Limestones in 
the upper plate. 
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A number of other "linears", or linear elements sugges- 
tive of faulting were recognized and are shown on the map, 
but were not checked this field season because of difficult 
access and time limitation. 

Lithologic Distinctions 

As previously mentioned, the Absaroka volcanic field is 
a highly dissected plateau of vent- and alluvial volcanic 
facies. Although compositions vary from andesite to basalt 
most rocks are andesitic 0 A wide range of textures occur in 
the volcanic pile. The stratigraphic nomemclature of the 
region has a long and complicated history which is treated 

by Smedes and Prostka (1972). 

The stratigraphy of the Yellowstone Plateau volcanics 
has been recently summarized by Christianson and Blank (1972). 
The more subdued topography and heavy vegetation of the 
plateau area combine to limit exposures of the Yellowstone 
volcanic sequence , but some recent flows can be mapped by 
recognition of flow-lines, fronts, topographic position, and 
pressure-ridges . 

Differentiation of volcanic, sedimentary, igneous and 
metamorphic lithologies was best accomplished by using a 
standard color-composite ERTS image. The Absaroka volcanics 
have a distinctily dark-brown color and appear much as they 
do in outcrop. The more mafic units have a characteristic 
dark blue tone on the color composite ERTS image. Similar 
drab colors identify the Yellowstone volcanic rocks when they 
are not obscured by vegetation. The plateau volcanics are 
most easily differentiated by their flow features, such as 
pressure ridges and steep fronts. Older crytalline rocks 
in the Beartooth and Teton Ranges have high reflectance in 
Band 7 and exhibit numerous joints. The crytalline blocks 
are often bounded by large faults. Sedimentary rocks in the 
Bighorn Basin are characterized by bedding, structure, and 
light tone. Figure F-4 shows the distribution of major rock 
types in the area 0 In some cases, the eastern contact of the 
volcanic rocks was difficult to define because large volcanic- 
covered pediments lap out over the sedimentary units of the 



Figure F-4 Major rock types mapped from ERTS-1 Imagery 
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Lithologic distinctions within the volcanic sequence 
were best accomplished with band 7 images used in conjunc- 
tion with a corresponding color composite. Initial field 
examinations were somewhat disappointing. Units readily 
recognized in the field could not be distinguished on the 
imagery. Careful examination of the problem revealed that 
most of the units were differentiated in the field on the 
basis of mineralogic texture and outcrop form. For example, 
some breccias were coarse and massive while others were fine 
and thinly bedded, and some were chaotic with interbedded 
lenses while others were uniform. Yet all of the mappable 
units are similar in composition and spectral response. In 
addition, many of the better exposures occur in near-vertical 
walls and are not imaged by vertical-looking sensors. The 
small-scale textural differences could not be resolved at 
ERTS scale. However, the cumulative effect of frequent 
ledges, "hoodoos," lenses, flows, etc. results in large 
changes that are detectable on the ERTS imagery. After the 
initial field check the meaningful distinctions reflecting 
different formation surface textures could be more confidently 
identified on the imagery and the interpretations improved 
accordingly . 

At the South Fork test site the Wapiti Formation was 
successfully subdivided into three units similar to those 
mapped by Pierce and Nelson (1969 ) j l) an. upper series 
of flows and flow breccias, 2) a middle mafic breccia, and 
3) a lower, light-colored volcanic sandstone, siltstone, 
and conglomerate. Again, the major difference between the 
flows and breccias is textural or topographic, but the 
fine-grained volcanic sandstone and siltstone has an 
anomalously high reflectance and appears light on the 
image. Spectral response of each rock type was measured 
in the field for the ERTS and EREP bands (Table F-l) . 

All of the values represent an average of at least 5 
readings over a large- outcrop. Since the field of view 
was small (100 sq. ft. or less), these values cannot 
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Latite dike 
(Wood R. ) 

17.0 

18.0 

17.5 

20.0 

17.5 

16.5 

19.0 

Latite dike 
(N. Fork) 

13.0 

20.0 

20.0 

21.0 

20.0 

21.0 

19.0 

dike 

(S. Fork) 

21.0 

25.0 

24.0 

22.0 

21.0 

24.0 

23.0 

Andesite 
( Funnel Mtn . ) 

12.0 

9.0 

11.0 

15.0 

11.0 

8.0 

9.5 


BRECCIAS AND FLOWS 


Early Basic 
Breccia 

17.5 

18.5 

Wiggins Fm. 

12.0 

11.0 

Crosby Breccia 

17.0 

18.0 

Crosby Breccia 
(green facies) 

19.0 

21.0 

Wapiti Fm. 

(Twb) "breccia" 

18.0 

20.0 

Wapiti Fm. 
(Tws) "flows" 

21.5 

27.0 

Wapiti Fm, 

( Tw j ) J im Mtn . 

20.0 

24.0 

Member - 
Trachyandesite 




21.0 

19.0 

20.0 

20.5 

20.0 

11.0 

10.0 

13.0 

13.0 

12.0 

19.0 

21.0 

18.0 

19.5 

20.0 

22.0 

16.0 

24.0 

21.0 

28.0 

20.0 

20.0 

16.0 

19.0 

21.0 

26.0 

22.0 

21.0 

25.0 

24.0 

22.0 

21.0 

19.0 

21.0 

21.0 


Table F— 1 • Filtered Photometer Readings 
From Major Lithologies of the 
Yellowstone-Absaroka Test Site 


: Values represent uncorrected photometer readings 

rather than absolute reflectances., Consequently 
they reflect only qualitative differences between 
the measured units. 
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reflect the overall effect of ledge-shadowing or vegetation. 
Consequently, the field-measurements minimize the weathering 
characteristics of particular rock types and probably 
maximize the compositional similarities. They also demon- 
strate that these rocks have a relatively rlat spectral 
response . 

Multiband photography was taken at the outcrop to 
further evaluate possible enhancement of volcanic litho- 
logic units. Figure F-5 shows a section of pyroclastic 
rocks in the Wood River test site. To the right of each 
image is a densitometric curve for that image. As can 
be seen from the curves, all of the rock units have a 
relatively uniform reflectance through this part of the 
spectrum. Consequently, multiband enhancement of the 
pyroclastic rocks is limited. 

Intrusive Rocks 

The volcanic flows, breccias and tuffs of the region 
are intruded by numerous intrusives, ranging from dikes a 
few feet thick to large plutons. In general, the larger 
intrusive bodies are much lighter in tone on all ERTS bands 
than the surrounding pyroclastics. Parsons (1963) recognized 
seven major vent areas in the Absaroka Range. Figure F-6 
shows the intrusives mapped from ERTS. Since some of the 
intrusions are covered with vegetation or occur on shadowed 
north sloped, they were not recognized on the ERTS inter- 
pretation. A number of intrusions not mapped on black-and- 
white ERTS images (Houston and Short, 1972) were readily 
distinguished on the color-composites. 

Most of the major vents are centers of extensive dike 
swarms. Unfortunately, most of the dikes could not be re- 
solved at ERTS scale. Several of the larger dikes in the 
North Fork Valley were mappable. 

Mineralized Zones 

Several mining districts in the Absarokas are associated 
with intrusive activity. The Kirwin-Wood River area was 
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Figure F-5 Multiband photography of Wood River test site 
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Figure F-6 Intrusives in the Yellowstone - Absaroka Provinces 


121 






122 


chosen as a test site for ERTS imagery interpretation. The 
Kirwin area is one of the major vent centers in the volcanic 
field. Numerous intrusive bodies are exposed in the immediate 
area. Several of these are mineralized. Recognition of the 
mineralized intrusions on ERTS imagery was possible because 
of reddish alteration zones which appear yellow on the color- 
composites. After targets of possible mineralization were 
picked from the ERTS imagery, high altitude aircraft imagery 
was employed in checking the interpretation at a larger scale. 
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Appendix Gs Geologic map of Carbon County , Wyoming interpreted 
” ’ from ERTS-1 imagery . 

from an unpublished reports 

Agard, S. S., 1974, Carbon County geologic map _ 
interpreted from ERTS-1 imagery; University of Wyoming, 
Department of Geology, unpublished special project 
report . 

Purpose 

The purpose of this study was to map the bedrock and sur- 
ficial geology of Carbon County, Wyoming, using ERTS-1 imagery, 
and to determine the accuracy and detail of the information 
obtained (Plate G-l). 

Location 

Carbon County encompasses nearly 8,400 square miles of 
central and southeastern Wyoming (Fig. G-l). It is a geo- 
logically diverse area with features ranging from large Pre- 
cambrian- cored uplifts to small complexely folded structures 
to large areas of flat— lying Tertiary formations. Regionally, 
it includes the northernmost part of the Southern Rocky 
Mountain Province and a major portion of the Wyoming Basin 
Province. The former is characterized by the broad anti- 
clinal uplifts of the Sierra Madre/Park Range and the Medi- 
cine Bow Mountains where large expanses of Precambrian are 
exposed. The Wyoming Basin in this area consists predomin- 
antly of Tertiary-filled basins interrupted and complica- 
ted by smaller uplifts such as the Rawlins, Ferris, Seminoe, 
Shirley and Freezeout Mountains. Older Tertiary formations 
are often involved in these structural disturbances while 
younger Cenozoic units are undisturbed <• Clastic Quaternary 
deposits, including large areas of sand dunes, playa lakes, 
evaporites, glacial drift and alluvium occupy much of the 
surface area. 

Methods 

Using a Bausch and Lomb zoom transfer scope (ZT-4) 
positive 9" by 9” ERTS transparencies were superimposed on 
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a base map (scale 1; 250, 000) compiled from the AMS 2° series. 

The images used were: 

1408 - 17235 4 Sep 73 

1353 - 17190 11 Jul 73 

1030 - 17242 22 Aug 72 MSS Color Composite 

1030 - 17244 22 Aug 72 MSS Color Composite 

diazo color composite of 1353 - 17190 
Color composites were used when available because the con- 
trasting colors of the units provide the best means for dif- 
ferentiating them. In areas of low color-contrast such as 
the outcrop areas of Upper Cretaceous and Tertiary units in 
the western portion of Carbon County , two identical color 
composite images were superimposed and used as one in order 
to increase the depth of color and contrast. Black and 
white images were used for areas in which the color composite 
coverage was obscured by cloud cover and for the eastern 
quarter of the county where color composites were not avail- 
able. Band 5 was used most extensively in these areas since 
it usually yields the greatest tonal contrast between litho- 
logic units. Band 7 has slightly less contrast but is par- 
ticularly useful in the heavily vegetated mountain areas which 
are too dark on the band 5 images. A diazo color composite 
of the eastern portion of Carbon County yielded good informa- 
tion when used in conjunction with the black and white images. 
However, the resolution of the diazo composite was too poor 
for accurate work in geologically complex areas □ 

Descriptions of Units 

The units described below are identified by both their 
map unit number and the name of the corresponding geologic 
formation. The photo map units were correlated with the 
formations using prior knowledge of the geology of the area 
and by relating the photomap units to the State Geologic 
Map of Wyoming (Love and others, 1955) and the stratigraphic 
nomenclature chart (Figure G-2). These correlations are 
tentative because they have not yet been field checked for 
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SOUTH CENTRAL WYOMING STRATIGRAPHIC NOMENCLATURE CHART 



Figure G-2. Stratigraphic nomenclature chart for the Carbon 
County area (after Brinker and Blackstone, 1951) 
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accuracy. The nomenclature used, is that of the 1955 Geologic 
map of Wyoming although some of the units have since been 
renamed or reclassified 0 

Xn general , the Paleozoic and Mesozoic units are exposed 
only on the flanks of the Precambrian-cored uplifts or in 
complexly faulted areas nearby o They are consistently thin 
units and are sometimes difficult to map in detail. They 
are usually distinct in character and color from adjacent 
units, but often the characteristics of a unit in one area 
do not carry through to other areas. The differences are 
primarily a function of locally different structural atti- 
tudes, elevations, vegetal cover, angle to sun, soil moisture 
and related factors. In such cases, the succession of units 
may be inferred from the location of a few consistently iden- 
tifiable units. The Upper Cretaceous and Tertiary units 
are more uniformly characteristic and could often be directly 
correlated from one area to another • 

The descriptions presented are based purely on the appear- 
ance of the units on the images used. Thus, they are tied to 
the particular times of year that the images were recorded 
and to the specific color processing used in creating the 
color composites. Different descriptions and somewhat dif- 
ferent results might be obtained using other images. 

Unit 1 - Precambrian is exposed in the core of the Sierra 
Madre, Medicine Bow, Ferris, Seminoe, Shirley, Pedro, and 
Elk Mountains. The Precambrian rocks generally have a crystal- 
line appearance, exhibit few sedimentary features and are 
heavily vegetated. They are characteristically deep brownish- 
red in color. Smaller features, like glacial drift, alluvium 
or sedimentary remnants within the area, cannot be distinguished 
against the dark background. 

Unit 2 - (Cambrian) is found in contact with Unit 1 on the 
flanks of the Rawlins, Ferris and Seminoe Mountains. The 
photomap unit may actually include some thin units of the 
younger Madison Formation. Xn the Ferris and Seminoe Moun- 
tains the Cambrian rocks have a deep red color (darker than 
Unit 1) in sharp contrast to lighter colored younger units. 
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In the Rawlins area. Unit 2 is a lighter, mottled, reddish- 
brown color. This color seems to reflect sparse tree cover. 

Unit 3 - (Carboniferous Madison Fm. ) occurs in the Raw- 
lins and Shirley Mountain areas. In the Rawlins region, it 
is the lightest colored unit in the area as contrasted to 
Unit 2 and the overlying Tertiary unit. In the complex 
structure of the Shirley Mountains, Unit 3 forms a rela- 
tively barren light colored marker unit contrasting with 
heavily vegetated Units 1 and Unit 4. 

Unit 4 - (Carboniferous Tensleep Fm. and Amsden Fm., ) is 
a marker unit used in correlating the sequence of lower units 
in different areas. It appears as a prominent, thin, white 
band in the Rawlins, Ferris and Seminoe Mountains and to the 
west of Elk Mountain . It is also striking, though somewhat 
darker in tone, due to vegetation in the Shirley and Freezout 
Mountains . 

Unit 5 - (Carboniferous undifferentiated) includes units 
3 and 4 described above where they could not be individually 
distinguished. These areas, such as the region north of the 
Medicine Bow Mountains, were heavily vegetated giving a red- 
dish hue to Unit 5 . 

Unit 6 - (Permian Phosphoria Fm. or Forelle and Satanka 
Fm.) is located in the Rawlins, Ferris, Seminoe and Freezout 
Mountains, and near Slate Ridge. In all areas, it has a 
yellowish-green color which is usually greener than surround- 
ing formations. Variation in actual hue and sharpness of 
the contact occurs from region to region. 

Unit 7 - (Triassic undifferentiated) represents the red- 
bed sequence of the Rawlins, Ferris and Freezeout Mountains. 

It has a characteristic yellow color with a slight green 
tint in some areas. Unit 7 includes the Chugwater, Dinwoody, 
and some parts of the Phosphoria Formation. 

Unit 8 - (Jurassic Undifferentiated) is located in the 
Ferris, Rawlins, and Slate Ridge uplifts. The light turquoise 
blue color is fairly even and forms a distinct contact with 


other "units. The Nugget Sandstone, and the Sundance and 
Gypsum Springs Formations make up this unit. 

Unit 9 — (Jurassic Morrison Fm. and Cretaceous Cleverly 
Fm.) is found in the same areas as Unit 8. Its color varies 
from a mottled, light greenish-yellow in the Rawlins area to 
light, blue-green elsewhere. It is generally lighter than 
surrounding units . 

Unit 10 - (Jurassic-Cretaceous Morrison, Cloverly & 

Sundance Fms.) includes Units 8 and 9 in areas where they are 
indistinguishable or too thin to be mappable. Unit 10 was 
mapped in the Freezeout Mountains and Como Ridge. It appears 
light turquoise in color. 

Unit 11 - (Cretaceous Muddy & Thermopolis Fms.) was 
mapped as two subunits in the Rawlins -Grenville Dome area. 

The older unit (11a) is dark blue while unit (lib) is a 
lighter blue-grey. The two subunits were indistinguishable 
in the Ferris and Freezeout Mountains and by Slate Ridge and 
Como Ridge. Here Unit 11 is light green in color. 

Unit 12 - (Cretaceous Frontier Fm.) also appears as two 
units in the Rawlins -Grenville Dome area„ Both are fairly 
consistent light blue with unit 12a being slightly darker. 
Elsewhere it is mapped as one unit, somewhat more blue in 
color than Unit 11. 

Unit 13 - (Cretaceous Niobrara Fm. ) was mapped in the 
northeast part of the county along the outer flanks of the 
Shirley Mountains, Slate Ridge and Como Ridge. It is highly 
mottled with brown, tan, and blue colors, and does not always 
form a distinct contact with younger units. 

Unit 14 - (Cretaceous Cody Fm. or Steele Fm. ) forms wide 
topographic low areas around the outer flanks of the uplifts. 

It is found throughout the county as a very light colored, 
sometimes mottled unit. It contrasts sharply with the gener- 
ally well defined bluish colored older units and the prominent 

brown ridges of Unit 15. 

Unit 15 - (Cretaceous Mesa Verde Fm.) occurs as a prominent 
ridge system around the major uplifts. It is easily traced 
from one region to another and so forms a good marker unit 
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for correlation. In some areas as many as six distinct subunits 
can be seen (as in the Sheep Mountain-Syverson Flats area) . 

These become less distinct and grade into one another as the 
local structure changes; steep dips decrease the number of 
mappable subunits . All subunits of the Mesaverde Formation 
are different shades of brown and grey. Unit (15c) is gener- 
ally darkest in color and most continuous in character from 
region to region. Other subunits are lighter* more mottled, 
and regionally more variable. 

Unit 16 - (Cretaceous Lewis Fm.) is found throughout the 
region in contact with Unit 15. It forms areas of low relief 
adjacent to the ridges formed by Unit 15 and is usually very 
light brown in color. 

Unit 17 - (Cretaceous Lance Fm. ) was mapped in the west- 
ern part of the county as a gently dipping hogback extending 
from near Lamont to Baggso It is mottled light brown and tan, 
and is sometimes indistinguishable from surrounding units. 

Unit 18 - (Cretaceous Medicine Bow Fm.) is located in the 
areas surrounding the Hanna Basin. To the east of the Hay- 
stack Mountains, it is a light tan color with a slight olive 
tint. In the Saddleback Hills area, it is a much darker 
brown -grey and more mottled in appearance. 

Unit 19 - (Paleozoic-Mesozoic undifferentiated) groups 
the above units (2 - 19) in areas too small to show detail or 
where individual units were not distinguishable. 

Unit 20 - (Tertiary Ferris Fm.) flanks the structural 
depression of the Hanna Basin. Its coloring is similar to 
that of Unit 19 but is generally darker than the underlying 
units. 

Unit 21 - (Tertiary Ferris Fm.) extends from Separation 
Flats south towards Baggs . In places as many as three dis- 
tinct subunits are seen. The lower unit is consistently 
light grey-brown through the area. The middle and upper units 
vary from a dark brown in the north to a bluish-green or 
bluish-brown in the south. The middle unit is usually darker 
than the upper unit but the distinction is gradational in 
places. 
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Unit 22 - (Tertiary Hanna Fm. ) is the most prominent 
unit of the Hanna Basin. Its appearance is highly mottled 
with blue, green and brown. Variations in color indicate 
subunits within Unit 22 and were used to trace possible 
structures within the area. Some of the apparent structures 
may be related to the local topography rather than true 
structure. Stereoscopic study of these areas might help 
solve this ambiguity . 

Unit 23a - (Tertiary Wind River Fm.) was mapped in the 
northeast part of the county around Bates Hole and the Shir- 
ley Mountains o It is a mottled blue grey and tan as con- 
trasted with the more even coloration of surrounding units. 

Unit 23b - (Tertiary Wasatch Fm. ) is located along the 
flank of the Washakie Basin forming a large flat lying 
region. It is mostly yellowish brown with several discon- 
tinuous units of darker blue and brown which may reflect 
subtle changes in lithology. 

Unit 24 - (Tertiary Green River Fm. - Tipton member) 
forms the outer unit rimming the Washakie Basin. It is a 
relatively dark blue grey with even coloring and texture. 

To the northwest it becomes somewhat lighter in color. 

Unit 25 - (Tertiary Wasatch Fm. - Cathedral Bluffs 
member) is also located on the flanks of Washakie Basin. 

The prominent bright yellow color is characteristic in the 
area. To the northwest the ridge appears to show two map- 
pable units, one bluish-grey and the other a mottled bright 
yellow. 

Unit 26a - (Tertiary Green River Fm. - Morrow Member) 
forms a highly dissected ridge (badlands?) in the Washakie 
Basin. The comparatively dark brownish-green color appears 
mottled due to local topography. 

Unit 26b - (Tertiary Green River Fm. - Laney member) is 
a slightly yellowish buff color in contrast to the much darker 
underlying units o The coloration is more even, suggesting 
less dissection of this part of the ridge. 
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Unit 27 - (Tertiary undifferentiated) includes most of 
the unfolded flat-lying units in the county. In the north 
unit 27 is dark grey with blue or red hues locally present 
along the flanks of the Ferris, Seminoe, Pedro, Shirley and 
F'reezeout Mountains. To the south. Unit 27 is a lighter, 
eyenly colored brown where it covers the structure in the 
lower units near Rawlins, Cedar Ridge and Pass Creek Ridge. 

Unit 2 8 - (Tertiary Browns Park Fm 0 ) covers a large area 
west and north of the Sierra Madre Mountains. The dark 
brown to brownish red unit exhibits low dips and covers the 
underlying structure except where streams have dissected it. 

It forms a prominent scarp along the north and western edges, 
particularly in the Miller Hill area. To the east it is over- 
lain by Unit 31 but the contact is not distinct. 

Unit 29 - (Tertiary North Park Fm.) is seen in the Sara- 
toga Valley. Its contact with the older units of the Sierra 
Madre and Medicine Bow Mountains is sometimes sharp and some- 
times gradational. It is generally a light brown or blue- 
grey color contrasting with the very dark red of Unit 1, 

Where Unit 29 is heavily vegetated the exact contact with 
Unit 1 is not apparent. 

Unit 30 - (Tertiary intrusives) is located west of the 
Sierra Madre Mountains near Battle Mountain. The intrusive 
bodies stands out above Unit 28 as darker colored circular 
buttes . 

Unit 31 - (Quaternary glacial drift) is located north of 
the Medicine Bow Mountains near Kennaday Peak. It is a red 
and brown mottled color and grades into surrounding units. 

Unit 32 - (Quaternary playa deposits) is found near the 
Rawlins uplift and Haystack Mountains 0 The playas are bluish- 
grey or white in color and are often associated with small 
lakes (also playas?) in the flatland areas of Unit 14 and 
other Quaternary units . 

Unit 33 and 33d - (Quaternary dunes) covers a large area 
between the Rawlins uplift and the Ferris-Seminoe Mountains. 
The striped areas on the map denote probable active sand and 


wind erosion. The dunes appear very white in color and show 
a dominant northeasterly trend. Other areas include stabilized 
dunes with some vegetation. Stabilized dune areas appear 
pockmarked like a karst surface and are generally dark tan or 
brown in color. 

Unit 34 - (Quaternary alluvium) is associated with modem 
streams throughout the county . The alluvium appears light tan 
or bluish-grey in some areas and red (due to vigorous vegeta- 
tion in others. Where Unit 34 is vegetated it is so designa- 
ted since it is often impossible to tell from the imagery if 
cropland areas along streams are actually underlain by bedrock 
or alluvium. This distinction was particularly difficult in 
the heavily farmed Saratoga Valley. Field checking is needed 
to determine which areas are actually alluvium. 


Structure 

In general, the broad structures apparent on the ERTS 
imagery were easily mapped. The major uplifts of the Rawlins, 
Ferris, Seminoe, Shirley and Freezout Mountains along with 
associated faults and structural complications were mapped 
in considerable detail. The major unconformities were Ter- 
tiary and Quaternary units overlie older units were readily 
discerned. However, structural detail interpreted from the 
ERTS imagery is not always consistent with resolution. In 
some areas , more detail could be seen than could be mapped at 
1 : 250,000 (i.e 0 the small area west of Elk Mountain). In 
other areas, very little structural detail could be seen (i.e. 
the blank area on the map east of Rock Mountain) . In the 
latter area, the color variations were discontinuous and the 
structure could not be interpreted. In some areas, a single 
unit could be subdivided many times whereas elsewhere it ap- 
peared only as one. 

This variation in mappability is controlled by a number 
of factors, including the attitude of the units, their rela- 
tion to sun angle, the elevation, vegetation cover and climate 
More subunits were generally distinguished in areas of shallow 
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dip; fewer were apparent in steeply dipping and structurally 
complex areas . In heavily vegetated areas , all units were 
characteristically red in color and difficult to distinguish. 

In areas sloping gently away from the sun azimuth enhancement 
served to increase the detail somewhat. 

Objectivity 

There are many built in biases in any geologic interpre- 
tation but they are perhaps more prevalent in a photo inter- 
pretive map. The factors described above tend to increase 
the detail interpreted in one area while limiting the inter- 
pretable detail in another. The imagery used was not uniform 
throughout the county. Some areas had the added advantage of 
color-composite presentation whereas the eastern part was 
available only black and white. Cloud cover in some areas 
resulted in piecemeal interpretations compiled from several 
images . 

The order in which the areas were mapped affected the 
end product. The Rawlins uplift and surrounding area was 
interpreted first. The interpreter was relatively unfamiliar 
with the map area, the data, and the interpretive procedures. 
Consequently, much time was spent mapping this small area. 

The initial attempt bogged down in detail, necessitating an 
overall assessment of the area for "mappable -units.'’ Also, 
the criterion and techniques used for interpreting the first 
areas mapped evolved as the study continued 0 

Perhaps the greatest bias in the interpretation was the 
correlation from the state map and the author's previous know- 
ledge as to what should be seen in each area. The sequence 
of units, their position and role in major structures of the 
county and their expected topographic expression were all known 
before the interpretation was begun. In general, the state 
map was not used for correlation of units until after the 
interpretation for each region was completed, but the struc- 
tural style and characteristics of each area checked cer- 
tainly influence interpretation of subsequent areas. 


Future Work 

The Carbon County photomap should be carefully checked 
against all available published data and detailed aerial 
photography. Areas of major discrepancies should be field- 
checked to determine the cause of the discrepancies. Until 
these essential checks are made, the information presented 
on the photo map and the correlations with geologic formations 
and structure must be considered tentative. 

Conclusions 

The ERTS-1 imagery of Carbon County yields good detail 
for large scale mapping. Photo linear elements that may 
represent faults and joints are shown in greater detail than 
on the Geologic Map of Wyoming. The areas of Tertiary out- 
crop exhibit subtle differences not always evident on the 
ground. Other areas showed less detail than has already 
been mapped. A thorough check of the photomap is the next 
essential step in determining validity of the interpretations 
and the utility of the ERTS imagery in this application. 
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Appendix H s Utility of high-altitude color photography for 
geologic mapping. 


from an unpublished master's theses s t 

Vargas H., Rodrigo, 1973, Photogeo logic map of 
the Jelm Mountain Quadrangle , Albany County , Wyoming ; _ 
University of Wyoming, Dept, of Geology, M. Sci. thesis, 
68p. 

Color infrared film positive (scale ls39600 approx.), 
color film positive (scale ls39600 approx.), and blach and 
white panchromatic prints (scale 1:27700 approx.) were used 
for photogeologic mapping of the Jelm Mountain quadrangle. 

The color and color infrared films correspond to two flight- 
lines oriented approximately N 35° W. The black and white 
prints are part of three N-S trending flight lines. The 
black and white photography was taken by the United States 
Geological Survey in July, 1947 and the color photography 
by the National Aeronautics and Space Administration in 
September, 1972 (mission 213). 

Prints for the color and color infrared photography were 
available but were not used because they contained less infor- 
mation than the corresponding transparency images, largely 
because of overexposure and the resulting loss of color con- 
trast. The same lithologic units were differentiated on the 
three photographic products used= But, it is necessary to 
keep in mind that: l) the black and white panchromatic photo- 
graphy was larger scale; 2) the black and white panchromatic 
film was flown at a different time than the color infrared 
photography and yields a somewhat different perspective on 
most slopes; and 3) the determination and delineation of 
these lithologic units were more quickly and easily made on 
the color infrared film positive than on either of the other 
photographic products. 

When comparing color and color infrared films with the 
black and white panchromatic prints, it is observed that the 
color mode allows an easier determination and delineation of 
the stratigraphic units present in the studied area. This 


i 


1 


i 


J 


137 


color mode is still more valuable on the color infrared film 

because of its higher contrast. 

Because color infrared film presents the greatest tonal 
and color contrast, it is considered the most valuable (in 
this particular case) of the three studied photographic 

types o 

The high tonal and color contrast allows faults and 
some lithologic units (Forelle Limestone, Morrison Formation) 
to be more readily and precisely delineated on the infrared 
photography „ Color contrast between the different lithologic 
units allows a ready observation of their offset along the 
fault line. Quaternary deposits are easily delineated by 
their geomorphic form, and their color is largely dependent 
on the color of the parent material from which they were de- 
rived. Pediment gravels are characteristically portrayed by 
a grayish green to dusky, yellow-green color with a spotty 
pattern due to the vegetation (sagebrush). This spotty pat- 
tern could be due to a "patterned ground" which is one of the 
most conspicuous and characteristic features of periglacial 
regions (Embleton, and King, 1968, p. 485). 

The color infrared photography allows ready identifica- 
tion of key stratigraphic units when used in geologic mapping. 
The key units in this area ares The red beds that are charac- 
teristically portrayed in yellow or yellow-green on the color 
infrared photography? the black shales which appear darker 
than any other sedimentary rock, and the light-colored sand- 
stones which are depicted in white on the imagery. 

Combining these key units with general stratigraphic 
knowledge of the area, a geologist can construct a reasonably 
accurate geologic map in a considerably shorter time than is 
required for completing the same work using standard photo- 
graphy and field mapping. Nevertheless the utility of the 
color infrared photography must not be over-emphasized. Its 
use does not allow the field work to be avoided but rather, 
increases the overall efficienty of the operation by reducing 
field work to a minimum. 
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Appendix I s Mapping of large glacial features using ERTS-1 
imagery . 

a summary of a paper presented March 5-9, 1973 at the 
Symposium on Significant Results Obtained from ERTS-1 s 
Breckenridge , R. M. , 1973, Glaciation of north- 
western Wyoming interpreted from ERTS-1 imagery; Sym- 
posium on significant results obtained from ERTS-1, 

March 5-9, 1973, Washington, Do C. , Report NASA-SP- 
327, NASA Scientific and Technical Information Office, 
Washington, D. C. , p. 363-370. 

During the Pleistocene epoch the largest mass of glacial 
ice outside of the continuous ice sheet in North America formed 
in northwestern Wyoming leaving a profound effect on the pre- 
sent day topography and scenery. In addition the sand and 
gravel deposits formed by glacial erosion are important sources 
of construction materials as well as being excellent aquifers. 
Glaciation also provided numerous natural reservoirs in the 
form of tarns and moraine- dammed lakes. Rocky Mountain 
glacial chronology is based on type deposits found in the 
Wind River Mountains of Wyoming. Although glaciation has lo- 
cally been studied in detail, the region as a whole has not 
been studied from a glacial standpoint. The large area in- 
volved, rugged topography, lack of control and general inac- 
cessibility have discouraged regional investigations. ERTS-1 
appeared to provide an excellent base for broad-scale geo- 
morphic studies of this type. 

Figure 1-1 is a small-scale photomosaic of ERTS-1 MSS-5 
imagery in the area studied. The darker area with moderate 
relief and the prominent lake in the northwest comer of the 
image is the Yellowstone Plateau. It is bordered on the east 
by the rugged Absaroka Range and on the south by the Teton 
and Gros Ventre Mountains . The large northwest-striking 
structure in the south-center of the mosaic is the Wind River 
Range. The north-trending ridges and valleys on the southwest 
margin of the image are the Wyoming Overthrust Belt. All of 
the images were taken in early fall (1972) and are cloud-free. 
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The imagery was analyzed using basic photo interpretive tech- 
niques together with color-additive and pseudo-stereoscopic 
viewing o 

Overlay maps were compiled outlining the extent of ice at 
glacial maximum. A simplified and reduced glacial map is pre- 
sented in Figure 1-2 . Glaciated areas as shown are based on 
the recognition of one or more glacial features. Glacial 
landforms can be grouped into erosional and depositional 
features both of which could be recognized on the imagery . Mor- 
aines are considered the most diagnostic depositional glacial 
landform. The Bull Lake and Pinedale type deposits (Wisconsin) 
were easily outlined on the imagery although they could not be 
separated where in close contact. These massive looping ter- 
minal moraines impound large glacial lakes on the flanks of 
the Wind River Range. 

Terminal and lateral moraines in mountain valleys could 
neither be reliably distinguished from stratified drift nor 
from other valleyside surficial deposits except in a few cases 
where characteristic geomorphic features were recognized. 

Most of the pre-Wisconsin glacial features have been subse- 
quently destroyed by erosion or covered by later deposits. 
Remaining exposures are too small to recognize on the imagery. 
Large gravel-capped mesas or pediments at the mountain flanks 
and in the basins probably represent erosional cycles of these 
early glaciations. These pediments or benches appear dark on 
all bands but are best recognized using a composite of MSS 5 
and 7. 

Where massive terminal moraines were absent or not recog- 
nizable, glaciated areas were outlined on the basis of glacial 
topography. Alpine glacial erosion is distinctive and stereo- 
scopic viewing of adjacent ERTS passes facilitated the recog- 
nition of U-shaped troughs, cirques, horns, aretes, cols and 
truncated spurs. Ice limits can be located approximately where 
valleys change from a U to V shape in cross valley profile. 

Most of the crests of Wyoming mountain ranges culminate in flat 
uplands often referred to as high-level erosion surfaces which 
bevel across the' rock structure. Glacial dissection of these 
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surfaces results in the classic biscuit-board topography. Mor- 
phology and topography recognized using stereo pairs allowed 
distinction of the unglaciated areas above the ice. Color com- 
posites enhanced the contrast between tundra vegetation on 
the erosion surfaces and the bare cirque walls. Glacial tarns 
distinguished on MSS-7 were useful as indicators of glaciated 
bedrock. 

Major flow directions (Fig. 1-2, arrows) were mapped using 
directional criteria such as slope direction and orientation 
of glacial troughs, U-shaped valleys, and cirque headwalls. 
Scoured bedrock served as an excellent indicator in some areas. 
An example of flow direction recognizable on ERTS imagery is 
in northern Yellowstone Park where ice moving from the Buffalo 
Plateau can be traced down the valley of the Yellowstone River 
to Gardiner, Montana. A combination of drumlins and lateral 
moraines together with scoured bedrock in the form of large 
grooves and oriented tarns show on the imagery as parallel 
lines. Linear features are expecially well developed on Mt. 
Everts near Mammoth where ice must have been over 2 , 000 feet 
deep in order to cover the top of the mountain. 
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Appendix J s Extimating available water contained in snowpack. 

an unpublished class report by: 

R 0 w 0 Marrs and J 0 Thompson, University of Wyoming, 

Water Resource Research Institute. 

The May 21, 1973 ERTS scene of the Wind River Range 
(1302-17362, Fig. J-l) was first interpreted to produce a map 
of snow-cover at the time of overpass. But the distribution 
of snow, by itself, is seldom an adequate index of available 
snowpack and cannot be used directly to estimate runoff. 

Snow depth and snow water-content data are required in order 
to make meaningful estimations . Accurate estimates of snow 
depth and density plus several years correlation between 
ground measurement of runoff and the snow-cover estimate 
made from ERTS scenes are needed. 

Determination of the areal distribution of snow from 
the ERTS scenes is straightforward except in scenes with 
cloud cover. The clouds can usually be distinguished from 
snow, but they do obscure the boundaries of the snow-covered 
areas . 

Color-additive presentation of snow scenes often allows 
the interpreter to identify areas where the snow is melting 
or has exceptionally high moisture content. However, it is 
not always possible to separate differences due to sun angle 
and slope variations and from differences reflecting moisture 
content or partial vegetation cover. Spectral comparisons 
made by color-additive viewing or by computer analysis can 
help in making these distinctions. Combinations of ERTS MSS 
bands 5 (red) and 7 (infrared) are often most useful. Both 
bands record the area of total snow cover , but the infrared 
reflectance of the snow seems to vary with moisture content 
more than does the red-band reflectance. Consequently, a 
false-color combination of these bands shows variations in 
color corresponding to varying moisture conditions on the 
surface of the snow layer. 

A series of snow— cover estimates and maps were compiled 
for the upper Green River drainage (Table J-l, Figs. J-2, J-3, 



An ERTS scene of the Wind River Mountains, north- 
ern Thrust Belt and Absaroka Range. The upper 
Green River drainage is in the center of the image 
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Table J-lt Snow cover and snow depth estimates made from 
ERTS-1 imagery 


DATE 

IMAGE 

SNOW 
COVER 
( total_L 

MELTING 

SNOW 

EST. 

DEPTH. 

CLOUDS 

15 Jan. '73 

1176-17355 

100% 

0 

20” 

0 

2 Feb. '73 

1194-17362 

100% 

0 

25” 

0 

20 Feb. *73 

1212-17363 

100% 

0 

38" 

0 

10 March *73 

1230-17364 

100% 

? 

20-40” 

60% 

15 April '73 

1266-17364 

100% 

5% 

10-40” 

0 

* 3 May '73 

1284-17363 

99% 

5% 

0-40" 

15% 

*21 May '73 

1302-17362 

65% 

46% 

0-30” 

0 

* 8 June '73 

1320-17361 

33% 

8% 

0-10” 

0 

26 June '73 

1338-17360 

5% 

5% 

0-5" 

60% 

1 Aug . '73 

1374-17353 

0 

0 

0 

5% 


♦Estimates made from maps (Figs. J-2, J-3, and J-4) 
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Figure J-4 June 8 snow-cover map of the upper Green River drainage 


Table J-2o Daily discharge records for the upper Green River drainage (from the Warren 
Bridge gage station)*, ( January-September 1973, discharge in cubic feet 
per second) . 


DAY 

JAN 

FEB 

MAR 

APR 

MY 

JUN 

JUL 

AUG 

SEP 

1 

108 

101 

135 

108 

242 

824 

1 , 840 

604 

271" 

2 

105 

105 

132 

110 

222 

1,080 

1,680 

568 

299 

3 

100 

110 

130 

115 

214 

1,070 

1,440 

562 

306 

4 

97 

113 

128 

120 

284 

944 

1,190 

616 

302 

5 

95 

115 

128 

130 

400 

796 

1,130 

664 

292 

6 

90 

115 

132 

135 

687 

775 

1,140 

712 

284 

7 

100 

114 

130 

130 

1,010 

936 

1,170 

705 

285 

8 

120 

110 

125 

125 

960 

-1,290 

1,190 

659 

308 

9 

130 

108 

120 

120 

568 

1,640 

1,150 

592 

363 

10 

140 

107 

125 

114 

430 

1,930 

1,080 

538 

383 

11 

150 

110 

130 

139 

445 

2,120 

984 

508 

447 

12 

150 

115 

132 

139 

475 

2,150 

920 

490 

600 

13 

154 

118 

125 

147 

586 

1,400 

1,010 

483 

719 

14 

158 

115 

120 

139 

747 

1,640 

1,080 

469 

664 

15 

155 

115 

122 

-147 

976 

1,780 

1,010 

451 

556 

16 

143 

117 

122 

155 

1,360 

1,680 

866 

435 

483 

17 

140 

120 

121 

134 

1,760 

1,390 

747 

430 

425 

18 

138 

122 

118 

144 

2,310 

1,060 

691 

461 

375 

19 

133 

124 

114 

142 

2,460 

838 

691 

461 

336 

20 

130 

125 

110 

114 

2,490 

700 

775 

468 

308 

21 

129 

123 

115 

110 

-2,550 

617 

1,000 

450 

292 

22 

128 

121 

115 

110 

2,000 

635 

1,240 

450 

276 
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Table J -2 ( cont . ) 


DAY 

JAN 

FEB 

MAR 

APR 

MAY 

JON 

JTJL 

AUG 

SEP 

23 

128 

121 

111 

119 

1,580 

778 

1,200 

465 

272 

24 

125 

125 

107 

132 

1,460 

997 

968 

440 

284 

25 

123 

120 

103 

169 

1,660 

1,210 

810 

411 

292 

26 

120 

125 

104 

187 

1,690 

1,340 

691 

410 

280 

27 

113 

130 

105 

194 

1,330 

1,510 

646 

380 

268 

28 

108 

135 

104 

250 

1,030 

1,670 

622 

350 

256 

29 

104 


100 

264 

904 

1,770 

628 

322 

250 

30 

102 


98 

250 

838 

1,840 

640 

298 

246 

31 

100 


105 


789 


634 

278 


TOTAL 

3,821 

3,279 3 

,666 

4,392 

34,457 

38,910 

30,843 

15,130 

10,721 

MEAN 

123 

117 

118 

146 

1,112 

1,297 

995 

488 

357 

MAX 

158 

135 

135 

264 

2,550 

2,150 

1,840 

712 

719 

MIN 

90 

101 

98 

108 

214 

617 

622 

278 

246 

AC-FT 

7,580 

6,500 7 

,270 

8,710 

68,350 

77,180 

61,180 

30,010 

21,270 

GAL YR 

1972 

TOTAL 248 

,615 

MEAN 679 

MAX 4, 

760 MIN 

91 AC-FT 

493,100 


WTR YR 

1973 

TOTAL 164 

,648 

MEAN 451 

MAX 2, 

550 MIN 

90 AC-FT 

326,600 


NOTE . — 

-No gage-height 

record Dec. 14 

to Mar. 

29. 





♦Records taken from U. S. Geological Survey, Water Resources Data for Wyoming, 
Part 1, Surface Water Records, 1973. 
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Table J-3. Surface water runoff forecasts for the Green River at Warren Brrdge- 
April to September, 1973 (1000 acre feet), (from U.S.G.S., 1973). 


MONTH 

FORECAST 

1973 

% OF 
ACREAGE 

FLOW 

1972 

1953-67 

AVERAGE 

FLOW 

1973 

DIFFERENCE 
FROM ACTUAL 

Feb. 

277 

86 

431 

323 

267 

+10 

March 

270 

84 

431 

323 

267 

+ 3 

April 

252 

78 

431 

323 

267 

-15 

May 

238 

74 

431 

323 

267 

-29 


Table J-4. 

Snow dep 

th and water 

content measurements. 

(from U.S.D.A. 

, 1973) o 



1973 



1972 


MONTH 

SNOW 

COURSE 

SURVEY 

DATE 

SNOW 

DEPTH (in) 

WATER 

CONTENT (in) 

WATER CONTENT 
1972 AVG. 

Feb. 

9F21 

1/29 

30 

7.0 


13.9 



10F15 

1/28 

26 

6.2 

Avg . depth 

15.8 



10F16 

1/30 

34 

7.2 

29.6 in. 

19.4 

12.0 






Avg. moisture 

1 A 1 



10F17 

1/29 

29 

6.8 

6.58 in. 




10F19 

1/28 

29 

5.7 


11.6 

8.6 

March 

9F21 

2/26 

35 

8.7 


16.9 



10F15 

2/25 

34 

7.4 

Avg o depth 

19.6 







35.6 in. 

o o cr 



10F16 

2/26 

41 

9.2 


23 • 5 







Avg. moisture 

1 Q ^ 



10F17 

2/26 

36 

8.3 

7.94 in. 




10F19 

2/25 

32 

6.1 


14 o 4 

10.8 
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Table J-4 (cont.) 


1973 1972 


MONTH 

SNOW 

SURVEY 

SNOW 

WATER 

WATER 

CONTENT 

COURSE 

DATE 

DEPTH (in) 

CONTENT 

(in) 

1972 

AVti • 

April 

9F21 

3/30 

43 

10.1 


19.3 



10F15 

3/28 

36 

7,8 

Avg . depth 

19.5 



10F16 

3/27 

46 

11.8 

42.2 in. 

27.1 



10F17 

3/30 

46 

10.5 

Avg. moisture 
9.76 in. 

21o4 



, 10F19 

3/30 

40 

8.6 


16.9 

12.6 

May 

9F21 

4/27 

32 

10.5 


15.6 



10F15 

4/28 

23 

8.8 

Avg . depth 

13.5 



10F16 

4/27 

41 

12.0 

31.6 in. 

26.8 



10F17 

4/29 

29 

9.5 

Avg . moisture 
9.68 in. 

16.5 



10F19 

4/28 

33 

7.6 


14.4 

11.8 


1 
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and J-4) using both the band 5 and band 7 ERTS images. These 
estimates were then compared to surface water runoff records 
(Table J-2) and runoff forecasts (Table J-3) to determine the 
potential for using snow-cover estimates interpreted from ERTS 
with snow co rse measurements (Table J-4) to estimate runnoff 
or available water in snowpack. A direct correlation could 
not be made in either of these applications. The surface 
water runoff measurements include contributions from snowmelt, 
rain, and from springs. Conversely, some of the meltwater 
evaporates and some sinks into the ground, so that only a por- 
tion of the water from melting snow goes directly into the 
surface system and is measured at the gage. Comparisons of 
the ERTS derived data with surface measurements indicate a 
very complex relationship • 

Work is presently underway to better define these rela- 
tionships through consideration of other factors such as temp- 
erature and precipitation and through correlations of ERTS 
image for 1972 and 1974 with appropriate records. Continued 
correlation of the repetitive ERTS data with seasonal and 
daily surface measurements may eventually reveal key relation- 
ships that are not immediately apparent in the results from a 
single season. 
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Appendix K s ERTS imagery applied to mapping of sand dimes. 

A detailed discussion of this study is presented in a 
special report under contract NAS 5-21799 and as part 
of a Geologic Survey of Wyoming report 

Kolm, K. E., 1974, ERTS MSS imagery applied to 
mapping and economic evaluation of sand dunes in Wyo- 
ming, National Technical Information Service, Spring- 
field, Virginia, 31 p. 

Kolm, K. E., 1974, ERTS MSS imagery applied to 
mapping of sand dunes in Wyoming; Geol. Survey of 
Wyoming, Report of Investigation 10, p. 34-38. 

Active and stabilized dune fields of regional extent 
were mapped using Earth Resources Technology Satellite 
(ERTS) imagery (Fig. K-l). Previously mapped dune fields 
(Ahlbrandt, 1973; Houston, 1974; Love, Weitz, and Hose, 

1955; Toehler, 1969) were confirmed by the ERTS image study, 
and some new dune fields were discovered. Additional 
confirmation was provided by field work and high— altitude 
aerial photographs which were available for some areas. 

Results indicate that color composite ERTS images are 
most helpful in locating active dune fields. This is 
attributed to subtle color differentiation between active 
dunes, clouds, sandy alluvium, pediments, alkali flats, 
snowf ields, and light-colored rock formations. The 
color composite images were also employed in mapping sta- 
bilized dune fields which usually exhibited a brown-green 
coloration characteristic of sparse and relatively dormant 
vegetation o The contrast between both stable and active 
dunes and the surroundings was most apparent on the image 
transparencies because the transparencies have higher 
resolution and greater flexibility in the intensity of 
illumination than do the prints. Color prints were most 
convenient for field use. Another general advantage of 
ERTS imagery is the display of regional linear patterns and 
topographic features , Crosscut relations between dune 
fields and surrounding rogk formations were also recognized. 

The disadvantages of ERTS imagery ares l) the imagery 
is of a small scale, and small dune fields were sometimes 


taken fnom 



SAND DUNE FIELD MAP OF WYOMING 

COMPILED USING 

ERTS MULTISPECTRAL SCANNER IMAGERY 


Figure K-l. 


Map of sand-dune fields shows migration directions of stabilized 
(solid arrows) and active (dotted arrows) dunes. Apparent wind 
directions are indicated by large, dashed arrows. 
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missed; and 2) there is sometimes a lacte of distinct 
color and tonal contrast between a stabilized dune field 
and its surroundings. 

The mapping with ERTS was evaluated by comparing the 
ERTS map with a similar map prepared using high-altitude 
aircraft photography. The comparison shows an impressive 
similarity. 

It was hypothesized that the differences between active 
and stabilized dunes might be. a record of a shift in wind 
directions from prehistoric to recent times. More research 
is needed to determine when this shift occurred. 
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Appendix L: Tectonics and linear elements interpreted 

from ERTS. 

A progress report summarizing applications of ERTS 
imagery to regional tectonic analysis, from a repor 
in preparation by J • Earle 0 

Several fracture pattern analyses "have been done on 
the Wyoming Province (Sales, 1968; Stone, 1969) and have 
been correlated with Moody and Hill's (1956) wrench fault 
fracture pattern study. A tectonic study of ERTS photo 
linear elements is being compiled for comparison with these 
published fracture-pattern and tectonic analyses. The 
linear features are compiled on a mosaic of black-and- 
white winter imagery, covering the region extending from 
latitudes 40°N to 45°N and longitudes 104°W to 112 W. 

The 1 inears were first mapped on black-and-white 
transparency images and then transferred to a mylar overlay 
of the mosaic. Large linears that traverse more than 
one image were mapped directly on the overlay. The azimuth 
and length of each linear was measured and rose diagrams 
of directional data are being compiled for interpretation 
and correlation of each province. A separate rose diagram 
, will be prepared for each major area of Precambrian outcrop, 
each major basin, the Absaroka volcanic plateau, and the 
Wyoming Thrust Belt. The linear trends in each area 
will be compared with the linear trends in other areas 
(Precambrian with’ Precambrian, basin with basin, Precambrian 
with basin, etc.). -A second overlay of the mosaic showing 
general or inferred linear trends obtained visually from 
previously compiled mosaics (summer imagery) has also been 
prepared. Linears on this mosaic will also be considered 
in the final interpretation. 

From a preliminary examination of the data, it appears 
that a valid correlation between linear trends mapped on 
ERTS imagery and empirically-derived tectonic trends is 
highly possible. Comparison of ERTS photo linear trends 
with those found on Skylab imagery, aircraft imagery, 
and geologic and topographic maps is planned. 
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Appendix 2d : Structural geology of south-central Wyoming. 

from a special report under contract NAS 5-21799. 

Barton, R., and Marrs, R. W., 1974, Location of 
geologic structures from interpretation of ERTS-1 
imagery. Carbon County, Wyoming; National Technical 
Information Service, Springfield, Virginia, Report 
E74-10333, 8 p. 

ERTS image 1066-17242 (Fig. M-l) covering Carbon 
County and adjacent areas was studied with a primary view 
to mapping structures in the sedimentary rocks of the Hanna 
Basin . 

A color-composite image was available for the area, 
so the interpretation was made both on the basis of tonal 
and textural contrast seen on the black-and-white ERTS-1 
image and color information available from the color composite. 
The overlay map (Fig. M-2) depicts structural and lithologic 
relations interpreted from the image set. After the interpre- 
tation was completed, the Geologic Map of Wyoming (Love and 
others, 1952) was used as a guide in naming the major 
lithologic units identified in the image interpretation. 

Direct comparison was then made of the image interpretation 
and the geologic map (generalized for the Hanna Basin area 
in Fig. M-3). The units mapped on the interpretation 
were readily correlated with major units on the state map. 
Several features (numbered 1 through 6 on the image interpre- 
tation were not indicated on the geologic map, but were 
marked by the interpreter as significant geologic features 
warranting further investigation. The numbered features 
described below are found on the photogeologic map (Fig. 

M-2) 'at the locations indicated by the corresponding numbers. 

1 . Although emphasis in the study was placed on 
lithologic relationships, some of the more prominant linear 
features were also noted. The east-west lineation west of 
Seminoe Reservoir, interpreted as a fault from the image, 
does not appear on the state geologic map° 

2 . To the northeast of the Washakie Basin the ERTS 
image indicates a structure, possibly domal, in the Lewis 



ERTS Image No. 1066-17242-6 covering Carbon 
County and adjacent areas. The present stud} 
was concerned with the eastern portion of 
the image, which consists mostly of folded 
sedimentary rocks of the Hanna Basin, 
basin is bounded to the north by the Sweet- 
water Uplift, to the west by the Rawlins 
Uplift, and to the southeast by Medicine 
Bow Mountains. The Sierra Madre Mountains 
occupy much of the southeastern part of the 
image. The Washakie Basin lies in the south 
west comer of the image area. 


Scale : miles 


Fault 

Medicine Bow Formation 
Lewis Shale 
Mesaverde Formation 
Precambrian 

Figure M-2. Geologic map showing structural features and 
lithologic relation between Medicine Bow, 

Lewis, and Mesaverde formations in the Carbon 
County area. Prepared from a portion of ERTS 
image 1066-17242 and colour composite 1066-17242 






Scale ' miles 



ruuiT 

Tertiary, undivided . 
Medicine Bow Formation 
Lewis Shale 
Mesaverde Formation 
Precambrian 


Geologic map showing relation between Medicine 
Bow, Lewis, and Mesaverde formations in Carbon 
County by J. L# Weitz and J • D» Love, 1952 
(Covers the same area as fig. M-2). 


Figure M-3. 
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shale. There is no evidence of this on the geologic 
map. 

3 0 Just north of center on the ERTS image, immediately 
east of Grenville Dome, occurs a structure which is interpre- 
ted as anticlinal. Again, this is not shown on the geologic 
mapo 

4 & 5. On the nose of the major synclinal structure 
east of Feature 3 and due west of Elk Mountain, in the 
Medicine Bow formation, two minor folds can be seen on 
the ERTS image; they are not shown on the geologic map. 

6. North northwest from Elk Mountain, the ERTS image 
suggests a major anticlinal fold in the Medicine Bow for- 
mation. The geologic map shows only Tertiary cover. 

Two of these features have since been checked against 
other data (nos. 6 and 3). Feature 6 was checked against 
work presented in a U. S. Geological Survey Bulletin (Dobbin 
and others, 1929). Attitudes reported for the rocks in 
the area discount the existence of an anticlinal structure. 
Comparison of the image interpretation with a topographic 
map of the area indicates the apparent structure is a 
combined effect of outcrop and topography. 

Feature 3 was field-checked. The investigator found 
the area covered by soils derived from the Cody shale (the 
dominant lithology) . A few scattered outcrops of the thin 
sands were found. The strike and dip data gathered from 
these outcrops were not excellent, but were sufficient 
to show that beds are conformable to the nearby Grenville 
Dome structure. There was no evidence of a separate dome in 
the anomalous area. Evidence was found suggesting that 
the anomalous area is sometimes covered by standing water— 
that it is a playa lake. The anomalous light-covered 
spot, interpreted as a circular outcrop pattern, corresponds 
to the relatively saline soil of the playa. 

The other four features have not yet been field checked. 

Generally favorable comparison of the Geologic Map 
and the ERTS-1 image interpretation demonstrates that the 


ERTS-1 imagery serves very well us reconnaissance mapping 
base. Topographically controlled outcrop patterns or soil 
anomalies can easily be misinterpreted as structure. 
Stereoscopic coverage would eliminate much of the ambi- 
guity that leads to such misinterpretation. Similar mapping 
efforts in other parts of Wyoming have demonstrated that 
color-composite imagery, and low-altitude aircraft photo- 
graphy greatly augment contruction of an accurate geologic 
reconnaisance map. The experience of the interpreter is also 
a very important factor. The experienced geologist/interpreter 
can use his knowledge of the regional structure and strati- 
graphy and his refined interpretive ability to maximize 
the information gained from the image interpretation and 
construct a more complete and accurate geologic map. 
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Appendix N : Mapping of major fractures with ERTS-1 imagery. 

from a special report on research under contract NAS 
5-21799. 

Parker, R. B., 1972, An evaluation of ERTS-1 
imagery for mapping major Earth fractures and related 
fractures; National Technical Information Service, 
Springfield, Virginia, Report E72-10349, 7 p. 

The primary objective of the study was to evaluate the 
suitability of ERTS imagery in mapping of regional fracture 
patterns in Precambrian crystalline rocks. Such large- 
scale features are commonly difficult to map in the 
field because of their great extent o The synoptic overview 
afforded by satellite imagery appeared promising for studying 
these features. The Wind River Mountains of Wyoming (Fig N-l) 
are an excellent test site inasmuch as an extensive Pre- 
cambrian terrane is well exposed over about 6,000 square 
•kilometers in a relatively unweathered condition. 

Aircraft imagery (Mission 184) was used for a regional 
reconnaissance of a strip approximately 6 km wide, crossing 
the Wind River Mountains from west to east. Broad distinc- 
tions were made among a group of crystalline rocks including 
gneisses, migmatites, and granites of diverse character. 

In addition, many linear features were noted which appeared 
to be fractures or faults of regional extent. 

The reconnaissance strip was field checked during August 
and September of 1972. Access to the area is limited to 
horseback or foot travel via a limited network of trails. 

All portions of the test strip could not be examined because 
of difficulty of access, but sufficient localities were 
visited to generally confirm the conclusion made from the 
aircraft imagery 

The linears interpreted from the aircraft images proved 
to be major zones of dislocation, cataclasis, and metasoma- 
tism is equally varied and ranges from weak recrystalli- 
zation mostly accompanied by epidotization to more completely 
altered rocks composed almost entirely of epidote. In 
two cases, replacement by pink potassium feldspar is complete 



ERTS-1 image 1014-17352-5 of the 
Wind River Mountains. Linear fea- 
tures were interpreted from this and 
other ERTS images of the same area. 
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with local development of rocks composed entirely of 
that mineral. 

The origin of these fractures is not clear at this 
point. Certainly movement has taken place along the fractures 
as evidenced by the cataclastic rocks in the fracture 
zones o They are not merely major joints. However, -the 
existence of major faults cannot be demonstrated because 
there are no rock units obviously displaced by the fractures. 

A more detailed study of the fractures should resolve this 
question. 

Figure N-2 shows the fracture pattern as mapped from 
the ERTS-1 imagery. For comparison, fractures mapped by 
previous workers which were confirmed by the ERTS-1 imagery 
are shown in Figure N-3. A number of fractures and dikes 
mapped by Granger and others (1971) could not be found on 
the satellite imagery. Conversely, their work did not 
indicate some very prominent linears observed on the ERTS-1 
imagery. The large linear indicated by the arrow in Figure 
N-l is the most striking feature in the entire range on 
the ERTS-1 imagery, yet is of uncertain significance at 
this time because it is neither a fault nor a lithologic 
contact on the map of Granger and others (1971). That 
area could not be visited in 1972 because of extreme 
inaccessibility, but a ground study is planned for 1973. 
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Fiaure N-2. Principal mapped faults, shear zones, 
and dikes also observed on ERTS-1 
imagery (from Parker, 1972, Figure 1). 
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from a special report on research under contract 

NAS 5-21799. _ , . 

Blackstone, D. L., 1973, Analysis of photo 
linear elements, Bighom-Pryor Mountains, Montana 
and Wyoming j National Technical Information Service, 
Springfield, Virginia, Report E73-10704, 15 p. 


In a continuing study of photo linear elements* 
in the exposed Precambrian rocks on the core of several 
Wyoming mountain ranges the following ERTS images were 
examined : 

1013 - 17291 - 5 5 August, 1972 (Fig. 2) 

1030 - 17235 - 5 22 August, 1972 (Fig. 3) 

1211 - 17302 - 5 19 February, 1973 (Fig. 4) 

1013 - 17291 - 4, 5, 7 Color Composite 

All well defined linear photo elements were plotted 
on overlay sheets at the same scale as the imagery (approxi- 
mately 1:1,000,000). The elements observed range from 
a few miles to tens of miles in length. The August imagery 
is essentially free of cloud cover whereas the February 
imagery has complete snow cover. Certain elements were 

accentuated by the snow cover . 

Fifty-two major linear elements were located on a 
map with some geographic data (Fig. 0-1 ). Included in the 
plots were monoclinal features involving the Lower Paleo- 
zoic sedimentary rocks on the assumption that they reflect 
basement features, as has been demonstrated by previous 
ground-mapping in both the Pryor and Bighorn Mountains. 

- No data were taken from the folded sediments in the 

Bighorn Basin proper. No ground check has been made 
subsequent to the examination of the ERTS imagery. However, 
the writer has had extensive field experience in these areas 

♦Note : The term photo linear element is used here to 

describe any markedly linear element observable 
on the imaaery irregardless of cause. In this 
particular' situation the linear features have a 
geomorphic expression. 
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The azimuth of fifty-two (52) linear photo elements 
were plotted on a rose diagram in order to visually portray 
the concentration and orientation of the elements o A 
scale of 1 cm equals two (2) percent of the total popula- 
tion was used. Figure 0-2a is the final plot. 

The plot indicates that the majority of the observed 
linear elements have a northeast-southwest orientation 
with a bimodal distribution » One concentration of azimuths 
lies between longitudes N 55° E and N 65° E and the other 
concentration centers near N 25° E. The two sets thus 
intersect at an acute angle of approximately 45 degrees . 

A lesser concentration of linear elements has a north- 
west-southeast orientation with the azimuths concentrated be- 
tween N 10° W and N 30° W. 

Many investigators have studied the problem of basement 
control upon the overlying sedimentary succession during 
the Laramide orogeny . Chamberlin (1945), in a classical 
paper, discussed the problem as it applies to the Bighorn 
Basin region. One facet of this problem is the role that 
fractures, shear zones, etc. in the Precambrian basement 
rocks might have played in Laramide deformation by being 
reactivated. 

If the observed photo linear elements are ancient 
Precambrian features, and if the pattern, in the basement 
rocks continues westward beneath the sedimentary rocks 
of the Bighorn Basin, then the possibility exists of testing 
the relationship between basement control and orientation 
of Laramide folds., 

"he writer utilized two maps prepared by the U. S. 
Geological Surveys 

Oil and Gas Investigations Preliminary Map No. 3. 
"Structure contour map of the Big Horn basin, 

Wyoming and Montana" by David A. Andrews, William 
G. Pierce and Jewell J. Kirby, 1944 

Oil and Gas Investigations Map OM 182. "Structure 
and contour map of the Tensleep Sandstone in the 
Bighorn Basin, Wyoming and Montana" by A. D. Zapp, 
1956. 
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From these maps the trend (azimuth) of 83 known and 
well defined anticlines in the Bighorn Basin was determined, 
and then plotted on a rose diagram in the same fashion 
that data concerning the photo linear elements were plotted. 

The plot of the azimuth of the axes of the folds is 
presented in Figure 0-2b. The dominant trend direction is 
northwest-southeast with a strong maximum (24%) concentrated 
between N 40° W and N 50° W* 

By comparison with Figure 0-2a, it is readily apparent 
that the trend of the anticlinal axes lies at an angle of 
about 80° to the principal concentration of trends of the 
photo linear elements. The only coincidence of trends is 
at a bearing of about + or - five degrees from true north. 
Approximately 7% of the photo linear elements and approxi- 
mately 9.6% of the fold axes have this trend. 

The preliminary conclusion is that the photo linear 
elements reflect major fracture systems in the Precambrian 
age basement rocks, but that they do not directly control 
the orientation of folds in the sedimentary rocks overlying 
the basement in the Bighorn Basin. 

A further conclusion is that the folds in the Bighorn 
Basin have a remarkably consistent trend of the anticlinal 
axes, that is related in some as yet undetermined fashion 
to the Laramide downfolding of the basin as a whole. 
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Appendix P ; Analysis of photo linear elements, Laramie 
Mountains , Wyoming . 

from a special report prepared under contract NAS 
5-21799. 

Blackstone, D. L., 1973, Analysis of photo linear 
elements, Laramie Mountains, Wyoming; National Techni- 
cal Information Service, Springfield, Virginia, Report 
E73-10663, 16 p. 

The essentially north-south trending Laramie Mountains 
developed during the Laramide orogeny, but have many associated 
northeast trending folds of considerable magnitude. The 
relationship of the strongly divergent structural trends 
is not readily apparent and is the subject of this study. 

Using both satellite and aircraft imagery the terrane 
was examined for linear structural features. The linear 
features are, for the most part, expressions of rock frac- 
tures rather than layering within the rock bodies. Re- 
cognition of linears was greatly facilitated by light snow 
cover. 

All such features were plotted on acetate overlay 
sheets, the strike (azimuth) determined, the features counted 
and the results displayed in rose diagrams as frequency 
of occurrence o No field check has been made to date* 
Correspondingly, the results must be considered as first 
order reconnaissance, subject to further check. 

The general distribution of lithologies of the Laramie 
Mountains is shown on Fig. P-1, modified from Condie (1969), 
Newhouse and Hagner (1957) and Darton, Blackwelder and 
Siebenthal (1910). All linear features taken from ERTS 
images are shown on Fig. P-2. Comparison of Figures 
P-1 and P-2 indicates that the underlying rock type does 
not control the distribution nor the trend of the linear 
features . 

The rose diagram plots of the azimuths of 720 linear 
features of the southern Laramie Range imagery is shown 
on Fig. P-3a. The plot defines three principal trend 
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Figure P-3a. Rose diagram showing 
dominant orientations of photo linear 
elements of the southern Laramie 
Range 


Figure P-3b. Rose diagram showing domi- 
nant orientations of photo linear elements 
of the northern Laramie Rasige. 
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directions of the linear features which are believed to 
be fractures. The three dominant directions ares 


Quadrant 


Concent rat ion 

Azimuth 

20 E. 


8.2% 

20° 

40 to 50 

E. 

14.5% 

0 

o 

in 

0 

+> 

0 

o 

70 to 80 

W. 

13.7% 

280° to 290° 


Granite and granite gneiss terranes totally dominate the 
northern Laramie Range, consequently there was no opportunity 
to establish any relationship between rock type and variation 
in orientation of linear features. 

The rose diagram plot of 160 linear features of the 
northern Laramie Range shown on Figure P-3b reveals a 
pronounced northeast orientation. The principal orientation 
is s 

Quadrant Concentration Azimuth 

N. 40 to 50 E. 31 o 8% 40° to 50° 

With a limited number of linear features having an orienta- 
tion of s 

Quadrant Concentration Azimuth 

No 20 W. 6.3% 340° 

In addition to the measurements described above, and 
made from interpretations of ERTS satellite imagery at 
an approximate scale of 1:000,000, a similar set of measure- 
ments was made from interpretations of aerial photography 
(color) at a scale of approximately 1:16,000. To accomplish 
this comparison, three small areas were selected to be 
analyzed through interpretation of the aerial photography . 

Data from the three small domains follows. 

Domain A - T. 17 N., Rs 0 71 and 72 W. 

167 Readings 

Domain B - T. 26 N., R. 72 W. 

169 Readings 

Area of granite and granitic gneiss 

Domain C - T. 30 N. , R. 75 W. 

280 Readings 

Area of granite and granitic gneiss 
No domain entirely within an area of schistose rocks was studied. 
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In Domain A the 1 observed data is summarized below. 


Ouadrant 
N. 45 E. 

Concentration 

Azimuth of Trends 

7.2% 

45° 

N. 80 E. 

10.1% 

80° 

N. 60 W. 

11.4% 

300° 

Data from Domain B 

is summarized below. 

Ouadrant 
N. 5 E. 

Concentration 

Azimuth of Trends 

2.9% 

5° 

N. 50 E. 

26.6% 

50° 

N. 35 E. 

4.1% 

85° 

N. 65 W. 

5.3% 

295° 

Data from Domain |C 
Ouadrant 
N. 15 E. 

are summarized 
Concentration 

below. 

Azimuth of Trends 

7 d% 

15° 

N. 50 E. 

17.5% 

50° 

N. 85 W. 

5.0% 

275° 

N. 60 W. 

11.1% 

300° 


Figure P-4 is a composite rose diagram of 616 readings 
in the restricted domains. 

A comparison of the data from the two sets of imagery 
(satellite and underflight) shows very little variation. 

The number of observations from the restricted domains is 
greater, indicating that smaller scale features can be 
seen from the underflight that are not detectable from the 
satellite imagery. It is also apparent that the minor 
fractures agree will in orientation with the major ones. 

Hamed Bekkar (1973) compared the orientation of the 
imagery features with the orientation of folds and faults 
for the entire Laramie Basin, and demonstrated the close 
agreement of the two sets of observations . Therefore, it is 
concluded that the northeast trending linear elements 
control the orientation of numerous Laramide folds in the 
areas marginal to the Laramie Mountain uplift. Examples of 
such folds are Como Bluff anticline, McGill anticline, and 
Flat Top anticline. 




The type of rock within the Precambrian terrane 
does not appear to influence the orientation of the 
fractures (assuming that the photo linear features are 
the expression of bedrock fracturing) • 
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Appendix 0 : Geology of photo linear elements. Great 

Divide Basin, Wyoming. 

from a special report for contract NAS 5-21799. 

Blackstone , D. L., 1973, Geology of photo linear _ 
elements. Great Divide Basin, Wyoming, National Techni- 
cal Information Service, Springfield, Virginia., 

Report E7 3-1095 8 /WR, 14 p 0 

The Great Divide basin is a broad structural depression 
lying east of the Rock Springs Uplift, south of the Wind River 
Mountains , and Sweetwater Uplift, and east of the Rawlins 
Uplift. The southern boundary is the poorly defined 
Wamsutter Arch. 

The structural basin does not fully coincide with the 
area of interior drainage that includes what has been called 
the Red Desert, Great Divide depression, etc. The Continental 
Divide separates into two strands in T. 27 N., R. 101 W., which 
diverge widely and enclose an area of about 4000 square 
miles before rejoining into a single divide near Divide 
Mountain in T 0 18 N., R. 81 W. The area of internal drainage 
is characterized by extensive playa lake deposits, sand 
dunes, silt and clay dunes, saline lakes, and very few 
permanent streams. In part this area of internal drainage 
may result from the blocking of eastward -flowing drainage 
by sand dune migration in the eastern part of the area. 

The features examined appear on ERTS images : 

1013-17300 (bands 4, 5, 7 and color composite) and 
1030-17242 (bands 4, 5, 7 and color composite). 

Figure Q-l outlines the major photo linear elements 
and can be compared with Fig. Q-2 to see the relationship to 
the major tectonic elements. 

Immediately north of U. S. Highway 1-80 (U. S. 30), 
at Tipton Siding, the imagery shows an area of strongly striped 
ground. The area is about 20 miles in length and about 
10 miles in width. First examination of the imagery suggests 
some type of fault control, but ground elevation examination 
reveals that the striped ground is the site of linear 
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Figure Q-1 




VLeucite Hills volcanic rocks Scale 1:1,000,000 

Major structural elements of the Great Divide Basin, 
Wyoming. 




Figure 0-2 




strips of stabilized sand dunes of moderate relief. Within 
each dune belt strip are blowouts and short transverse 
dunes. 

The dunes are stabilized by greasewood ( Sarcobatug 
sp.) rabbit brush ( Chrvsothamus sp.) and other minor species 
which have a contrasting color, and also are markedly different 
in color value from the sagebrush ( Artemesia sp.) which 
is the common shrubbery cover of the region. A more de- 
tailed treatment of the relationship between Wyoming sand 
dimes and the stabilizing vegetation is presented in separate 
ERTS-1 special report (Kolm, 1974). 

The alignment of the dunes is about N. 75 E. This 
trend is more northeasterly than the prominent and presently 
active Killpecker dune field in the northern part of the 
basin. The distribution of the dimes is limited on the 
east by the rising series of escarpments encircling the 
northern flank of the Washakie Basin, The dune area 
essentially terminates along the eastern end in the Red 
Desert playa • 

A second, more subtle, photo linear feature (location 12, 
Fig. Q-l ) crosses the area in a N. 80 E. direction beginning 
near Point of Rocks siding and continuing northeasterly 
to a point about due north of Creston, Wyoming. The linear 
element appears on the color positive transparency image 
as a faint line which at times appears lighter than the 
surrounding terrain and at other places appears darker. 

Several traverses were made in a north-south direction 
across the linear element in the area between Table Rock 
and Creston Junction, north of U. S. Highway 1-80. Specific 
ground control was sought to establish the nature of the 
linear elements. 

Particular attention was given to geomorphic form, 
stream alignment, observable bedrock geology which might 
reveal faulting, vegetational changes, distribution of 
playa deposits , and aeolian features . 

The result of the traverses was negative with regard 
to critically defining the nature of the Point of Rocks 


photo linear element. No positive faulting was found. Minor 
faulting with a trend of approximately N. 65 E. has been reported 
by Pipiringos (1961) and by Masursky (1962) in the Great Divide 
basin. These faults are reported to have displacements up to 
70 feet. Unfortunately this system of fractures does not coin- 
cide in position with the feature under consideration. The 
low relief and widespread Quaternary deposits made detailed 
examination of the bedrock very difficult. Previous mapping 
is undoubtedly detailed enough to define the bedrock structure. 

The Great Divide basin area is characterized by internal 
drainage and the phenomena associated with such drainage situ- 
ations <> To date, no single paper has discussed geomorphology 
and origin of the internal drainage in this basin. In fact, 
only in 1864 did the topographic base map of Wyoming reflect 
the dual divide. The most complete presentation of the geology 
is that in the Hydrologic Atlas (Welder and McGreevy, 1966). 

The geomorphic history is complex, and warrants study as to 
its origin. 

The playa-lake flats are the site of a limited amount of 
lacustrine deposition, but the depressions in which the play as 
exist must be due in part to an earlier cycle of deflation. 

The material removed from these sites must be accounted for 
even if transported considerable distances. 

The total effects of wind action, both deflational and 
depositional are very large in this region. The extensive 
Killpecker dune field is still partially active, and is per- 
haps the most striking example; though there are extensive* areas 
of active dunes in the eastern part of Separation Flats. The 
results of wind action in other parts of the Great Divide 
basin are less pronounced though several areas of stabilized 
dunes are present (location 1, Fig. Q-l). 

In many areas thin strands of completely stabilized linear 
dunes can be seen 0 These appear on the ground as low ridges 
covered by sagebrush ( Artemesia sp.); though other plants, such 
as rabbit brush ( Chrvsothamus sp 0 )^ favor this environments 
The ridges are from two to ten feet in hieght and, 10 to 30 
feet in width. The dune ridges vary in length from tens to 
hundreds of feet. The dunes are unobtrusive on first 
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inspection, and become evident only as the observer learns to 
discriminate these forms from the bedrock vegetation associa- 
tion. 

The total effect of the distribution of these smaller, 
finer linear dunes is to give a northeast trending grain to 
certain areas of topography which might be interpreted as a 
bedrock structural grain. 

The scale of the imagery precludes final judgment as to 
how much of the photo linear aspect is due to aeolian action, 
rather than bedrock control. This investigator believes that 
wind created phenomena are far more pervasive than previously 
considered. 
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Appendix R : Tectonic analysis of southwestern Wyoming from 

ERTS-1 imagery. 

from a special report on work performed under contract 

NAS 5-21799 s 

Blackstone, D. L 0 , 1972, Tectonic analysis of south- 
western Wyoming from ERTS-1 imagery; National Technical 

Information Service, Springfield, Virginia, Report E72- 

10377, 8p. 

Examination of the positive transparency no. 1013-17300-5, 
August 5, 1972, Rock Springs Uplift area, Wyoming and those 
immediately adjacent reveals two previously unreported tectonic 
features (Figure R-l). 

1. A linear element trends about M 40° W (320°) and passes 
approximately through the common comer of Townships 

25 - 26, N., and Ranges 105-106 W., 6th Principal Meri- 
dian, Wyoming. The element passes approximately 4 miles 
northeast of the community of Farson, Sweetwater County, 
Wyoming and will be referred to by the name Farson ele- 
ment.. 

2. A linear element trending N. 65° E. (65°) and aligned 
northeastward from the town of Rock Springs, Sweet- 
water County, Wyoming, and thence across the Red 
Desert to a point near Lamont, Carbon County, Wyoming, 
will be called by the name Red Desert element. 

The two linear elements are inferred at this time to repre- 
sent the surface trace of deep seated structural elements, 
probably faults, which increase in magnitude of displacement 
with depth. The element passing near Farson is subparallel to 
the zone of normal faulting previously mapped at Tabernacle 
Butte (McGrew & Berman, 1955). The linear element continues 
southeastward into the area of the Leucite Hills volcanic field 
(Kemp & Knight, 1903). The northeast flank of the Rock Springs 
Uplift as far south as the Black Buttes station on the U.P.R.R. 
is also aligned with this trend. 

Two interpretations of the data are mades 

a. The Leucite Hills volcanic vents are associated 
with and possibly controlled by the Farson linear 
element , 
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b. The Farson structural element was active in Lara- 
mide time delimiting the northeast flank of the 
Rock Springs Uplift; and must have continued ac- 
tive until the eruption of the Leucite Hills vol- 
canic products in Pliocene time. The volcanic 
rocks are dated as 1.25 m. y. by Bradley (1961). 

The second structural element trends northeastward from 
Rock Springs, crosses the Rock Springs Uplift and extends across 
the Red Desert/Great Divide Basin to the east. The east flank 
of the Rock Springs Uplift is ruptured by numerous high-angle 
reverse faults (Schultz, 1909) with northeast strikes. It 
should be noted that previously mapped and documented struc- 
tural elements having the same structural trends exist in 
southeast Wyoming • The Como Bluff (Dunbar, 1944) to Pinto 
Creek faulted anticlinal trend (Kaabar, 1971) is a very promin- 
ent feature displacing the Precambrian basement rocks along its 
eastern extent in the Laramie Mountains. The Mullen Creek 
part of the Mullen Creek - Nash Fork shear zone (Houston, 1968) 
in the northern Medicine Bow Mountains, Wyoming, also has the 
same orientation. The rather poorly defined Wyoming lineament 
(Ransome, 1915, Blackstone, 1951) is a zone of structural dis- 
continuity that is also roughly parallel to the Rock Springs 
element . 

The Precambrian rock exposures of the northern Laramie 
Mountains terminate along a faulted syncline at Deer Creek 
southeast of Casper, Wyoming which is exactly on the north- 
eastward extension of the Rock Springs - Red Desert linear ele- 
ment. There is no evidence of structural disturbance where the 
projected alignment crosses the Ferris - Seminoe Mountains in 
central Wyoming. 

The significance of the Rock Springs - Red Desert element 
is that it further documents a system of northeast trending 
structural elements in SE Wyoming which have had a profound 
effect upon the delineation of the tectonic pattern* 

In terms of interpreting the Laramide tectonic pattern of 
Wyoming the new data reinforces the concept that the pattern is 
best compared to a mosaic of blocks; and that the blocks are 





bounded by structural elements which have a limited number of 

orientations, probably not more than six. 

Regional tectonic implications for origins of the tectonic 
patterns in Wyoming have been suggested by several investiga- 
tors (Blackstone, 1963; Malahoff and Moberly, 1968; Sales, 1968). 
The basic question to be resolved is whether the major tectonic 
elevations and depressions of the crust in the Rocky Mountain 
foreland are due to a compressive, horizontally directed stress 
field or a stress field reflecting dominantly vertical motion. 

The existence of major subparallel linear elements, extend- 
ind for tens of miles suggests a compressive stress field over 
large areas, at least in the early stages of the regional defor- 
mation • 
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Appendix S > Comparison of ERTS imagery with Skylab S-190A, 

~ S-190B, and aerial photographs for lineation 

mapping . 

from a special report under contracts NAS 5-21799 and 
Nas 9-13298 s 

Houston, R. S., Marrs, R. W., and Tomes, B. J 0 , 1974, 
Some illustrations of the advantage of improved resolu- 
tion in geologic studies ; National Technical Information 
Service, Springfield, Virginia, Report E74-10628, p. 56-67. 

A portion of the Wind River Mountains of Wyoming has been 
studied by use of ERTS images, Skylab 190A and 190B photographs, 
aircraft photographs, and on the ground. The Wind River Moun- 
tains are a northwest striking, Precambrian cored, major uplift 
located in west central Wyoming. The uplift is approximately 
100 miles long and over 30 miles wide at its greatest width. 

The uplift has the form of an assymetric anticline with steep 
overturned and faulted limb on the southwest and great gently 
dipping northwest striking hogbacks of Paleozoic and Mesozoic 
sedimentary rocks on the northeast (Berg, 1962). The Precam- 
brian rocks of the core of the structure are primarily Early 
Precambrian felsic gneisses and schists and various types of 
foliated granite cut by basalt dikes of different Precambrian 
ages. In the extreme southeast part of the uplift (Bayley, 
1965a), a major greenstone belt cuts through the range. 

The Precambrian core of the Wind River Mountains is cut 
by numerous fractures ranging from major faults and shear 
zones to minor faults and joints. These fractures have been 
mapped in a small part of the uplift (Barrus, 1968; Granger 
and others, 1971; Parker, 1972; Pearson and others, 1971; 
and Richmond, 1945) and they range from shear zones with exten- 
sive catalysis and mylonization that probably developed during 
Precambrian episodes of high-rank metamorphism to discrete 
fractures of Tertiary age. Many of the Precambrian fractures 
were filled with basaltic magma and are now expressed as basalt 
and diabase dikes that cut the gneisses and granitic terrane. 

The most striking aspect of ERTS and Skylab images of the 
Wind River Mountains is the multitude of lineations that are 
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primarily expressions of the fractures and dikes, that show on 
the images. The initial study of the Wind River lineations 
was made by Parker (1972, Appendix N) who mapped the major 
fractures of the mountains using ERTS MSS-5 and 7 images and 
compared these lineations seen on aircraft images and recog- 
nized on the ground by a field survey. 

Figures S-l and S-2 show lineations mapped from ERTS 
images and Skylab 190A photographs of the western part of the 
Wind River Mountains. In this study all lineations were mapped 
that could be recognized, rather than simply the major linea- 
tions that were mapped by Parker. Three factors affected the 
mapping -- resolution, sun angle, and sun azimuth. ERTS is 
sun synchronous so that images are ordinarily generated in the 
morning with a sun azimuth that established a bias towards 
northeast-striking linear features; whereas Skylab passes give 
different perspectives and may result in the emphasis of a 
completely different set of lineations as shown in Figure S-3. 
Rose diagrams in Figure S-3 show that there are not only more 
lineations recognized on the Skylab 190A photograph but that 
the major set strikes north-northwest rather than northeast as 
is the case for ERTS. This is a sun azimuth effect because 
the major linear features are approximately perpendicular to 
the sun azimuth in both examples. It is primarily due to the 
fact that many lineations are expressed topographically and 
sun shadowing due to topography is more pronounced at right 
angles to the direction of illumination. Sun elevation will 
also play a role in lineation mapping because low sun eleva- 
tions cause distinct shadowing where lineations are topographi- 
cally expressed. Therefore, a few more lineations might have 
been mapped at the 42° sun elevation of ERTS than at the 49° 
sun elevation of Skylab, but this cannot be quantitized because 
of the variation in azimuth. 

A judgment of the affect of resolution cannot be made 
di(rectly because of the role played by sun azimuth. For exam- 
ple, Parker's map of major mapped fractures in the Wind River 
Mountains (Appendix N) suggests that more fractures strike 
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Figure S— 



.Comparison of linear features mapped using ERTS images and Skylab S- 1 90 B photographs 
of a part of the Wind River Mountains showing sun-azimuth bias. 
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Fiqure S-l. Ma P showing all lineations that can 
be detected on ERTS images of the 
Wind River Mountains. 
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Figure S-2. Map showing all lineations that can be detected on Skylab 
S-190A photography of the Wind River Mountains. 
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northwest than northeast. If this is true throughout the range, 
the sun azimuth of Skylab would tend to bring out lineations 
that result from these fractures better than that of ERTS so 
that more lineations should show on Skylab photographs regard- 
less of resolution. However, if we select a lineation direction 
such as east-west that should not be biased by sun azimuth the 
number of lineations mapped on a Skylab photograph is about 
three times that of ERTS for the area shown in Figure S-3. 

This cannot be taken as a quantitative measure of the exact 
resolution difference for lineation measurements because the 
total number of lineations measured by use of Skylab photo- 
graphs is only an order of magnitude greater than ERTS for 
this area (Fig. S-3), but it is clear, nonetheless, that more 
lineations can be recognized on the higher resolution Skylab 
190B photograph. 

If we compare known faults (as mapped by field geologists) 
with linears mapped using Skylab and ERTS, we find a similar 
correspondence. In the northwestern part of the Wind River 
Mountains Richmond (1945) mapped a series of faults and shear 
zones and this mapped area is covered by both Skylab and ERTS. 
Results of study of linears from this area are shown belowi 


Total number of linears 

Total length of linears 

Total number of linears that 
correspond to mapped faults 
and shear zones 

Total length of corresponding 
linears 

Total number of corresponding 
linears unbiased by sun azimuth 


Skylab (190A) ERTS 

202 104 

600 kilometers 520 

(kilometers) 

17 13 

95 kilometers 95 

(kilometers) . 
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Appendix T s Vegetation analysis in the Laramie Basin, Wyoming 
from ERTS-1 imagery. 

from a special report on research conducted under contract 
NAS 5-21799: 

Evans, M 0 A., and Redfem, R. F. , 1973, Vegetation 
analysis in the Laramie Basin, Wyoming from ERTS-1 imagery; 
National Technical Information Service, Springfield, Vir- 
ginia, Report E74-10145, 54p. 

A vegetation map of the Laramie Basin test area was con- 
structed from ERTS image interpretations in order to test the 
utility of the ERTS imagery for mapping and analysis of range 
vegetation. Image 1029-17148, bands 5 and 7, was used in the 
interpretation. The image interpretation was supplemented by 
an isodensity contour map of the band 5 image. 

Some vegetation types can be discriminated from the red 
band alone. Heavily vegetated floodplains are easily separated 
from drier areas supporting other vegetation types. Flood- 
plain vegetation shows a variety of contrasting radiance pat- 
terns on both bands 5 and 7 . These patterns are readily inter- 
preted by comparison of density contour maps derived from 
ERTS bands 5 and 7. Ligh- toned areas surrounded by floodplain 
vegetation may be elevated grassland areas or simply alkaline 
flats devoid of heavy vegetation cover. Alkaline areas show 
up distinctly on band 5 but are easily confused with dense 
mesic vegetation on the infrared bands. 

The Laramie Basin is predominantly a shortgrass prairie. 

This grassland exhibits a mottled appearance on the ERTS imagery. 
In some areas slightly alkaline conditions produce light-toned 
regions on the band 5 imagery. Straight-line tone changes often 
represent fence lines and corresponding differences in grazing 
intensities. Some of the tonal contrasts may reflect differ- 
ences in species composition of the grassland, especially 
where grazing pressure affects species composition and vege- 
tative biomass. 

Reflectance characteristics of various range vegetation 
types were observed qualitatively. Plant absorption in the 
visible spectrum is due primarily to quantum absorption which 
is a function of the density of pigments. Thus, absorption in 
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the green and blue bands indicates the amount of pigments pre- 
sent, especially chlorophyll. The chlorophyll content, in 
turn, can be related to vegetation type, cover, biomass, water 
availability, leaf area index, available nutrients, etc. In- 
frared reflectance, which is substantial in many species, is 
closely related to development of the leaf structure, particu- 
larly in the mesophyll layer, as well as to leaf area index and 
cover, (Gates et al, 1965; Hoffer and Johannsen, 1969; Gates, 
1970; and Meyers, 1970). Thus, broad-leaf vegetation in mesic 
environments, with profuse and well-developed leaves, should 
exhibit considerable absorption of the visible energy and have 
high infrared reflectance. This holds true for the Laramie 
Basin. Drier vegetation shows less visible absorption and 
less infrared reflectance. Xeric plains vegetation, with its ^ 
poorer leaf structure development and sparse cover, shows little 
difference in visible and infrared reflectance. 

The effect of water on leaf structure and quantity of 
vegetation is apparent when comparing crops in dry and irri- 
gated fields and dry plains with wet hay meadows. For example, 
irrigated crops appear very dark in the visible ERTS-1 bands 
and extremely light in the infrared, whereas the effect is 
less pronounced in dry-land crops. Differences in crop species 
do not affect reflectance variations as severely as does water 
availability. Similarly, the leaf structure of the wet meadow 
vegetation is well developed and with a large leaf area index, 
thus, much visible radiation is absorbed and infrared is re- 
flected. Dry grasslands show little or no difference between 

infrared and visible reflectances. 

The ERTS images are definitely of value in photogeographi- 
cal studies. Readily discriminated with a large degree of 
confidence are deciduous forests, coniferous forests, and 
floodplain vegetation, irrigated farming, dry-land farming, 
alkaline areas, and grassland-shrubland. 

The two ERTS-1 bands that appear most useful are 5 and 7, 

Band 4 is similar to band 5 with regard to vegetation reflec- 
tance, whereas band 6, being a combination of deep red and 
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very near infrared, represents a frequency range in which many 
plants show high absorption in the red and high reflectance 
in the infrared. As such, band 6 is generally no more effec- 
tive in vegetation analysis than an average of bands 5 and 7. 
However, this does not mean that bands 4 and 6 are never im- 
portant in interpretation. Indeed, their value may be signifi- 
cant in quantitative vegetation analysis. 

There are no coniferous forests in the area of the prelim- 
inary vegetation map. However, in nearby mountain areas, 
conifers are readily mapped from the imagery because they are 
relatively dark in both bands 5 and 7. Clearcut and burned- 
over areas can be identified and studied quantitatively. Cloud 
cover is common in high mountain areas and often obscures 
large areas of coniferous forest. A small test area lying 
immediately south and east of the main Laramie Plains test 
area was chosen to evaluate the ERTS data for mapping of vege- 
tation in mountainous areas. The area was mapped both by 
interpretation of aerial photographs and of ERTS 1 imagery. 
Again, isodensity contour maps were used as an interpretation 
aid. 

A field check of the vegetation zones mapped from the 
ERTS imagery was conducted in late June and early July 1973. 

Both the Laramie Basin and Laramie Range study areas were 
field checked. In both areas a number of each of the mapped 
zones were checked, and a species account was made. The dom- 
inant species in each area were determined and a list was made 
of all species encountered in each area. 

The ERTS imagery is seen to provide resolution approaching 
theoretical limits through the clear and dry Wyoming atmosphere. 
Consequently, comparatively small stands of vegetation can be 
identified under optimum conditions: they must be somewhat 

wider than one resolution element. A minimum sized stand 
should be on the order of .16km by .16km, a size equal to 
four resolution elements, to obtain a value approaching the 
actual radiance (less atmospheric attenuation). The larger the 
stand, the less the worry with edge effects and scattering into 
the instantaneous field of view, and extensive communities of 
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the major vegetation types (conifer, brush, grassland, etc.) 
should be delimitable from ERTS data.. Atmospheric scattering, 
however, does not appear to be an important problem in the 
Laramie Basin during optimum imaging weather. 

Vegetation identification by ERTS imagery is often limited 
by factors such as sun angle, topography, complexity of the 
vegetated landscape, and substrate. A scene-radiance which 
would be correlated to brush on light soils with a tolerable 
degree of accuracy could not be so correlated on darker soils, 
where grassland might give the same radiance. Vegetation on 
slopes sheltered from full sunlight will reflect levels of 
radiation dissimilar to corresponding vegetation on slopes 
bathed in full sunlight. Changes in sun angle can alter re- 
flectance as well as radiance due to the effects of shadows 
and light incidence angle. In view of such dissimilarities, 
the concept of "signature" in ecological remote sensing is 
not useful in a strict sense. Thus automatic mapping and 
monitoring of "natural" vegetation is obviously far more 
difficult than that of crops, which has had its own share of 
problems (Hoffer and Goodrick, 1971). 

There are several possible solutions to the problems of 
satellite-based remote sensing of natural vegetation. Clus- 
tering techniques (Smedes et al«, 1971), which group like 
radiances for interpretation, can be useful but not sufficient 
in themselves for identification as indicated by the present 
study. With computer analysis of several radiation bands, it 
is possible to effect various mathematical transforms of the 
radiance data (Smedes, 1971 and Kridgler et al., 1970). Such 
transformations promise to be particularly useful in molli- 
fying effects of atmosphere and topography. However, if it is 
deemed useful to map and monitor vegetation from orbit, it 
appears clear that much low-level and ground— truth data will 
be needed for vegetation maps of any accuracy • Apparently it 
will be easier and more accurate to map vegetation by predomin- 
antly manual methods rather than automatically if maps amenable 
to botanists are to be produced. In addition to proximate data. 
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models which, comprehend factors such as topography, sun angle, 
and substrate are required to automatically monitor vegetation „ 
Of course, the volume of augmentative data is a function of the 
accuracy desired. It appears possible to map broad vegetation 
types such as conifers, broad-leafed trees and shrubs, brush- 
land, grassland, desert, etc., with present ERTS imagery alone. 
Increasing resolution will eventually require more sophisticated 
imaging coupled with sophisticated techniques of analysis, in- 
cluding modeling. 

A rough inverse correlation exists between radiance types 
and the amount of vegetative biomass present. Pearson and 
Miller (1972), working in the shortgrass prairie of the Pawnee 
National Grassland, have found good correlation between the 
ratio of infrared to red reflectance and green biomass. With 
a proper model, the estimation of biomass over large areas 
appears possible. 

Certain wavelengths seem to be indicators of potential 
productivity in vegetation. For example, Odum (1971, p. 62) 
suggests that the ' v yel low-green" index of Margalef, which is 
inversely correlated with the ratio of productivity to res- 
piration, might be used as a remote sensing tool. Myers (1970, 
p. 264) suggests that reflectance at about .54 um is a function 
of the chlorophyll concentration and that increased reflec- 
tance of this wavelength possibly indicates nutrient deficien- 
cies. In addition such things as the ratio of reflectance in 
a near infrared band to reflectance in a band on the red might 
indicate the potential rate of photosynthesis. Also, remote 
sensing techniques can measure absorption of photosynthetic 
wavelengths by vegetation, especially when problems of atmos- 
pheric scattering are solved. If suitable remote sensing 
indices of ecological factors can be found, they can be aug- 
mented by meteorological and other data and possibly used to 
extimate productivity and respiration. 

A difficult but important ecological parameter to measure 
is evapotranspiration. Measurement of this parameter can be 
apprehended by energy budget studies. If such studies can be 
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done for given elements of terrain from space platforms, this 
parameter, an important factor in water cycling as well as an 
indicator of net primary productivity (Rosenzweig, 1968), might 
be estimated and applied to vegetation analysis. 
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Appendix U i Identification of irrigated crops. 

from a special report on research conducted for contract 

NAS 5-21799: . . 

Evans, M. A., and Marrs, R. W., 1974, Identification 
of irrigated crop types from ERTS-1 density contour maps 
and color infrared photography; National Technical Infor- 
mation Service, Springfield, Virginia, Report E74-10312, 
12p. 

The use of ERTS-1 imagery for identification and mapping 
of irrigated and dry-land crops was investigated in an agricul- 
tural area north of Wheatland, Wyoming. Initial attempts to 
identify crop types by direct photointerpretation were only 
marginally successful. On the other hand, a trained interpre- 
ter could readily identify the various crops on 1:120, 000- 
scale color infrared photography of that same area. The con- 
trast apparent from the color infrared aerial photography and 
the ability of the interpreter to delineate the field patterns 
on the ERTS imagery suggest that image-enhancement procedures 
might improve the success of crop-identification from ERTS. 

The lack of textural information and detail on the ERTS imagery 
might be compensated by a more accurate estimate of the rela- 
tive brightness of the fields on the different ERTS bands. 
Consequently we elected to construct density contour maps of 
the test area from the ERTS imagery, in an effort to enhance 
the critical contrasts. 

\ The test area is a mixed dry-land/irrigated agricultural 
^irea approximately 2- by 2 l/4-mi. in size and located about 
4 miles northeast of Wheatland, Wyoming. The area comprises 
both irrigated and dry-land farms, rangeland, and heavily- 
vegetated floodplains. Each of these three major zones is 
represented by a distinct tone, texture, and pattern which can 
be readily identified on the ERTS-1 imagery. Further subdi- 
vision of the zones is more difficult. 

Fourteen crop classes were recognized from the August, 

1973 color infrared photography of the test area. They include 
com, sugar beets, alfalfa, grains (oats and barley), grain 
used as a nurse crop for alfalfa, grass/hay and pasture, 
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beans, plowed ground, natural vegetation, natural vegetation 
supported by irrigation-water runoff, water (ponds and reser- 
voirs), uncultivated floodplain, cultivated floodplain (al- 
falfa) and farm yards. 

Density contour maps were constructed from the ERTS imagery 
to give an enlarged presentation of the test area along with 
absolute film density information. The' large-scale (1:19,760) 
density map permits accurate delineation of the contrasting 
fields, while the density contours allow the interpreter a 
means of making objective decisions about contrasts and similar- 
ities of various areas. He can also use calibrated density 
contour maps as a means of comparing various bands and time- 
sequential images. The interpreter must use these contour 
maps in much the same way as he would a low-altitude photo- 
graph, keeping in mind the mechanical variables as well as 
the natural variables (irrigation, climatology, etc.). 

In this study, three time-sequential ERTS images were 
selected on the basis of image quality and time distribution. 
The red-band, 9-in. positive transparencies were chosen for 
time-comparison density slices. The instrument used was a 
Joyce Loebel/Tech Ops isodensitracer. 

For control, an outline map of the crop types was con- 
structed from the aerial photograph. The crop types were 
positively identified from field checks and were the same 
crops throughout the summer. The settings on the micro- 
densitometer were the same for each density contour map. 

Figures U-l, U-2, and U-3 are black-and-white duplicates 
of the red-band isodensitracings (density contour maps). 

In original form, the various density patterns are displayed 
in color. A crop map, interpreted from a 1:120,000 color 
| infrared aerial photograph of the area has been superimposed 

| on each of the isodensitracings. 

The density contour map shown in Figure U-l was taken 

from ERTS-1 image 1334-17134-5, 22 June 1973. This image was 

taken early in the growing season and the natural vegetation 

areas show up dark, because of early growth of grasses and 

l sagebrush. The com, beans and beets are still young and 
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show as darker areas suggesting various amounts of vegetative 
cover • 

The density contour map of Figure U-2 taken from ERTS-1 
image 1352-17132-5, 10 July, shows the natural vegetation as 
lighter than in June, indicating dryer conditions or grazing 
pressure. The beets are darker and can be distinguished 
from the com. The alfalfa and grain are generally darker 
than in June. 

The isodens it racing shown in Figure U-3, from ERTS-1 
image 1388-17125-5, 15 August, shows the com and beets as 
darker than before. The alfalfa in fields 4, 7, and 8 is 
much lighter, after having been harvested. Similar observa- 
tions drawn from other images and density contour maps serve 
as correlative data for interpretations regarding the vege- 
tativetype and condition • 

An additional application of the isodensity map lies in 
the enhancement of boundary and mix problems. Examination of 
the three density contour maps (Figures U-l, U-2, and U-3) 
shows that large, evenly vegetated fields are quite distinct, 
but, as the fields become progressively smaller or less even, 
or as contrast decreases , the crops become more and more 
difficult to map and identify. The density contour maps give 
the interpreter a good indication of the magnitude of this un 
certainty. 
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Appendix V » ERTS-1 imagery as a tool for mapping major vege- 
tation zones in Wyoming o 

a previously unpublished report by? 

M. A. Evans and R. W. Marrs on research performed 

under contract NAS 5-21799. 

A false-color infrared ERTS-1 mosaic of Wyoming was used 
as an interpretation base for compiling general vegetation 
maps of Wyoming. Individual ERTS-1 images were then selected 
for more detailed interpretation of vegetative subzones in 
various kinds of terrain., 

The earliest report concerning plant zones of Wyoming is 
by Nelson (1896). Later, Cary (1917) reported on the life 
zones of Wyoming while Kirkwood (1922) and Rydberg (1928) dis- 
cussed forest distribution and life zones of the Rocky Moun- 
tain region. More recently, Costello (1954), Daubenmire (1943), 
and Porter (1962) have discussed vegetation zones of the Rocky 
Mountain region and Wyoming. Yet, there are still no compre- 
hensive vegetation maps for Wyoming. The chief value of 
these previous works lies in the descriptions of the various 
vegetation zones observed. These descriptions provide much 
of the information necessary as background for mapping large 
areas with ERTS imagery. Individual locations where the vege- 
tation and geology are known can be pinpointed on an ERTS-1 
image and correlated with available field data to identify 
typical relationships between vegetative communities and spectral 
response. The texture, hue and tone of each zone has definite 
characteristics. Some zone boundaries contrast sharply and 
are easy to define. Other zones blend together and are diffi- 
cult or impossible to distinguish. Experience and personal 
knowledge of the area and zone are a considerable asset in map- 
ping these indistinct boundaries. 

This report presents maps of the vegetation zones of Wyo- 
ming compiled from ERTS-1 images and illustrates the utility 
of ERTS-1 imagery in regional vegetation mapping. The project 
summarized in this report, including field work, required about 
two months total time. Although mapping vegetation from ERTS 
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imagery has advantages in speed and efficiency , some natural 
detail is missed. The resolution of the imagery is limited 
to features 70 meters or greater in diameter. Thus, to 
purists of the natural sciences, ERTS image interpretations 
may he practically useless. But, to many, ERTS offers a 
unique apportunity to obtain up-to-date general information 
previously difficult or impossible to obtain. 

The Wyoming index map (Figure V-l) shows major drainage 
systems, water bodies, county boundarxes and major cities 
for use as reference with the Wyoming vegetation maps (Figs. 

V-2, V-3, V-4 and V-5). The imagery used in compiling the 
regional vegetation maps was a false-color composite mosaic 
of Wyoming (see Fig. 1 this report) prepared by General 
Electric Photographic Engineering Laboratories. The vegeta- 
tion zones were chosen to correspond to those defined by 
Porter (1962) and currently in use by Wyoming botanists 
(Dorn, 1972). The maps compiled from the ERTS image inter- 
pretations were partially field checked during summer and 
fall, 1973. 

Vegetation zones of Wyoming are primarily controlled 
by altitude. The major vegetation zones descriminable 
from the ERTS-1 color composite mosaic ares 1) alpine, 

2) timbered mountain slopes, 3) foothills scrub, 4) grass- 
land including eastern plains and foothills (Great Plains) 
and interior grassland plains, 5) desert and basins, and 
6) river bottoms. The major agricultural areas of Wyoming 
were also mapped. Major lakes are shown by double outlines. 

A list of the species found in each zone is given by 
Porter (1962). 

The alpine zone of Wyoming includes the high mountain 
peaks and ridges above timberline, usually above 10,500 
feet in elevation. Typical terrain includes rocky summits, 
scree slopes, lakes, alpine tundra, meadows, snowbanks, 
and in some' places , glaciers • On the ERTS imagery , white 
to bluish areas are characteristic of the snow fields and 
bare rocky summits. Lakes appear as small black dots and 
alpine tundra and meadbws appear light pink. 
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Vegetation of the alpine zone is sparse, consisting 
chiefly of grasses and low-growing herbaceous flowering 
plants that flourish briefly in the short growing season. 

Some dwarf or scrubby woody shrubs are present. 

Timbered mountain slopes (Fig. V-3) extend from the 
alpine zone to an elevation of about 6,000 feet. Conifers 
are the chief visible component of the vegetation community. 
Aspen and grassy meadows are interspersed throughout 
this zone. On the ERTS imagery, zone 2 appears dark red 
to black in color with small areas of bright red and pink 
representing meadows and aspen groves. Figures V-6 and 
V-8 are ERTS images illustrating vegetation zones 1 and 2. 
Figure V-6 corresponds to area IV, shown in Figure V-2, 
while Figure V-8 corresponds to area I of Figure V-2 . 

Figures V-7, V-9, and V-10 are vegetation maps compiled 
from interpretations of the areas shown in Figures V-6 and 
V-8. 

The Columbia Plateau element is found only in extreme 
northwestern Wyoming, primarily in the Yellowstone National 
Park and Jackson Hole area* Little data was available 
on this element. Consequently, it is mapped as open parks 

and meadows in Figure V-7 • 

The Black Hills element of northeastern Wyoming is 
mapped as timbered mountain slopes. On the ERTS false- 
color mosaic, this area has a mottled red to dark red 
appearance. The vegetation here is more related to the 
midwestern deciduous flora, than to the Rocky Mountain 
flora. Figure V-ll is an ERTS image which shows the Black 
Hills element and the zone. Figure V-12 is a detailed 
zone map of this region. 

Figure V-3 shows zone 3 (foothills scrub), zone 4 
(grassland) and zone 5 (deserts and basins) , Zone 3, 
foothills scrub, is found along the lower slopes of 
mountains.’ The vegetation of this zone is mostly shrubs 
and small trees with some grasses. The foothills scrub 
zone is often rocky and arid. Mountain mohogany, sagebrush. 



Figure V-6. ERTS-1 image 1015-17404-4,5,7, 7 August 1972, showing the 
Yellowstone Park and Columbia Plateau areas of Wyoming, 
Idaho, and Montana. 
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Figure V-7 . A detailed map of Zones I and 2 taken from the original of Figure V-7 . 





Figure V-8. 


ERTS-1 image 1013-17294-4,5,7, 5 August 1972, showing the 
Wind River Mountains and Wind River Basin of central 
Wyoming. 
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Figure V-12. A detailed map of Zone I ( BlacK Hills Element of Timbered Mountain 
Slopes) and Zones 3, 4, and 6 taKen from the original of Figure V-ll. 




scattered limber pine, and juniper are common shrubs and 
trees. On the ERTS mosaic this zone appears brown to 
reddish where more moisture is available. 

Zone 4, grassland, is characteristic of the rolling 
plains and hills at lower elevations (3,000 to 6,000 feet). 

The grasslands of Wyoming can be divided into two major 
types, l) eastern plains and foothills, and 2) interior 
grassland plains. In Figure V-3, only the eastern plains 
and foothills type are mapped as grassland. This general 
type is characterized by "relatively tall grasses and 
forbs" and a truly grassland type of flora. But this 
flora grades into a more arid basin northward (Powder 
River Basin). Much of the grassland in eastern Wyoming 
has been converted to agricultural uses (Fig. V-5). On the 
ERTS imagery, the eastern plains and foothills grasslands 
appear pale pink to green and are limited to a strip east of 
the Laramie Range and southwest of the Black Hills. 

The interior grasslands are closely related to desert 
and basin flora. This type is a shortgrass prairie domi- 
nated by blue grama, bluegrass, and rabbit brush with some 
sagebrush present. The type is defined on ERTS as a mottled, 

off-white to light -brown color. 

Zone 5, desert and basin, includes flat, undulating 
areas at elevations from 4,000 to 6,000 feet. These areas are 
very arid and often alkaline. Dry streambeds and playa 
lakes are also very common. The sparse vegetation is mostly 
xerophytic and saline tolerant. On the ERTS mosaic, these 
vast areas appear cream to white in color . In the Great 
Divide Basin area. Figure V-13, (Area III, Fig. V-3) the 
variety of tones and colors is due to change in soil and 
rock types and changes in density of vegetation. The blue 
-to blue-green color represents areas covered by sagebrush. 

Figure V-4 (Map III) shows zone 6, major river bottoms 
and irrigated croplands. This zone is limited by the effects 
of water. Natural river bottoms are covered by lush grasslands 
and dense growths of cottonwood and willows . In very arid 
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areas there is no mappable river-bottom flora. Where 
topography permits, the water is often used to irrigate 
native hay meadows, alfalfa, and other crops. On the ERTS 
mosaic, zone 6 appears as ribbons or bands of light red 
to red. 

Figure V-5 (Map XV) shows the major agriculture areas 
of Wyoming. A previous paper (Evans and Marrs, 1974) discusses 
the application of ERTS-1 imagery to mapping of irrigated 
croplands • 

This mapping program demonstrates that ERTS-1 images 
are of considerable value in vegetation mapping. Tonal 
changes, color, brightness, hue, and texture of the ERTS 
mosaic can be related to vegetative zones and patterns. 

Often the ERTS imagery defines detail which can not be 
shown in the general mapping units used here. To map such 
detail, one must first conduct more intense field checking 
and make careful studies of smaller areas. He may then 
be able to expand the techniques to larger areas. Also, 
the ideal mapping base would be a mosaic compiled from 
images obtained on consecutive days or a series of such 
mosaics of the area representing different periods throughout 
the year or growing season. 

The desert and basin and grasslands zones are the most 
difficult to map accurately due to the great variety of 
tones and changes in color. For the purpose of this program, 
the maps for these areas were generalized as described 
above. More detailed maps can be produced from the ERTS-1 
imagery, but such mapping requires the careful analysis 
of the spectral reflectance values recorded by ERTS and 
correlation with a ground control such as that used by 
Gordon (Appendix W, this report). 
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Appendix W : The use of ERTS data for range mapping and 

estimation of above-ground, green biomass. 

from a report on research performed under contracts 
NAS 5-21799 and NAS 9-13298. 

Gordon, R. C. 1974, Multispectral imagery for 
range vegetation type mapping and estimation of above- 
ground green biomass. National Technical Information 
Service, Springfield, Virginia, 45 p. 

Objectives 

The objectives of this research were: (l) to determine 
whether or not the multispectral imagery available from 
the ERTS-1 satellite was useful for mapping the rangeland 
resource, and (2) to estimate above-ground green biomass 
by relating it to total scene reflectance, thereby diversi- 
fying the applications of remotely sensed data for use in 
the field by the resource manager. 

Neither of these goals were set with the idea of 
eliminating field work but, to assist it by determining 
changes in the rangeland and its above-ground green 
biomass on a regular basis. The development and advance- 
ment of remote sensing techniques for range management will 
allow the conservationist to locate areas which require 
his attention. 

Study Area 

The study was conducted on 68,000 acres of rangeland 
in the southeastern portion of the Red Desert west of Baggs, 
Wyoming . Extending through the area are the southern 
structural rim of the Washakie Basin and a complex system 
of east-west oriented faults, synclines and anticlines 
(Gibbens, 1972). The irregular topography extends from 
6,200 to 7,200 feet above sea level. 

Shales, sandstone or alluvial mixtures of the two 
make up the major soil parent materials. The often cal- 
careous , brown to pale-brown loams make up the majority of 
the soils described by the Soil Conservation Service, 

These loams have developed mainly from the Browns Park 
formation and the montmorillanitic clay soils from resid- 
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ual or transported shales • Both are neutral to very 
basic in reaction (Gibbens, 1972). 

The climate consists of cold winters, broken by a 
short hot summer. Precipitation ranges from about 8 to 
11 inches per annum with the majority of winter moisture 
in the form of snow. Gibbens (1972) states that the grow- 
ing season is approximately 96 days. 

The vegetation types are usually dominated by big 
sagebrush ( Artemisia tridentata Nutt.) or Utah juniper 
( Juniperus osteosperma (Torr.) Little). These provide 
good cover for antelope, mule deer, and sage grouse 
populations . 

Under the Bureau of Land Management supervision, the 
study area was divided into a four pasture rest-rotation 
unit, a two pasture deferred unit and one season-long 
grazing unit. Summer use of the pastures was confined 
to cattle, with sheep grazing the four-pasture rest-rotation 
unit during the winter. 

Procedures 

U. S. Geological Survey 1 s 24, 000-scale topographic 
maps were used as a mapping base. Information was trans- 
ferred, using a Bausch and Lomb ZT-4 zoom transfer scope, 
from colour and colour infrared aerial photography obtained 
from NASA flights 213 and 248 (Figs. W-la and W-lb) . The 
scales of these images were 1:43,400 and 1:104, 500. Accuracy 
of the vegetation type lines was determined by field 
reconaissance • 

The Skylab vegetation maps were prepared from colour 
photographs number 104-S190A and 35-S190B Track 48, taken 
on August 4, 1973, enlarged on the aforementioned zoom 
transfer to a scale of 1:250,000. This map was then enlarged 
with a Saltzman Projector system to 1:24,000 for direct 
comparison to the* base map (Figs, W-2a and W-2b) • Skylab 
images S190A and S190B, nine-inch format, were at scales 
of 1:712,917 and 1:477,979 respectively. 
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The third set of vegetation type maps was prepared 
from ERTS image number 1318-17253, band 5, taken on June 
6, 1973 (Figs. W-3a and W-3b) . This image was enlarged 
photographically to about 1*250,000, the information mapped 
directly on the print, and transferred to the Saltzman 
Projector for enlargement to 1*24,000. Band 5 was chosen 
over the other bands and the colour composites because of 
its high contrast and the difficulties encountered in making 
colour enlargements • Both the ERTS map and the Skylab 
range type maps covered the entire 68,000 acre study area. 

Vegetation transects were established in the study 
area to adequately type the vegetation. Crown cover and 
basal cover were obtained for shrubs and for grasses and 
forbs respectively on these transects using the same pro- 
cedure followed by Gibbens (1972). Duplication of this 
method allowed direct comparison of data. 

Vegetation types were named for the dominant species 
in the area of the transect and other transects with similar 
composition and cover characteristics were grouped to 
form subtypes of a major type. The major range vegetation 
types and subtypes were designated on all three vegetation 
maps . 

Prior to the first ERTS overpass of the 1973 study 
season, 10, two- by three-meter reference plots were es- 
tablished both in the Poison Buttes pasture and 10, in 
Pasture C, of the rotation pasture system. Meter readings 
of relative reflectance were obtained from each of these 
20 plots during the five 1973 summer overpasses using 
a portable photometer manufactured by Science and Mechanics. 

In addition to the 20 reference plots, 350 plots were 
established within which above-ground green biomass was 
determined. Above-ground green biomass included all 
green herbaceous material and leaves of shrubs but not 
perennial branches . The Pasture C and Poison Buttes bio- 
mass sites each had 175 of these plots. 
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Relative reflectance measurements were obtained for 
each plot under conditions similar to those occurring 
during the overpass and using the same techniques <> Meas- 
urements were taken as soon after the overpass as possible. 

Above-ground green biomass was sampled after each 
overpass using a double sampling technique described by 
Pechanec and Pickford (1937). After the July 30, 1973, 
overpass 50 plots were clipped in entirety to have an accurate 
above-ground green biomass weight for comparison with the 
estimates made on July 31, 1974. 

Soil moisture was determined by collecting samples 
of the top inch of soil during each overpass in the 1973 
field season and during the 1974 field season, when the 
comparative transects were being measured for relative 
reflectance. 

Relative reflectance information was obtained with a 
calibrated photometer. The photometer was calibrated 
using the procedure followed by Raines and Lee (1974). 
Reference grey cards of known reflectance were then measured 
and the values were used as a standard for the field meter 
readings. Meter readings recorded in the field were 
converted to percent relative reflectance by plotting 
them on the graph of reflectance versus meter reading. 


Results 

Mapping with the ERTS imagery, although hampered by 
colour bias, cloud cover, loss of detail through enlargement 
and the resultant poof location association caused by the 
latter, was successful. Both ERTS and the Skylab images 
provided useful maps of the range vegetation types. The 
ERTS-1 map contained representatives of all the vegetation 
types except one and it was the smallest type • The type 
not included was a seeded patch of Agropyron cristatum . 

Type 0. The image used was band 5 (red, .6 to .7 ) and 

was taken in June, 1973. The map produced from this 
image would serve as a useful survey base for the resource 

manager • 


Techniques for estimation of above-ground green bio- 
mass have not been confirmed due to the mechanical failure 
of the photometer o However, graphs of log transmittance 
versus above-ground green biomass (Fig. W-4) indicate two 
things. First, it appears possible, for the Poison Buttes 
site, to estimate above-ground green biomass by relating 
it to log transmittance. Secondly, for Pasture C, it may 
be possible to estimate the green biomass produced on this 
location by relating it directly to an area whose above- 
ground green biomass can be estimated from the ERTS imagery. 
These two conclusions may be tested by relating the 1974 
summer imagery with the 1974 production data. 
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Figure W-la Range vegetation map compiled from 

low level aerial photography supplemented 
with ground reconnaissance 
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Figure W-lb Range vegetation map compiled from 

low level aerial photography supplemented 
with ground reconnaissance 
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FigureW-2aRonge vegetation map interpreted 
from SKYLAB photographs 
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Figure W-2tiRange vegetation map interpreted 
from SKYLAB photographs 
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FigureW-3aRange vegetation map developed 
from ERTS-I imagery 
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from ERTS-I imagery 
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Appendix X : Mapping of seleniferous soil units. 

a summary of research for contract NAS 9-13298: 

Kolm, K. E., 1975, The Location and Mapping of 

Seleniferous Soil Units in the Wasatch Fm. , Powder 

River Basin, Wyoming, using Sky lab Photographs; 

Dept, of Geology, University of Wyoming, Special 

Report, (in progress). 

/ 

Within the Wasatch Fm. , Powder River Basin, Wyoming, 
sandy soils contain considerable amounts of selenium and 
support the toxic plant species of Astragalus bisulcatus 
and Oxytropis lambertii . The principal objective of this 
investigation was to evaluate the success of Skylab photo- 
graphy in preparing maps of these toxic vegetation communities 
and their associated soils. 

In order to accurately map from Skylab photographs, the 
toxic vegetation community found on the sandstones and sand- 
stone-derived soils in the Wasatch Fm. , it was first necessary 
to conduct a reconnaissance study of an area known to contain 
these plants. The reference area chosen is located just north 
of the site of AMAX* s isouth Belle Ayr Mine. 

Field work conducted at the AMAX site resulted in l) 
identification of toxic plant types present and their dis- 
tribution, 2) soil and rock type identification and their re- 
lationship to supported plant communities, and 3) a soils and 
toxic vegetation map prepared from low-altitude aerial photo- 
graphy. 

It was observed that seleniferous plants do not display 
a heavy population density, and secondary non-seleniferous 
vegetation overlapped both the toxic and not-toxic soils. 
Therefore, it was impossible to map these toxic plants using 
Skylab photographs because objects of this size are unresol- 
vable. However, a very high correlation between seleniferous 
plants and sandy soils was observed. Using the techniques 
of density analysis and photographic enlargement, a toxic 
soils map was prepared for the area chosen at the AMAX mine 
site using Skylab S190B color photographs. Track 59, Pass 28 „ 
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A regional map of these seleniferous soils in the Wasatch 
Fm. , Powder River Basin, Wyoming, is now being prepared and 
will be field checked. A similar toxic soils map, using ERTS 
imagery, will not be attempted due to the seleniferous soil 
lenses being too small to be resolved at the ERTS scale. 
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Appendix Y : Programs employed in Q-mode factor analysis of ERTS-1 multispectral data. 


DIMENSION DATA < 3 » 10 ) »C ( 55 ) *2(3*10) »S IG( 10 ) , AVGU0 » .RI100 ) 

DIMENSION V(10tlO)tE(lOI,TV<5n*H(10i.FFUO)#B(10) 

DIMENSION F( 3,10) »F$TAR(3,10) ,VSTAR( 10,10) ,C0NTRL(4) 

DIMENSION INDEX! 10) ,0(55) 

10 FORMAT (30F1.0) 

READ (105,10) { (DATAU ,J) ,J*1, 10) ,1=1,3) 

C0NTRL ( 1 ) =2. 0 
CONTRL ( 4 ) =• 995 
DO 2 1=1,2 
DO 3 J= 1 , 2 
CONTRL ( 2 ) = I 
CONTRL ( 3 ) = J 

CALL SQFAC (CONTRL, 3, 10, 10, 3, DAT A, AVG,SIG,Z,C,E, V,F,VSTAR,FSTAR,NC, 
1 TV, H, FF, B, IER,D, INDEX, R) 

3 CONTINUE 

2 CONTINUE 

STOP 
END 
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SUBROUTINE SQFAC { CONTRL.NVAR, NLOC , N, NSEC, DATA, AVG, SIG, Z *C, E, V, F, 
l VSTAR* FSTAR ,NC , TV, H,FF ,B» I ER,D» INDEX, R) 

MnnF-1 FOR R-MODE, MODE=2 FOR Q-MODE FACTOR ANALYSIS 
MTYPE*1 FOR CORRELATION MATRIX, MTYPE=2 FOR COSINE MATRIX 
DATA=NVAR X NLOC DATA MATRIX 
NVAR=NUMBER OF VARIABLES 
NLOC^NUMBER OF LOCATIONS 

C* NXN COSINE OR CORRELATION MATRIX STORED AS A SECTOR ARRAY 
WITH C( 1 ) = ( It U ELEMENT, C(2)*I2,1) ELEMENT, C(3)=(2,2) ELEMENT 
C(4)=<3,1) ELEMENT, ETC. 

M*N*(N+l>/2 

N2=N**2 

NSEC-NLOC IF N*NVAR, NSEC=NVAR IF N*NLOC 

ZeNVAR^X^NLOC^NORMAL I ZED DATA MATRIX. COMPONENTS OF Z-COMPONENTS 
OF DATA MINUS AVG (FOR R-MODE) DIVIDED BY SIG 

SIG=N COMPONENT VECTOR OF STD DEV (FOR R-MODE) OR RMS LENGTH (FOR 

AVG°N E COMPONENT VECTOR OF AVERAGES FOR R-MODE AND ZERO FOR Q-MODE 
V-NXN WORK MATRIX (EIGENVECTORS) 

E=N COMPONENT WORK VECTOR (EIGENVALUES) . 

TV.H.F.B.IER ARE DEFINED BY COMMENTS IN SUBR VARMX 

Is DE=F LEFT MATRIX IS TO BE ROTATED AND RIGHT MATRIX USEC 

FOR LOADINGS, ISIDE*2 IF RIGHT MATRIX IS TO BE ROTATED AND LEFT 

MATRIX USED FOR LOADINGS 

F=UNROTATED FACTOR MATRIX 

FSTAR=ROTATEO FACTOR MATRIX 

VSTAR= ROTATED EIGENVECTOR MATRIX 

FF=FR ACT ION ^ OF EIGENVALUE SUM TO BE USED IN FACTOR REPRESENTATION 

DIMENSION CONTRL(4),DATA<NVAR,NLOC),AVG(l),SIGai.Z(N\/AR.NLOCI 

DIMENSION C(1),E(1),V(N,N) ,F ( NSEC , N) , VSTAR (N, N ), FSTAR ( NSEC ,N ) 
DIMENSION TV ( 51 ) *H(1) ,FF(1) ,B(1),D(1) , INDEX ( 1 ) ,R ( 1 ) 

M*N*(N*i)/2 
MODE=CONTRL( 1 ) 

MTYPE=CONTRL( 2) 

I S I DE-CON TRL ( 3 ) 

EF*CONTRL(4) 

CALL COICONTRL.NVAR, NLOC, N, DATA, AVG, SIG,Z,C) 

CALL WRICO(CONTRL, NVAR, NLOC, N, DATA, AVG, SIG, Z,C, INDEX! 


BE ROTATED AND LEFT 


FACTOR REPRESENTATION 


gg 

fig 



0041 


DO 1 1=1, M 

0042 


0(IK(I) 

0043 

1 

CONTINUE 

0044 


CALL EIGEN(D,V»N,0) 

0045 


KK=0 

0046 


DO 2 J=1,N 

0047 


DO 3 1 = 1, N 

0048 


If tl.GT.J) GO TO 3 

0049 


KK=KK+1 

0050 


IF tl.EQ.J) E(I)*D(KK) 

0051 

3 

CONTINUE 

0052 

2 

CONTINUE 

0053 


DO 60 L=1,N 

0054 


IF (E(L).LT.O.O) WRITE 

0055 

61 

FORMAT ( 1X,25HWARNING- 

0056 

1 

0 SO. RT. EXISTS) 

0057 


E(L)=ABS(E(L) ) 

0058 

60 

CONTINUE 

0059 


NEV = 0 

0060 


S=0.0 

0061 


FN=N 

006 2 


EP=EF*FN 

0063 


DO 5 K*1 * N 

0064 


IF (S.GT.EP) GO TO 5 

0065 


S=S*E l K ) 

0066 


NEV=NE V+ 1 

0067 

5 

CONTINUE 

0068 


IF CMTYPE.EQ.l) P*NSEC 

0069 


IF (MTYPE.EQ.2) P=1.0 

0070 


IF IM0DE.EQ.2) GO TO 6 

0071 


DO 7 1=1, NEV 

0072 


DO 8 J=1,NL0C 

0073 


S*0.0 

0074 

• f" 

DO 9 K=i , NVAR 

0075 


S*S+V (K,I)*Z(K,J) 

0076 

9 

CONTINUE 

0077 


F( J, I ) =S/SQRT ( E ( I ) ) 

0078 

8 

CONTINUE 

0079 

7 

CONTINUE 

0080 


GO TO 10 

0081 

6 

no 11 1=1, NVAR 

0082 


DO 12 J=l , NE V 


(108,61) E(L) 

NEG. EIGENVALUE! »F12*5»33H) MADE POSITIVE S 
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0083 


0084 


0085 


0086 

13 

0087 


0088 

12 

0089 

11 

0090 

10 

0091 


0092 


0093 


0094 


0095 


0096 


0097 


0098 


0099 


0100 

17 

0101 


0102 


0103 


0104 

15 

0105 

14 

0106 


0107 


0108 


0109 


0110 


0111 


0112 

21 

0113 

20 

0114 


0115 


0116 


0117 


0118 


0119 


0120 


0121 

24 

0122 


0123 

23 

0124 

22 


EMAX*ABS(Z< I » J) ) 


S*0.0 

DO 13 K*1 , NLOC 
$*S*Z( I,K)*V(K,J) 

CONTINUE 

F C I • J )*S/SQRT IE ( J) > 

CONTINUE 

CONTINUE 

CONTINUE 

CALL WREIGt V»E|F *N f NSEC ? NEV* INDEX ) 

ERMS*0*0 
EHAX=0.0 
DO 14 1*1, NVAR 
DO 15 J*l» NLOC 
$* 0.0 

00 17 K=l» NEV 

IF ( MQDE.EQ. 1 ) S*S-«-Vl I *K) *SQRT( E( KU *F( J«KI 
IF (M0DE.EQ.2) S=S*FIl »K> *SQRT ( E( K ») *V U* K» 
CONTINUE 
Z( It J 1 *Z< I . J>-S 

ERHS = ERMS4ZU.J)**2 

IF (ABS(Z( It JH.GT.EMAX) E HA X* A BS ( Z ( I * J ) ) 
CONTINUE 

CONTINUE . „ 

CALL WRERR(ERHS,EMAX,Z, NVAR, NLOC, INDEX) 

IF (M0DE.EQ.2) GO TO 18 
IF (ISIDE.EQ.2) GO TO 19 
DO 20 J*1»NEV 
DO 21 1*1, NLOC 

FSTAR ( I , J)*F( I » J)*SQRT (E( J) ) 

CONTINUE 

CONTINUE . 

CALL VARMX ( NLOC, NEV, FSTAR, NC, TV, H,FF,B,IER) 

DO 22 1*1, NEV 
DO 23 j*l , NEV 
Mrt*NEV* ( J“1 ) ♦I 

s*o.c 

DO 24 K*1 , NLOC 

S*S*F STARK, J)*F(K,I ) 

CONTINUE 

R(MM)*S/(SQRT(E< 1 1 )*P) 

CONTINUE 

CONTINUE 



to 

in 





0125 


CALL WRTRA(NEV*R* INDEX) 

0126 


00 25 I *1, NVAR 

0127 


DO 26 J*1 » NEV 

0128 


$*o«o 

0129 


DO 27 K*1 » NEV 

0130 


MM«NEV*t J-1>*K 

0131 


S*S*VU,K)*R(MM) 

0132 

27 

CONTINUE 

0133 


VSTAR( I , J ) *S 

0134 

26 

CONTINUE 

0135 

25 

CONTINUE 

0136 


GO TO 28 

0137 

19 

CONTINUE 

0138 


DO 29 J*l,NEV 

0139 


DO 30 I *1* NVAR 

0140 


VSTARI I,J)*V(I, J)*SQRT(EC J)l 

0141 

30 

CONTINUE 

0142 

29 

CONTINUE 

0143 


CALL VARMXINVAR, NEV, VSTAR ,NC ,TV ,H, FF , B, IER) 

0144 


DO 31 1*1, NEV 

0145 


DO 32 J*l , NE V 

0146 


MM*NEV* I J- 1 ) ♦ I 

0147 


$* 0.0 

0148 


00 33 K*l , NVAR 

0149 


S*S+V(K,I)*VSTAR(K,J) 

0150 

33 

CONTINUE 

0151 


R|MM)*$/SQRTIEII > 1 

0152 

32 

CONTINUE 

0153 

31 

CONTINUE 

0154 


CALL WRTRAINEV.R, INDEX) 

0155 


DO 34 I *1, NLOC 

0156 


DO 35 J*1 » NE V 

0157 


S*0,0 

0158 


DO 36 K*i » NEV 

0159 


MM*NEV*( J-D+K 

0160 


S*S*FI I ,K)*R(MM) 

0161 

36 

CONTINUE 

0162 


ESTAR I I » J ) = S 

0163 

35 

CONTINUE 

0164 

34 

CONTINUE 

0165 


GO TO 28 

0166 

18 

CONTINUE 
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167 


IF I ISI0E.EQ.21 GO TO 37 

166 


00 38 J=1*NE V 

169 


DO 39 1*1 * NLOC 

170 


V START I , J )=V( I , J ) *SQRT(E ( J ) ) 

171 

39 

CONTINUE 

172 

38 

CONTINUE 

173 


CALL VARMX ( NLOC » NEV* VSTAR , NC»TV ,H, FF, B, IER > 

176 


DO 60 1*1, NEV 

175 


DO 61 J*1 , NEV 

176 


MM*NEV*(J~1)+I 

177 


S*0.0 

170 


DO 62 K*l, NLOC 

179 


S*S+V(K, I )*VSTAR(K, J) 

180 

62 

CONTINUE 

181 


R{KM)*S/SQRT(E( I ) ) 

182 

61 

CONTINUE 

183 

60 

CONTINUE 

186 


CALL WRTRAINEVtR, INDEX) 

185 


DO 63 1=1* NVAR 

186 


DO 66 J=1 , NEV 

187 


S*0.0 

188 


DO 65 K= 1 , NEV 

189 


MM*NEV* (J-l ) +K 

190 


S*S^F(I,K)*R(MM) 

191 

65 

CONTINUE 

192 


FSTARC I * J)*S 

193 

66 

CONTINUE 

196 

63 

CONTINUE 

195 


GO TO 28 

196 

37 

CONTINUE 

197 


DO 66 J = l » NE V 

198 


DO 67 1*1, NVAR 

199 


FSTARl I,J)*F( I,J)*SQRT(E( J) ) 

200 

67 

CONTINUE 

201 

66 

CONTINUE 

202 


CALL VARMX ( NVAR , NEV, FSTAR.NC* TV *H,FF,B* IER) 

203 


DO 68 1*1, NEV 

206 


00 69 J*1 , NEV 

205 


MN*NEV* ( J- 1 ) ♦ I 

206 


S*0.0 

207 


DO 50 K*1 , NVAR 

208 


S*$+F(K, I)*FSTAR(X,J) 


NJ 

-J 


0209 

50 

CONTINUE 

0210 


R(MM)*S/( SORT CEI I )I*P) 

0211 

49 

CONTINUE 

0212 

48 

CONTINUE 

0213 


CALL WRTRA(NEV t R, INDEX) 

0214 


DO 51 1*1 * NLOC 

0215 


DO 52 J*1 » NE V 

0216 


S*0.0 

021? 


DO 53 K*1 1 NE V 

0218 


MM*NEV* ( J-l ) ♦K 

0219 


$*S+V( IrK»*R(MH) 

0220 

53 

CONTINUE 

0221 


VSTAR ( 1 1 J ) *S 

0222 

52 

CONTINUE 

0223 

51 

CONTINUE 

0224 

28 

CONTINUE 

0225 


CALL HRSTR (VSTAR tE ST AR.Nf NSEC tNEVf INDEX) 

0226 


RETURN 

0227 


END 
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0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
0011 
0Q12 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 
0021 
0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 


101 

100 

102 

103 


SUBROUTINE MR ICO ( CONTRL » N VAR, NLOC « N, DATA. AVG,SIG,Z,C, INDEX) 

FORMAT U8H OUTPUT OF SUBR CO) 

FORMAT (8H MODE = f I2»10X,9H MTYPE = »I2*8H NVAR = »I5»8H NLOC = » 
15, 9H ISIOE - , 1 5, 6H EF * ,F10.5) 

FORMAT ( IX , 12H DATA MATRIX) 

FORMAT (IX) 

FORMAT (IX, 19H NORMED DATA MATRIX) 

FORMAT (6X,20H VECTOR LENGTH) 

FORMAT (IX, I5.F20.5) ^ , 

FORMAT (6X,42H MEAN STANOARO OEVIATION) 

FORMAT (IX, I 5, 2F20. 5 ) 

FORMAT ( IX, 19H CORRELATION MATRIX) 

FORMAT (IX, 1 AH COSINE MATRIX) 

FORMAT (IH1) „ , , 

DIMENSION DATA (NVAR, NLOC) , Z(NVAR, NLOC ),C(l),SIG(l), AVG( 1) , INDEX( I ) 

DIMENSION CONTRL ( l ) 

MODE-CONTRL < l ) 

MTYPE*CONTRL ( 2 ) 

ISIDE=C0NTRL(3) 

F:F=CONTRL ( 4) 

WRITE (108,12) 

WRITE (108,1) 

WRITE (108,2) MODE, MTYPE, NVAR.NLOC, I SIDE,EF 
WRITE (108,3) 

CALL WRMAT(DATA, NVAR, NLOC, INDEX) 

WRITE (108,4) 

IF (M0DE.EQ.2) GO TO 100 
WRITE (108,8) 

00 101 I- 1 , NVAR 

WRITE (108,9) I * AVG ( I ) »SI G ( I ) 

CONTINUE 
GO TO 103 
WRITE (108,6) 

DO 102 1*1, NLOC 
WRITE (108,7) I , SIG( I ) 

CONTINUE 
CONTINUE 
WRITE (108,4) 

WRITE <108.5) 
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0040 

CALL WRMAT ( Z»NVAR , NLOC 

, INDEX) 

0041 

WRITE (108,4) 


0042 

IF (MTYPE.EQ.l) WRITE 

( 108,10) 

0043 

IF (MTYPE.EQ.2) WRITE 

(108,11) 

0044 

CALL WRTRKC.N, INDEX) 


0045 

RETURN 


0046 

END 



0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
0011 
0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 
0021 
0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 


% 


1 

2 

3 

4 

5 

6 

7 

8 
9 

1C 


50 


FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 


EV) 


SUBROUTINE WRE1G( V*E*F ,N, NSEC, NEV, INDEX) 

FORMAT MHD ^ _ 11U Cltu 

{ IX 1 30H INDEX EIGENVALUES CUM. SUM 

(1X,I6,2F12.5) 

( IX ) 

(1X,23H EIGENVECTORS (COLUMNS)) 

(1X,18H FACTORS (COLUMNS)) 
(1X*7HE-VALUE,5(I3,F10.5)) 

(IX, 3110) 

( IX, 8HE-VECT0R *5 ( I 3,F10.5) ) 

( 1X,6HFACT0R,5( I3,F10.5) ) 

DIMENSION V(N,N),E(1)*F(NSEC*N) ,INDEX(1) 

WRITE (108,1) 

WRITE (108,2) 

S=0.0 

DO 50 1=1, N 
S*S+E( I) 

WRITE (108,3) I * E ( I ) * S 

CONTINUE 

WRITE (108,4) 

WRITE (108,5) 

CALL WRMAT(V,N,NEV, INDEX) 

WRITE (108,1) 

WRITE (108,6) 

CALL WRMAT(F,NSEC,NEV, INDEX) 

WRITE (106,8) N, NSEC *NEV 

WRITE (106,7) ((I, E(I) ),I*l»NEV) 

WRITE (106,9) ( ( (J,Vd , J ) ) ,!=1»N) , J-UNEV) 

WRITE (106,10) ( ( ( J,F( I , J ) ) » 1=1 »NSEC ) , J-l »NEV) 

RETURN 

END 




0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
0011 
0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 


SUBROUTINE WRSTRI VSTAR, F STAR, N, NSEC* NEV* INDEX I 
C VSTAR»N X NEV MATRIX OF VARI MAX ROTATED EIGENVECTORS 

C FSTAR=NSEC X NEV MATRIX OF VARI MAX ADJUSTED FACTORS 

C N=NUMBER OF COMPONENTS IN EIGENVECTORS 

C NSEC=OTHER DIMENSION FROM N (THAT IS, NSEC=NVAR IF N=NLOC ANO 

C NSEC=NLOC IF N=NVAR) 

C NEV=»NUMBER OF EIGENVALUES AND EIGENVECTORS USED 

C INDEX=WORK VECTOR OF MAX COMPONENTS, MAX=PAX ( NVAR » NLOC ) 

1 FORMAT ( 1 H 1 > 

3 FORMAT ( IX , 1 3H VSTAR MATRIX) 

A FORMAT ( IX , 13H FSTAR MATRIX) 

DIMENSION VSTAR (N,N) ,F STAR ( NSEC ,N ) ,INDEX(1) 

WRITE (108,1) 

WRITE (108,3) 

CALL WRMATl VSTAR ,N, NEV, INDEX) 

WRITE (108,1) 

WRITE (108,4) 

CALL WRMAT(FSTAR, NSEC, NEV, INDEX) 

RETURN 

END 
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0001 


SUBROUT INE WRERR<ERMS*EMA X,Z*NVAR, NLOC « I NOEX I 

0002 

1 

FORMAT (IHl) 

0003 

2 

FORMAT ( 1 X 1 6 H ERMS* » F20* 5 » 6H EMAX=,F20*5) 

0004 

3 

FORMAT (IX) 

0005 

4 

FORMAT (1X.24H ERROR OF REPRESENTATION) 

0006 


DIMENSION Z ( NVAR , NLOC ) , I NDEX ( 1 ) 

0007 


WRITE (108,1) 

0008 


WRITE (108,2) ERMS»EMAX 

0009 


WRITE (108,3) 

0010 


WRITE (108,4) 

0011 


CALL WRMAT(Z,NVAR, NLOC, INDEX) 

0012 


RETURN 

0013 


END 
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0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
0011 
0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 
0021 
0022 

0023 

0024 

0025 

0026 


i 

SUBROUTINE WRMATI A,N*M* INDEX) 

PRINT A NX* MATRIX WITH COLUMN AND ROW NUMBERING. A*MATRIX 

INDEX* M-OI MENS IONAL WORK VECTOR j 

DIMENSION A(N,M> , INDEX(M) 

1 FORMAT ( IX, 1 13,9112) 

2 FORMAT ( IX , 1 5 , 10F12. 5 ) 

9 FORMAT (IX) 

MR*M+10 
ICB*-9 
DO 3 K*1 , M 
INDEX ( K ) *K 

3 CONTINUE 

4 MR*MR- 10 
ICB*ICB*10 

IF I MR-10 ) 5,5,6 

5 ICE*ICB+MR-1 
GO TO 7 

6 lCE*ICB+9 

7 WRITE (108,9) 

WRITE (108,1) ( INDEX(K) fK*ICB, ICE) 

DO 8 I R*1 , N 

WRITE (108,2) IR,(A(IR,IC),IC*ICB,ICE) 

8 CONTINUE 

IF (MR.GT.10) GO TO 4 

RETURN 

END 
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0001 

0002 C 

0003 C 

000 A c 

0005 C 

0006 

0007 

0008 

0009 

0010 
0011 
0012 
0013 
0016 

0015 

0016 

0017 

0018 

0019 

0020 
0021 
0022 
0023 
0026 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 
0036 


SUBROUTINE WRTRIC A *N, INDEX) ci cmpnt 

A*TRIANGULAR MATRIX WRITTEN IN VECTOR SEQUENCE, fj 1 ***^} 

A( 21*12*1) ELEMENT, A(3)=I2»2) ELEMENT, A (61* (3, l I ELEMENT, E C. 

N*ORDER OF MATRIX 
M=N*( N+1 1/2 

FORMAT (IX, 113,9112) 

FORMAT ( IX, 15, 10F12. 5) 

FORMAT (IX) 

DIMENSION A( 1) , INDEX! 1 ) 

NC*0 
MR*N^10 



ICB*-9 

00 l K= 1 , N 

INDEX ( K 1 

CONTINUE 

CONTINUE 

NC*NC*1 

MRsMR-10 

ICBsICB+10 

IF (MR-10) 2,2*3 

IGE*ICB*MR-1 


GO TO 6 
lCE*ICB+9 
WRITE (108,22) 

WRITE (108,20) ( INDEX(K) ,K*ICB, ICE) 



L*0 

DO 5 IR*ICB,N 
L*L*l 

KRB*i + ( I R- 1 ) * I R / 2 ♦ ( NC-1 » * TO 

IF (L.LE.10) KRE*KRB+L-l 

IF (L.GT.10) KRE=KRB*9 

WRITE (108,21) IR, ( A(K) ,K=KRB,KRE) 

CONTINUE 

IF (MR.GT.10) GO TO 6 


0035 RETURN 

0036 END 
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0001 SUBROUTINE WRTRAINE V »R * INDEX ) 

0002 DIMENSION RINEV.NEV) , INDEXtl) 

0003 1 FORMAT (1H1.23H VARIMAX TRANSFORMATION) 

0004 WRITE ( 108 1 1 ) 

0005 CALL WRMAT<R,NEV,NEV, INDEX) 

0006 RETURN 

0007 END 




0002 C 

0003 C 

0004 C 

0005 C 

0006 C 

0007 C 
6008 C 

0009 C 

0010 C 

0011 
0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 
0021 
0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 


SUBROUTINE CO(CONTRL *NVAR » NIOC . N, DATA ,AVG*$IG*Z*C1 
NsNVAR *F MODE* l* N=NLOC IF M0DE=2 _ ^ 

M=N*<N*h/2 = NUMBER OF ELEMENTS IN UPPER TRIANGULAR PART OF R 

MODE* 1 FOR R-MODE ANALYSIS, M0DE=2 F0R 1 onus ARE 
DATA*NVAR X NLOC DATA MATRIX, COLUMNS ARE LOCATIONS, ROWS ARE 

VARIABLES, 

R*NXN CORRELATION MATRIX IF MTYPE*1 

^IS N ST0RED E AS A A R VECT0R, T THE UPPER TRIANGULAR PART IS STORED 
CGLUMN^^ iqjy CONTRL ( 4 } , DATA! N VAR, NLOC ) , AVG( 1 ) ,SIG( 1 ) , ZI NVAR, 


DIMENSION C ( 1 ) 
MODE=CONTRL( l) 
MTYPE*CONTRL (2) 

IF (M0DE.EQ,2) GO TO 1 
DO 2 1*1, NVAR 
AVG( I ) =0.0 

2 CONTINUE 

IF (MTYPE.EQ.2) GO TO 14 

DO 3 J=1,NL0C 

DO 4 1*1, NVAR 

AVGl I ) = A VG ( I )+DATA ( I , J) 

4 CONTINUE 

3 CONTINUE 
FNLOC*NLOC 

00 5 1*1, NVAR 
AVG< II =AVG (I) /FNLOC 

5 CONTINUE 

DO 6 J = l , NLOC 
DO 7 1*1, NVAR 
Z( I, J)*DATA( I, J )-AVG( I ) 

7 CONTINUE 

6 CONTINUE 
GO TO 17 

l CONTINUE 

DO 8 J*1 , NLOC 
AVGl J) =0,0 

8 CONTINUE 

IF (MTYPE.EQ.2) GO TO 14 
DO 9 1=1, NVAR 
00 10 J = l» NLOC 


IN 

NLOC ) 
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0042 

-■ ■ .. V •" 

AVG IJ ) *AVG1 J ! ♦DATA! I ,J 1 

0043 

10 

CONTINUE 

0044 

9 

CONTINUE 

0045 


FNVAR=NVAR 

0046 


DO 11 J* l » NLOC 

0047 


AVG< J»*AVG( JJ/FNVAR 

0048 

11 

CONTINUE 

0049 


DO 12 1*1 » NVAR 

0050 


DO 13 J* l * NLOC 

0051 


Z1I#J)*DATA( I v JI-AVGUi 

0052 

13 

CONTINUE 

0053 

12 

CONTINUE 

0054 


GO TO 17 

0055 

14 

DO 15 1*1# NVAR 

0056 


DO 16 J*lt NLOC 

0057 


Z( I, J)*DATA{ I, J) 

0058 

16 

CONTINUE 

0059 

15 

CONTINUE 

0060 

17 

IF (MODE. EQ. 21 GO TO 18 

0061 


DO 19 I *1 » NVAR 

0062 


S IGC I ) *0.0 

0063 

19 

CONTINUE 

0064 


DO 20 J*1 * NLOC 

0065 


DO 21 1*1# NVAR 

0066 


SI.G( 1 l*SIG( I ) +2( I * J)**2 

0067 

21 

CONTINUE 

0068 

20 

CONTINUE 

0069 


DO 22 1*1# NVAR 

0070 


IF (MTYPE.EQ.l) SIG( I)*SIG( II/FNLOC 

0071 


SIGI I )=SQRT( SIG( I ) ) 

0072 

22 

CONTINUE 

0073 


00 23 J*i» NLQC 

0074 

•* 

DO 24 I *1 # NVAR 

0075 


Z(I#J)*Z(I«J)/SIG(I) 

0076 

24 

CONTINUE 

0077 

23 

CONTINUE 

0078 


GO TO 25 

0079 

18 

DO 26 J*1 ♦ NLOC 

0080 


S I G ( J ) *0.0 
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0081 

26 

CONTINUE 

0082 


00 27 I*1,NVAR 

0083 


DO 28 J*1 * NLOC 

0084 


SIG< JI=SIG( J)*Z(It J)**2 

0085 

28 

CONTINUE 

0086 

27 

CONTINUE 

0087 


FNVAR*NVAR 

0088 


00 29 J*1 , NLGC 

0089 


IF (MTYPE.EQ.l) SIGt JI*SIG( JI/FNVAR 

0090 


SIG( J I *SQRT( SIG< J I I 

0091 

29 

CONTINUE 

0092 


DO 30 I = 1 » N V AR 

0093 


DO 31 J* 1 i NLOC 

0094 


Z 1 1 • J I =Z ( I»J)/SIG»JI 

0095 

31 

CONTINUE 

0096 

30 

CONTINUE 

0097 

25 

CONTINUE 

0098 


IF (M0DE.EQ.2) GO TO 32 

0099 


KK*0 

0100 


DO 33 1=1, N 

0101 


DO 34 J*l,I 

0102 


KK*KK*l 

0103 


S*0.0 

0104 


DO 35 JJ=1,NL0C 

0105 


S*S*Z(I f JJ)*Z( J, JJl 

0106 

35 

CONTINUE 

0107 


IF ( MTYPE. EQ. 1 '# S-S/FNLOC 

0108 


C ( KK I *S 

01Q9 

34 

CONTINUE 

0110 

33 

CONTINUE 

0111 


GO TO 39 

0112 

32 

CONTINUE 

0113 


KK*0 

0114 


DO 36 1*1, N 

0115 


DO 37 J=i,I 

0116 


KK*KK+ 1 

0117 


S*0.0 

0118 


DO 38 I I = 1,NVAR 

0119 


S*S+Z(II,I)+Z(II,4! 

0120 

36 

CONTINUE 

0121 


IF (MTYPE.EQ.il S=S/FNVAR 

0122 


C IKK 1 * S 


0123 

37 

CONTINUE 

0124 

36 

CONTINUE 

0125 

39 

CONTINUE 

0126 


RETURN 

0127 


END 


000 1 c 

0002 C 

0003 C 

0004 C 

0005 C 

0006 C 

0007 C 

0008 C 

0009 C 

0010 C 

0011 C 

0012 C 

0013 C 

0014 C 

0015 C 

0016 C 

0017 C 

0018 C 

0019 C 

0020 C 

0021 C 

0022 C 

0023 C 

0024 C 

0025 C 

0026 C 

0027 C 

0028 C 

0029 C 

0030 C 

0031 C 

0032 C 

0033 C 

0034 C 

0035 C 

0036 C 

0037 C 

0038 C 

0039 C 

0040 C 

0041 


SUBROUTINE EIGEN 
PURPOSE 

COMPUTE EIGENVALUES AND EIGENVECTORS OF A REAL SYMMETRIC 
MATRIX 

USAGE 

CALL EIGEN(A,R,N,MV) 

DESCRIPTION OF PARAMETERS 

A - ORIGINAL MATRIX ( SYMMETRIC) * DESTROYED IN COMPUTATION. 
RESULTANT EIGENVALUES ARE DEVELOPED IN DIAGONAL OF 
MATRIX A IN DESCENDING ORDER. 

R - RESULTANT MATRIX OF EIGENVECTORS (STORED COLUMNWISE, 

IN SAME SEQUENCE AS EIGENVALUES) 

N - ORDER OF MATRICES A AND R 
MV- INPUT CODE 

0 COMPUTE EIGENVALUES AND EIGENVECTORS 

1 COMPUTE EIGENVALUES ONLY (R NEED NOT BE 
DIMENSIONED BUT MUST STILL APPEAR IN CALLING 
SEQUENCE) 

REMARKS 

ORIGINAL MATRIX A MUST BE REAL SYMMETRIC (STORAGE MODE*!) 
MATRIX A CANNOT BE IN THE SAME LOCATION AS MATRIX R 

SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED 
NONE 

METHOD 

DIAGONAL I ZATION METHOD ORIGIN*. ED BY JACOBI AND ADAPTEC 
BY VON NEUMANN FOR LARGE COMPUTERS AS FOUND IN •MATHEMATICAL 
METHODS FOR DIGITAL COMPUTERS*, EDITED BY A. RALSTON AND 
H.S. WILF, JOHN WILEY AND SONS, NEW YORK, 1962, CHAPTER 7 


SUBROUTINE El GEN (A ,R,N»MV ) 


EIGE 

10 

EIGE 

20 

EIGE 

30 

EIGE 

40 

EIGE 

50 

EIGE 

60 

EIGE 

70 

EIGE 

80 

EIGE 

90 

EIGE 

100 

EIGE 

110 

EIGE 

120 

EIGE 

130 

EIGE 

140 

EIGE 

150 

EIGE 

160 

EIGE 

170 

EIGE 

180 

EIGE 

190 

EIGE 

200 

EIGE 

2iC 

EIGE 

220 

EIGE 

230 

EIGE 

240 

EIGE 

250 

EIGE 

260 

EIGE 

270 

EIGE 

280 

EIGE 

290 

EIGE 

300 

EIGE 

310 

EIGE 

320 

EIGE 

330 

EIGE 

340 

EIGE 

350 

EIGE 

360 

EIGE 

370 

EIGE 

380 

EIGE 

390 

EIGE 

400 

EIGE 

410 
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0042 

0043 

0044 

C 

r 


DIMENSION Ml)tR(l) 

EIGE 

EIGE 

EIGE 

420 

430 

440 

0045 

c 



EIGE 

450 

0046 

c 


IF A DOUBLE PRECISION VERSION OF THIS ROUTINE IS DESIRED, THE 

EIGE 

460 

0047 

c 


C IN COLUMN 1 SHOULD BE REMOVED FROM THE DOUBLE PRECISION 

EIGE 

470 

0048 

c 


STATEMENT WHICH FOLLOWS. 

EIGE 

480 

0049 

c 



EIGE 

490 

0050 

c 


OOUBLE PRECISION A, R , ANORM , ANRMX, THR, X , Y, SINX , S I NX2, COSX , 

EIGE 

500 

0051 

c 

1 

C0SX2,SINCS, RANGE 

EIGE 

510 

0052 

c 



EIGE 

520 

0053 

c 


THE C MUST ALSO BE REMOVED FROM DOUBLE PRECISION STATEMENTS 

EIGE 

530 

0054 

c 


APPEARING IN OTHER ROUTINES USED IN CONJUNCTION WITH THIS 

EIGE 

540 

0055 

c 


ROUTINE. 

EIGE 

550 

0056 

c 



EIGE 

560 

0057 

c 


THE DOUBLE PRECISION VERSION OF THIS SUBROUTINE MUST ALSO 

EIGE 

570 

0058 

c 


CONTAIN DOUBLE PRECISION FORTRAN FUNCTIONS. SORT IN STATEMENTS 

EIGE 

580 

0059 

c 


40, 68, 75, AND 78 MUST BE CHANGED TO DSQRT. ABS IN STATEMENT 

EIGE 

590 

0060 

c 


62 MUST BE CHANGED TO DABS. THE CONSTANT IN STATEMENT 5 SHOULD 

EIGE 

600 

0061 

c 


BE CHANGED TO l.OD-12. 

EIGE 

610 

0062 

c 



EIGE 

620 

0063 

r 



EIGE 

630 

0064 

C 


• 

EIGE 

640 

0065 

C 


GENERATE IDENTITY MATRIX 

EIGE 

650 

0066 

C 



EIGE 

660 

0067 


5 

RANGE=1.0E-6 

EIGE 

670 

0068 



IFIMV-l) 10,25,10 

EIGE 

680 

0069 


10 

IQs-N 

EIGE 

690 

0070 



DO 20 j*l,N 

EIGE 

700 

0071 



IQ=IQ*N 

EIGE 

710 

0072 



DO 20 I *1 » N 

EIGE 

720 

0073 



I J*IQ+I 

EIGE 

730 

0074 



R( I J 1=0.0 

EIGE 

740 

0075 



IFU-J) 20,15,20 

EIGE 

750 

0076 


15 

R< I J > = 1.0 

EIGE 

760 

0077 


20 

CONTINUE 

EIGE 

770 

0078 

c 



EIGE 

780 

0079 

c 


COMPUTE INITIAL AND FINAL NORMS (ANORM AND ANORMX ) 

EIGE 

790 

0080 

c 



EIGE 

800 

0081 


25 

ANORM=0.0 

EIGE 

810 

0082 



DO 35 1*1, N 

EIGE 

820 

0083 



DO 35 J*I,N ' 

EIGE 

830 
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00 8 A 



IF € I-JJ 30.35,30 

0085 


30 

IA*i+( J*j-J)/2 

0086 



ANORM= ANORM+A ( I A ) * A ( I A » 

0087 


35 

CONTINUE 

0088 



IF(ANORM) 165, 165, AO 

0089 


AO 

ANORM^l • A1 A*SQRT I ANORM ) 

0090 



ANRMX=ANORM*RANGE/FLOAT(N) 

0091 

C 



0092 

c 


INITIALIZE INDICATORS ANO COMPUTE THRESHOLD, 

0093 

c 



009A 



IND*0 

0095 



THR*ANORM 

0096 


A5 

THR“THR/FLOAT( N) 

0097 


50 

L*1 

0098 


55 

M*L + l 

0099 

c 



0100 

c 


COMPUTE SIN AND COS 

0101 

c 



0102 


60 

MQ*{M*M-M)/2 

0103 



LQ*IL*L-L)/2 

010A 



LM=L*MQ 

0105 


62 

I F ( ABSIA(LM) )-THR> 130,65,65 

0106 


65 

IND=1 

0107 



LL-L^LQ 

0108 



MM*M*MQ 

0109 



X=0» 5* (A(LL)-A(MM) ) 

0110 


68 

Y=-A ( LM ) / SORT ( A ( LM) *A ( LM ) +X*X ) 

0111 



I F I X ) 70,75,75 

0112 


70 

Y=-Y 

0113 


75 

S I NX=Y / SORT ( 2«0*( 1«0>( SQRT(1.0-Y*Y ))) ) 

011A 



SINX2=SINX*S I NX 

0115 


78 

COSX* SQRTI1.0-SINX2) 

0116 



C0SX2=C0SX*C0SX 

0117 



SINCS =SINX*COSX 

0118 

c 



0119 

c 


ROTATE L AND M COLUMNS 

0120 

c 



0121 



ILQ*N* ( L-l ) 

0122 



IMQ-N* ( M- 1 ) 

0123 



DO 125 1*1* N 

012A 



IQ* ( I*I-I>/2 

0125 



IFTl-LI 80,115,80 


E1GE 8 AO 
EIGE 850 
EIGE 860 
EIGE 870 
EIGE 880 
EIGE 890 
EIGE 900 
EIGE 910 
EIGE 920 
EIGE 930 
EIGE 9A0 
EIGE 950 
EIGE 960 
EIGE 970 
EIGE 980 
EIGE 990 
EIGE1000 
EIGE 10 10 
EIGE1C20 
EIGE1030 
EIGELOAO 
EIGE1C50 
* E IGE1060 
EIGE 1070 
EIGE1080 
EIGE1090 
EIGE 1100 
EIGE 11 10 
EIGE 1120 
EIGE1130 
EIGE11A0 
EIGE1 150 
EIGE 1160 
EIGE1170 
EIGE 1180 
E1GE1190 
EIGE1200 
El GE 1 210 
EIGE 1220 
EIGE1230 
EIGE12A0 
EIGE1250 
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0126 

80 

IFtl-M) 85, 115 » 90 

E1GE1260 

0127 

85 

IM*I*MQ 

EIGE1270 

0128 


GO TO 95 

EIGE1280 

0129 

90 

IM=M+IQ 

El GE 1290 

0130 

95 

IF ( I-L ) 100,105,105 

EIGE1300 

0131 

100 

IL-I+LQ 

EIGE1310 

0132 


GO TO no 

EIGE1320 

0133 

105 

IL=L+I Q 

EIGE1330 

0134 

110 

X=*A( IL )*COSX-A( IM)*SINX 

EIGE1340 

0135 


A ( IM )=A ( I L ) *S I NX* A ( I M) *COSX 

EIGE1350 

0136 


AtIU-X 

EIGE1360 

0137 

115 

IF(MV-l) 120,125,120 

EIGE1370 

0138 

120 

ILR*ILQ+I 

EIGE1380 

0139 


IMR= I MQ+ I 

EIGE1390 

0140 


X=Rt ILR)*COSX-RtIMR)*SINX 

E1GE14C0 

0141 


R t IMR ) =R 1 1 LR ) *SI NX+R 1 1 MR ) *COSX 

EIGE1410 

0142 


R ( I LR ) *X 

EIGE1420 

0143 

125 

CONTINUE 

EIGE1430 

0144 


X*2*0*A ( LM ) *SI NCS 

E IGE1440 

0145 


Y=At LL >*C0SX2*A<MM)*SINX2-X 

EIGE1450 

0146 


X-AtLLl*SINX2 + A(MM)*COSX2«-X 

EIGE1460 

0147 


At LM) * t A t LL ) -A t MM) ) *SINCS+A(LM) *( C0SX2-SINX2) 

EIGE1470 

0148 


At LL ) *Y 

EIGE1480 

0149 


At MM)*X 

EIGE1490 

0150 

C 


EIGE1500 

0151 

C 

TESTS FOR COMPLETION 

EIGE1510 

0152 

C 


EIGE1520 

0153 

C 

TEST FOR M * LAST COLUMN 

EIGE1530 

0154 

c 


EIGE1540 

0155 

130 

IF t M-N ) 135,140,135 

E IGE1550 

0156 

135 

M*M*l 

EIGE1560 

0157 


GO TO 60 

E IGE1570 

0158 

C 


E IGE1580 

0159 

c 

TEST FOR L * SECOND FROM LAST COLUMN 

E IGE1590 

0160 

c 


EIGE1600 

0161 

140 

iFCL-tN-m 145,150,145 

E IGE1610 

0162 

145 

L*L*i 

E I GE 1620 

0163 


GO TO 55 

EIGE1630 

0164 

150 

IFtIND-l) 160,155,160 

EIGE1640 

0165 

155 

IND»0 

EIGE1650 

0166 


GO TO 50 

E I GE 1660 

0167 

c 

- ■ ■ 'V 

EIGE1670 
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0168 

0169 

0170 

0171 

0172 

0173 

0174 

0175 

0176 

0177 

0178 

0179 

0180 
0181 
0182 

0183 

0184 

0185 

0186 

0187 

0188 

0189 

0190 

0191 

0192 

0193 
019* 
0195 


C COMPARE THRESHOLD WITH FINAL NORM 

C 

160 IF ( THR-ANRMX I 165.165.45 

C 

C SORT EIGENVALUES AND EIGENVECTORS 

C 

165 IQ*-N 

DO 185 I-l.N 
IQ* IQ+N 

LL-I+I I4I-II/2 
JQ*N*( 1-2) 

DO 185 J=I,N 
JQ*JQ«-N 

MM*J*< J*J-J>/2 
IF C A ( LL )-A { HP ) ) 170,185.185 
170 X*A(LU 

AI LL )*A (MM) 

A ( MM ) *X 

IF (MV- l ) 175,185,175 
175 DO 180 K*1 , N 

‘ ILR*IQ+K 
IMR* JQ-t-K 
X*RC ILR ) 

R ( ILR )*R( I MR ) 

180 R(IMR)*X 

185 CONTINUE 

RETURN 
END 


E IGE l 663 

EIGE1690 

E IGE1700 

EIGE1710 

EIGE1720 

E1GE1730 

E IGE 1740 

EIGE1750 

EIGE1760 

EIGE1770 

EIGE1780 

EIGE1790 

EIGE1800 

E IGE 1810 

EIGE1820 

El GE 1830 

EIGE1840 

EIGE1850 

E IGE 1860 

EIGE1870 

EIGE1880 

EIGE1890 

EIGE1900 

EIGE1910 

E IGE 1920 

E IGE 1930 

E IGE 1940 

E IGE 19 50 






uasas*.**cr**r»* ' h'i, ji ^ 
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N) 

O' 

O' 


0001 C VARM 10 

0002 C •••••••«• VARM 20 

0003 C VARM 30 

0004 C SUBROUTINE VARMX VARM AO 

0005 C VARM 50 

0006 C PURPOSE VARM 60 

0007 C PERFORM ORTHOGONAL ROTATIONS OF A FACTOR MATRIX. THIS VARM 70 

0008 C SUBROUTINE NORMALLY OCCURS IN A SEQUENCE OF CALLS TO SUB- VARM 83 

0009 C ROUTINES CQRRE, EIGEN, TRACE, LOAO, VARMX IN THE PERFORMANCE VARM 90 


0010 

0011 

C 

c 

OF A 

FACTOR ANALYSIS. 

0012 

c 

USAGE 


0013 

0014 

c 

c 

CALL 

VARMX (M»K,A,NC»TV*H»F,D*IERI 

0015 

c 

DESCRIPTION OF PARAMETERS 

0016 

c 

M 

- NUMBER OF VARIABLES AND NUMBER OF ROWS OF MATRIX A. 

0017 

c 

K 

- NUMBER OF FACTORS. 

0018 

c 

A 

- INPUT IS THE ORIGINAL FACTOR MATRIX, AND OUTPUT IS 

0019 

c 


THE ROTATED FACTOR MATRIX. THE ORDER OF MATRIX A 

0020 

c 


IS M X K. 

0021 

c 

NC 

- OUTPUT VARIABLE CONTAINING THE NUMBER OF ITERATION 

0022 

c 


CYCLES PERFORMED. 

0023 

c 

TV 

- OUTPUT VECTOR CONTAINING THE VARIANCE OF THE FACTOR 

0024 

c 


MATRIX FOR EACH ITERATION CYCLE. THE VARIANCE PRIOR 

0025 

c 


TO THE FIRST ITERATION CYCLE IS ALSO CALCULATED. 

0026 

c 


THIS MEANS THAT NC+1 VARIANCES ARE STORED IN VECTOR 

0027 

c 


TV. MAXIMUM NUMBER OF ITERATION CYCLES ALLOWED IN 

0028 

c 


THIS SUBROUTINE IS 50. THEREFORE, THE LENGTH OF 

0029 

c 


VECTOR TV IS 51. 

0030 

c 

H 

- OUTPUT VECTOR OF LENGTH M CONTAINING THE ORIGINAL 

0031 

c 


COMMUNAL I TI ES. 

0032 

c 

F 

- OUTPUT VECTOR OF LENGTH M CONTAINING THE FINAL 

0033 

c 


COMMONALITIES. 

0034 

c 

D 

- OUTPUT VECTOR OF LENGTH M CONTAINING THE DIFFERENCES 

0035 

c 


BETWEEN THE ORIGINAL AND FINAL COMMUNAL IT IES . 

0036 

c 

IER 

- ERROR INDICATOR 

0037 

c 


IER=0 - NO ERROR 

0038 

c 


IER-1 - CONVERGENCE WAS NOT ACHIEVED IN 50 CYCLES 

0039 

c 


OF ROTATION 


VARM 100 
VARM 110 
VARM 120 
VARM 130 
VARM 140 
VARM 150 
VARM 160 
VARM 170 
VARM 180 
VARM 190 
VARM 200 
VARM 210 
VARM 220 
VARM 230 
VARM 240 
VARM 250 
VARM 260 
VARM 270 
VARM 280 
VARM 290 
VARM 300 
VARM 310 
VARM 320 
VARM 330 
VARM 340 
VARM 350 
VARM 351 
VARM 352 
VARM 353 
VARM 354 


0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 

0058 

0059 

0060 
0061 
0062 

0063 

0064 

0065 

0066 

0067 

0068 

0069 

0070 

0071 

0072 

0073 

0074 

0075 

0076 

0077 

0078 

0079 

0080 
0081 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


ROTATION. 

SUBROUTINES AND FUNCTION SUBPROGRAMS REOUIRED 
NONE 

MET KMSER.S VARIMAX 

educational and psychological measurement, VOL XIX, NO. 3, 
1959. 

SUBROUTINE VARMX < M,K, A, NC ,TV , H ,F , D, I ER > 

DIMENSION All), TV ID, H(i), F(l), 0(1) 

ROUBLE PRECISION * S^COMT^sIhStS’sIMT 

THP f MUST ALSO BE REMOVED FROM DOUBLE PRECISION STATEMENTS 

APPEARING IN OTHER ROUTINES USED IN CONJUNCTION 

ROUTINE. 

STATEMENTS 280, 320, AND 375 MUST BE CHANGED TO DABS. 


INITIALIZATION 
I £R=0 


VARM 360 
VARM 370 
VARM 38C 
VARM 390 
VARM 400 
VARM 410 
VARM 420 
VARM 430 
VARM 440 
VARM 450 
VARM 460 
VARM 470 
VARM 480 
VARM 490 
VARM 500 
VARM 510 
VARM 520 
VARM 530 
VARM 540 
VARM 55C 
VARM 560 
VARM 570 
VARM 580 
VARM 590 
VARM 600 
VARM 610 
VARM 620 
VARM 630 
VARM 640 
VARM 650 
VARM 660 
VARM 670 
VARM 680 
VARM 690 
VARM 700 
VARM 710 
VARM 720 
VARM 730 
VARM 740 
VARM 750 
VARM 760 
VARM 761 
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0082 



FPS*0. 00116 

0083 



TVLT*G.O 

0084 



LL*K-1 

0085 



NV~1 

0086 



NC*0 

0087 



FN*M 

0086 



FFN*FN*FN 

0089 



CONS-O* 7071066 

0090 

C 



0091 

c 


CALCULATE ORIGINAL COMMUNALI TIES 

0092 

c 



0093 



DO 110 I = 1 » M 

0094 



HI I )=0.0 

0095 



DO 110 J=1 ,K 

0096 



L*M*( J-1I+I 

0097 


110 

H ( I )*H< I ) ♦ A { L ) *A < L 1 

0098 

c 



0099 

c 


CALCULATE NORMALIZED FACTOR MATRIX 

0100 

c 



0101 



DO 120 1*1 1 M 

0102 


115 

H( I i* SQRTlHi I ) ) 

0103 



DO 120 J*i,K 

0104 



L*M*( J-1KI 

0105 


120 

AtL»*AIL)/H( I ) 

0106 



GO TO 132 

0107 

c 



0108 

c 


CALCULATE VARIANCE FOR FACTOR MATRIX 

0109 

c 



0110 


130 

NV=NV*1 

0111 



TVLT=TV(NV-1> 

0112 


132 

TV ( NV ) =0*0 

0113 



DO 150 J*1 » K 

0114 



AA=0.0 

0115 



88*0*0 

0116 



LB*M* ( J- 1 ) 

0117 



00 140 1*1, M 

0118 



L*LB+I 

0119 



CC*A ( L ) *A ( L ) 

0120 



AA*AA*CC 

0121 


140 

BB*BB+CC*CC 

0122 


150 

TV C NV ) =TV ( NV ) ♦ ( Ftf^BB— AA*AA ) /FFN 

0123 



IF! NV-*5 1 >160,155,155 


VARM 770 
VARM 780 
VARM 790 
VARM 800 
VARM 810 
VARM 820 
VARM 830 
VARM 840 
VARM 850 
VARM 860 
VARM 870 
VARM 880 
- VARM 890 
VARM 900 
VARM 910 
VARM 920 
VARM 930 
VARM 940 
VARM 950 
VARM 960 
VARM 970 
VARM 980 
VARM 990 
VARM1000 
VARM1010 
VARM1020 
VARM1030 
VARM1040 
VARM1050 
VARM1060 
VARM1070 
VARM 10 BO 
VARM1090 
VARM 1100 
VARM 11 10 
VARM1120 
VARM1130 
VARM 11 40 
VARM1150 
VARM 1160 
VARM1170 
VARM 11 80 
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0124 

0125 

0126 

0127 

0128 
0125 

0130 

0131 

0132 

0133 

0134 

0135 

0136 

0137 

0138 

0139 

0140 

0141 

0142 

0143 

0144 

0145 

0146 

0147 

0148 

0149 

0150 

0151 

0152 

0153 

0154 

0155 

0156 

0157 

0158 

0159 

0160 
0161 
0162 

0163 

0164 

0165 


155 1ER«1 

GO TO 430 
C 

C PERFORM CONVERGENCE TEST 

C 160 ! F ( ( TV ( NV I -TVLT ) - ( 1 • E-7 ) I 170, 170, 190 


170 NC*NC+l 

I F ( NC-3 ) 190, 190, 430 

C 

C ROTATION OF TWO FACTORS CONTINUES UP TO 
C THE STATEMENT 120. 

c 

190 00 420 J*1 , LL 

Li=M*( J-l) 

II*J+1 

C 

c CALCULATE MUM AND DEN 

C 

DO 420 Kl=I I »K 
L2*M* ( K 1-1 ) 

AA=0.0 

BB=0.0 

CC*0.0 

DD-Q.O 

DO 230 1*1, M 

L3*L1+I 

L4*L2+I 

U= ( A ( L 3 ) ♦ A ( L4 ) ) * C A ( L3 1 *A C L4) 1 
T*AI L3 ) *A(L4 ) 

T»T ♦ T 

CC=CC* ( U*T ) *( U-T ) 

DD*DD*2.0*U*T 
AA*AA*U 
230 BB*BB*T 

T*DD-2.0*AA*BB/FN 
B«GC-<AA*AA-BB*BB) /FN 
C 

C COMPARISON OF NUM AND DEN 

C 

IFIT-B) 280, 240, 320 
240 IFHT»B)-EPS) 420, 250, 250 

C 



VARM1181 
VARM1182 
VARM 1190 
VARM1200 
VARM1210 
VARM1220 

VARM1230 

VARM1240 

VARM1250 

VARM 1260 

VARM1270 

VARM1280 

VARM1290 

VARM1300 

VARM 13 10 

VARM1320 

VARM1330 

VARM1340 

VARM1350 

VARM1360 

VARM1370 

VARM l 380 

VARM 1390 

VARM1400 

VARM 14 10 

VARM 1420 

VARM1430 

VARM1440 

VARM 1450 

VARM1460 

VARM1470 

VARM1480 

VARM 1490 

VARM1500 

VARM1510 

VARM1520 

VARM 15 30 

VARM1540 

VARM1550 

VARM1560 

VARM 1570 

VARM 1580 




-$* v'zciV.Kix tft sm 




1 • 
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0166 

C 

MUM ♦ DEN IS GREATER THAN OR EQUAL 

0167 

C 

TOLERANCE FACTOR 

0168 

C 


0169 

250 

C0S4T=C0NS 

0170 


SIN4T=C0N$ 

on i 


GQ TO 350 

0172 

C 


0173 

c 

NUM IS LESS THAN DEN 

02 74 ■ 

c 


0175 

280 

TAN4T* A8S ( T ) / ABS(B) 

0176 


IFITAN4T-EPS) 300, 290, 290 

0177 

290 

C0S4T=1.G/ SQRT ( 1 • 0*TAN4T*TAN4T ) 

0178 


SIN4T*TAN4T*C0S4T 

0179 


GO TO 350 

0180 

300 

IF (61 310, 420, 420 

0181 

31C 

SINP*CONS 

0182 


COSP=CONS 

0183 


GO TO 400 

0184 

C 


0185 

c 

NUM IS GREATER THAN DEN 

0186 

c 


0187 

32C 

CTN4T* ABSIT/B) 

0188 


IFICTN4T-EPS) 340, 330, 330 

0189 

3 30 

SIN4T=1.0/ SQRTU.0+CTN4T*CTN4T) 

0190 


C0S4T=CTN4T*SIN4T 

0191 


GO TO 350 

0192 

340 

COS4T*0.0 

0193 


SIN4T=1.0 

0194 

C 


0195 

c 

DETERMINE COS THETA AND SIN THETA 

0196 

c 


0197 

350 

C0S2T* SQRT( ( 1.0+C0S4T)/2.0) 

0198 


SIN2T=SIN4T/(2.0*C0S2T) 

0199 

355 

COST* SQRTt ( 1 •0+C0S2T) 72* 0 ) 

0200 


SINT*SIN2T/( 2. ONCOST) 

0201 

C 


0202 

C 

DETERMINE COS PHI AND SIN PHI 

0203 

c; 


0204 


IF ( B I 370, 370, 360 

0205 

360 

COSP*COST 

0206 


SINP=SINT 

0207 


GO TO 380 


VARM1590 

VARM1600 

V ARM 16 10 

VARM1620 

VARM1630 

VARM1640 

VARM1650 

VARM166D 

VARM 1670 

VARM1680 

VARM 1690 

VARM1700 

VARM1710 

VARM1720 

VARM 1730 

VARMI740 

VARM 1750 

VARM1760 

VARM1770 

VARM1780 

VARM 1790 

VARM1800 

VARM 1810 

VARM l 820 

VARM1830 

VARM1840 

VARM1850 

VARM1860 

VARM 1870 

VARM1880 

VARM1890 

VARM 1900 

VARM 19 10 

VARM1920 

VARM1930 

VARM1940 

VARM1950 

VARM1960 

VARM1970 

VARM 1980 

VARM1990 

VARM2000 
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0208 

0209 

0210 
0211 
0212 
0213 
021 4 

0215 

0216 

0217 

0218 

0219 

0220 
0221 
0222 

0223 

0224 

0225 

0226 

0227 

0228 

0229 

0230 

0231 

0232 

0233 

0234 

0235 

0236 

0237 

0238 

0239 

0240 

0241 

0242 

0243 


370 COSP=CONS*COST+GONS*SINT 

375 SIMP* ABS(C0NS*C0ST-C0NS*SINT> 

380 IF(T) 390, 390, 400 

390 SINP=-SINP 

C 

C PERFORM ROTATION 

C 

400 00 410 1*1, M 

13*11+1 
L4*L2+I 

AA*A(L3)*C0SP«-A<L4)*SINP 
A(L4)*-A(L3>*SINP+A(L4)*CQSP 
410 A(L3>*AA 

420 CONTINUE 

GO TO 130 
C 

C OENORMALUE VARI MAX LOADINGS 

C 

43C 00 440 |*l, M 

00 440 J*i,K 
L=M* ( J- 1 ) ♦ I 

44C AlL)=A(U*H(n 

C 

C CHECK ON COMMUNAL I TIES 

C 

NC=NV-1 
DO 450 1*1, M 
450 H(I)=H(II*HII) 

DO 470 I=1»M 
F(I)*0.0 
DO 460 J*1,K 
L*M*(J-i)+I 

460 F(n*Fin + A(L)*A(LI 

470 D( I l*H( I)-F1 I ) 

RETURN 

END 






VARM2G10 

VARM2020 

VARM2C30 

VARM2040 

VARM2050 

VARM2C60 

VARM2070 

VARM2080 

VARM2090 

VARM21G0 

VARM2110 

VARM2120 

VARM2130 

VARM2140 

VARM2150 

VARM2160 

VARM2170 

VARM2180 

VARM2190 

VARM2200 

VARM2210 

VARM2220 

VARM2230 

VARM2240 

VARM2250 

VARM2260 

VARM2270 

VARM2280 

VARM2290 

VARM2300 

VARM2310 

VARM2320 

VARM2330 

VARM2340 

VARM2350 

VARM2360 
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OUTPUT OF SUBR CO 
MODE » 2 MTYPE » 2 NVAR « 3 NLOC * 10 ISIOE » 


DATA MATRIX 

I 2 

1 7.00000 A. 00000 

2 A. 00000 1.00000 

3 3. 00000 8.00000 

VECTOR LENGTH 
* par l 8.60232 

IM ~ 2 9.00000 

. S 3 8.36660 

p4 A 10.0A988 

2? 5 11.74734 

£ 6 11.57584 

7 6.55744 

£0 8 13.03840 

9 9.48683 

5g 10 8.48528 

NORMEO DATA MATRIX 

1 2 

1 .81373 .44444 

2 .46499 .11111 

3 .34874 .88889 

COSINE MATRIX 

1 2 

1 1.00000 

2 .72332 1.00000 

3 .95870 ,88978 

4 4 .96007 .50857 

5 .95988 .87963 

6 .84355 .97905 

7 .99275 .79638 

8 .95399 .89479 

9 .98029 .84327 

10 .95899 .65473 


3 4 5 

6.00000 8.00000 8.00000 

3.00000 6.00000 5.00000 

5.00000 1.00000 7.00000 


3 

4 

5 

.71714 

.79603 

.68101 

.35857 

.59702 

•42563 

.59761 

.09950 

.59588 


3 

4 

5 

1.00000 

.84440 

1.00000 


.99710 

.85550 

1.00000 

.96024 

.66188 

.94863 

.98426 

.92562 

.98660 

.99920 

.83947 

.99891 

.99531 

.89153 

.99601 

.90150 

•91468 

.88283 


I 

I 


1 EF « .99500 


6 7 8 9 10 

7.00000 5.00000 9.00000 7.00000 8.GGC09 

2.00000 3.00000 5.00000 4.00000 2.00000 

9.00000 3.00000 8.00000 5.00000 2.C0000 


6 

7 

8 

9 

10 

60471 

.76249 

•69027 

.73786 

.94281 

17277 

.45750 

.38348 

.42164 

,23570 

77748 

.45750 

.61357 

.52705 

.23570 

6 

7 

8 

9 

10 


1.00000 

.89582 1.00000 

.96071 .98247 1.00000 

.92881 .99664 .99440 1.00000 

.79410 .93455 .88580 .91927 


1.00000 
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* 


INDEX 

EIGENVALUES 

CUM. SUM EV 


1 

9.15653 

9.15653 


2 

.74683 

9.90336 


3 

.09660 

9.99996 


4 

.00000 

9.99996 


5 

.00000 

9.99996 


6 

.00000 

9.99996 


7 

•00000 

9.99996 


8 

.00000 

9.99996 


9 

.00000 

9.99996 


10 

.00000 

9.99996 


EIGENVECTORS (COLUMNS) 



1 

2 

3 

1 

.32285 

-.24670 

-.03507 

2 

.28223 

.60175 

.04605 

3 

.32957 

.08500 

-.02579 

4 

•29088 

-.53799 

-.30689 

5 

.32873 

.06962 

-.26704 

6 

.31009 

.39336 

.20233 

7 

.32846 

-.11522 

-.15149 

8 

.32886 

.10142 

-.14579 

9 

.33015 

-.02084 

-.12839 

10 

.30598 

-.31112 

.85373 


aabva trvrasrao 


»■« rg m 


FACTORS < COLUMNS 1 

12 3 

.75247 -.39561 .52657 

.38153 -.38986 -.83810 

.53685 .83157 -.14242 


ij 

sj 
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ERMS- .00000 EMAX* 


.00001 




* 



ERROR OF REPRESENTATION 









1 2 

1 .00001 .ooooo 

2 .00000 ,00000 

3 .OOOOO ,00001 

3 

.00001 

-.00000 

.00001 

4 

.00001 

.00000 

•ooooo 

5 

•00001 

.00000 

.00001 

6 

.00001 

.ooooc 

.00000 

7 

•OOOOO 
.0000 1 
•ooooo 

8 

•ooooo 

.ccooo 

•ooooo 

9 

•ooooo 

.00001 

•ooooo 

10 

•ooooo 

•ooooo 

•ooooc 



VAR I MAX TRANSFORMATION 

I 2 

1 .72733 .68439 

2 -.68428 .72887 

3 -.05225 -.01854 


3 

.05078 

-.02227 

.99846 
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s§ 


O tz 


VSTAR MATRIX 
1 

1 *85702 

2 • 26457 

3 .67550 

4 .96332 

5 .68667 

6 .44658 

7 .79350 

8 .66618 

9 .74104 

10 .84356 


FSTAR MATRIX 

1 

1 .79049 

2 .58806 

3 -.17112 



2 

3 

51343 

.04348 

96326 

.04608 

73621 

.04101 

26530 

-.04018 

72618 

-.03369 

88880 

.10287 

60853 

.00568 

74578 

.00334 

67135 

.01130 

43279 

.31794 


2 3 

.21688 .57278 

.00749 -.80875 

.97616 -.13345 
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